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Abstract 


Quantitative  Feedback  Theory  developed  by  Dr.  Isaac 
Horowitz  of  the  University  of  Colorado  is  used  to  design 
the  control  laws  for  a  Short  Take  Off  and  Landi ng ( STOL ) 
aircraft.  Compensators  are  presented  for  two  longitudinal 
variaoles,  angle  of  attack  and  forward  velocity,  which  are 
controlled  via  the  use  of  five  separate  control  surfaces: 
canard,  stabilator,  ailerons,  upper  and  lower  thrust 
reversing  vanes.  The  final  design  must  exhibit  robust 
qualities  over  three  flight  conditions  despite  surface 
failures. 

The  state-space  matrix  representation  of  the  aircraft  is 
developed  from  pertubation  equations  using  linearized 
aerodynamic  data.  Transfer  functions  relating  servo  input 
signals  to  aircraft  outputs  are  obtained  from  the  state-space 
equations.  The  original  output  set  included  the  flight  path 
angle  and  velocity;  however,  the  non-minimum  phase 
cnaracteri sti cs  of  the  flight  path  angle  precluded  its  use  by 
the  type  of  Quantitative  Feedback  Theory  used  in  this  thesis 
since  unstaole  plants  can  arise.  Instead,  the  minimum  phase 
variaoles  angle  of  attack  and  velocity  are  controlled.  The 
ten  separate  transfer  functions  relating  the  two  output 
variables  to  the  five  input  commands  form  a  5  X  2  plant 
transfer  function  matrix.  These  separate  transfer  functions 


are  combined  using  a  weighting  vector  into  a  2X2  minimum 
phase  plant  matrix  for  each  flight  condition/failure 
combination.  Quantitative  Feedback  Theory  is  applied  to  the 
resulting  plants  to  yield  robust  control  . 

A  single  set  of  fixed  compensators  and  prefilters  are 
designed  to  handle  the  entire  plant  set,  consisting  of  three 
single-surface  failures  and  two  dual-surface  failures  at  each 
flight  condition.  For  these  failures  neither  Fault 
Detecti on/ Identi f icati on  ,  nor  scheduled  compensation,  is 
required.  Surfaces  are  assumed  locked  at  zero  degrees 
deflection  after  failure,  generating  no  net  moment  after 
failure.  Digital  simulations  have  shown  the  control  to  be 
robust  over  the  three  flight  conditions  and  surface  failures. 
Loop  bandwidths  for  the  velocity  and  angle  of  attack  loops 
are  35  and  12  rad/sec  respectively.  Control  surface  rates  and 
deflections  are  shown  to  saturate  only  for  the  double  failure 
cases . 

Quantitative  Feedback  Theory  effectively  controls  the 
aircraft  despite  large  uncertainty  due  to  flight  condition 
changes  and/or  control  surface  failures  without  identification 
Application  of  QFT  eliminates  the  use  of  identification  to 
acnieve  robustness  and  the  associated  false  alarm  and  missed 
detection  problems.  Efforts  to  expand  upon  the  base  of  flight 
control  design  using  this  method  are  recommended,  especially 
direct  design  in  the  discrete  domain.  Research  should  also 
continue  on  developing  a  computer-aided  design  program  to 
expedite  the  synthesis  of  controllers  using  QFT. 


ROBUST  FLIGHT  CONTROL  SYSTE'-I  FOR  A  STOL  AIRCRAFT 
DESIGN  USING  QUANTITATIVE  FEEDBACK  THEORY 


I .  Introduction 

Future  aircraft  may  incorporate  many  control  surfaces 
to  meet  stringent  performance  and  aerodynamic  efficiency 
requ i rements .  These  many  surfaces  increase  the  survivability 
of  tne  aircraft  by  providing  redundant  control  in  the  event 
of  surface  failures.  The  design  of  a  flight  control  system 
that  automatically  redi stributes  control  authority  among  the 
remaining  surfaces  in  the  event  of  failure  is  a  significant 
challenge.  A  control  law  sufficiently  robust  to  encompass  the 
dynamic  uncertainty,  as  well  as  surface  failures,  is  desired. 
Failure  identification  is  useful  but  there  is  a  trade  off 
between  the  reliability  of  the  failure  identification  and  the 
time  critical ity  of  the  failure.  Thus  a  robust  (non¬ 
identification)  design  provides  control  without  the  need  of 
reliable  identification. 

Quantitative  Feedback  Tneory(QFT)  developed  by  Dr.  Isaac 
Horowitz  promises  to  yield  roDUStness  witnoot  identification. 
Tnis  tneory  has  been  successfully  applied  to  a  number  of 

^  v 

difficult  problems[l ,2 ,3 ,4 j .  ~QFT  inherently  includes 
uncertainty  and  control  system  failures  within  the  design 
procedure.  One  a  priori  designs  for  acceptaole  system 
responses  with  uncertainty  and  failures,  raking  QFT  w e 1 1 


thesis  uses  QFT  in  the  design  of  a  flight  control  system  for 
a  future  Air  Force  Short-Take-Off-and-Landing(STOL) 
experimental  aircraft. The  aircraft  nas  additional  control 
surfaces  not  found  on  current  aircraft.  Hence  it  provides  an 
excellent  platform  for  demonstrating  reconf i gureabl e  flignt 
control  system  design. 

Tne  remainder  of  cnapter  one  provides  additional 
oacxground  information,  sets  the  scope  of  the  problem,  lists 
tne  assumptions  made,  and  outlines  the  approach  used. 

I . 1  Background 

The  Air  Force  Flight  Dynamics  Laboratory  is  currently 
investigating  the  development  of  reconf i gurabl e  flight 
control  systems.  The  primary  goals  are  better  reliability, 
maintainability,  survivability  ,  simplicity,  and  reduced  life 
cycle  costs[5]. 

a.  Reliability,  maintainability,  and  economy  will  be 
enhanced  by  the  use  of  simpler,  smaller,  and  less  redundant 
control  system  components.  Identical  components  simplify 

,  .id  i  nta  i  nance  problems  by  reducing  the  variety  of  parts 
required  in  inventory.  Training  and  maintainance  problems  are 
reduced  since  personnel  have  to  remember  less  aoout  each 
system  due  to  system  simplicity.  Larger  numoers  of  each 
particular  servo  could  reduce  costs  from  economies  of  scale. 

b.  The  existence  of  functionally  redundant  control 


surfaces  on  the  airframe  reduces  the  sensitivity  of  the 
control  system  to  the  loss  of  surfaces.  The  dependancy  of 
aircraft  flight  control  on  any  one  surface  is  reduced.  This 


leads  to  increased  survi vabi 1 i ty  in  the  event  of  damage, 
c.  Many  current  aircraft  have  large  control  surfaces  driven 
by  redundant,  expensive  actuators.  Less  complex  actuators 
promise  to  reduce  procurement  costs  and  maintenance  of  flight 
control  systems.  The  current  complex  actuators  driving 
critical  control  surfaces  must  withstand  multiple  failures 
and  still  function.  Such  actuators  are  expensive  and 
difficult  to  maintain.  It  is  estimated  that  simpler  control 
actuators  may  decrease  the  life  cycle  cost(LCC)  of  the  flight 
control  system(FCS)  by  30  percent[5j.  Flignt  control  surfaces 
utilized  in  tnis  thesis  include: 

a.  Canard 

b.  Stabilator 

c.  Ailerons 

d.  Thrust  reversing  vanes 

As  previously  noted,  QFT  is  inherently  suited  to  design 
flight  control  systems  for  reconf i gurabl e  aircraft.  In  QFT 
the  uncertainties  and  failures  are  considered  beforehand,  and 
the  design  acnieves  the  desired  response.  QFT  also  allows  the 
designer  to  see  clearly  tne  trade-offs  between  the  extent  of 
tne  resulting  plant  uncertainties,  the  narrowness  of  the 
performance  tolerances,  and  the  resulting  bandwidth(  the 
"cost  of  feedoack").  For  instance,  certain  flight  conditions 
and  failure  sets  may  lead  to  compensation  elements  with 
unrea 1 i s t i ca 11 y  large  bandwidths.  QFT  allows  the  designer  to 
identify  these  problems  at  the  beginning  of  the  design  cycle. 


QFT  employs  well  developed  frequency  domain  design  tecnniques 
that  are  simple  to  apply  and  give  the  designer  a  "feel”  for 
the  problem  and  its  trade-offs. 

1 . 2  Problem 

The  focus  of  this  thesis  is  applying  QFT  to  design  a 
reconfigurable  aircraft  control  system.  Performance 
tolerances  on  two  longitudnal  outputs  must  oe  satisfied  over 
a  range  of  flight  conditions  including  multiple  control 
surface  failures.  Possiole  operating  conditions  are: 

a.  All  surfaces  operating  normally. 

b.  Canard  failed,  fixed  at  trim  position. 

c.  Ailerons  failed,  fixed  at  trim  position. 

d.  Elevators  failed,  fixed  at  trim  position. 

e.  Thrust  vectoring  failed,  fixed  at  trim. 

f.  Combinations  of  the  above. 

1 . 3  Assumptions 

Tne  assumptions  made  in  this  thesis  include: 
a.  The  aircraft  equations  of  motion  can  be  linearized 
about  an  equ i 1 i or i urn ( tr im )  position  so  that 
small  perturoation  models  are  valid, 
o.  Commanas  and  aircraft  responses  do  not  invalidate  tne 
assumed  linear  model. 

c.  Mass  remains  constant  during  the  command  sequence. 

d.  Tne  control  surfaces  failures  are  assumed 


symmetrical.  This  eliminates  any  cross  coupling  into 


lateral  modes.  This  thesis  does  not  investigate 
-cross  coupling  because  there  are  already  7  control 
surfaces  in  the  model(  canard,  ailerons,  stabilator, 
top  and  bottom  vane  pair).  Without  this  constraint  the 
problem  would  be  beyond  the  scope  of  a  single  thesis. 
Arnold[l]  has  applied  QFT  for  the  design  of  a 
reconf i gurabl e  FCS  with  cross  coupling, 
f.  Actuator  models  given  by  McDonnell  Douglas  are 

assumed  for  tne  QFT  design.  They  are  third  order  for 
tne  canards,  ailerons,  and  staoilators,  and  first 
order  for  tne  tnrust  reversing  vanes. 

These  assumptions,  as  well  as  others  introduced  later, 
are  explained  in  the  body  of  the  tnesis.  An  additional 
assumption  is  that  the  reader  is  familiar  with  state-space 
representai on  of  dynamic  systems,  frequency  response 
characteristics,  transfer  functions,  matrix  representations, 
and  matri x  al gebra . 

I . 4  Design  Outline 

The  folowing  steps  outline  the  control  system  design 
process  used  in  this  thesis: 

a.  The  control  structure  and  variables  are  defined, 
o.  Plant  matrices  are  derived  from  aerodynamic 
oata.  Tnese  matrices  are  deriveu  at  eacn  flight 
condition  and  failure  case. 


c.  Transfer  functions  relating  each  of  the  five 
inputs  to  each  of  the  two  outputs  are  derived 


from  the  plant  matrices. 

d.  The  ten  individual  transfer  functions  P  are 

i  j 

combined  into  an  equivalent  two  input  and  two 

output  plant  transfer  function  matrix, 

P  =  [p  ],  by  combining  the  individual  surfaces 
i  J 

through  the  use  of  a  weighting  vector  L  .  This 
vector  is  originally  chosen  to  divide  the 
control  authority  among  the  surfaces  in  a  set 
percentage.  Tne  resulting  non-minimum  pnase 
plants  can  lead  to  instability  using  QFT,  so 
the  A  vector  is  moaified  such  that  the 
resulting  P  are  minimun  phase. 

e.  Reconf i gurabl e  terms  u.  are  developed  to  make 

i  j 

the  plant  matrix  diagonally  dominant.  Diagonally 
dominant  systems  exhibit  less  cross  coupling 
between  channels  which  result  in  decreased 
loop  transmi ssi om  bandwidths.  The  u.  also 

ij 

increases  robustness  under  failures  by  feeding 
command  authority  normally  reserved  for  one 
control  channel  into  the  other  in  case  of  a 
failure. 

f.  Tne  resulting  plant  matrices  are  inverted, 

-1 

£  =lP*  j.tnen  eacn  individual  element  inverted, 

i  j 

to  for.n  tne  ^  matrix,  wnere 


£=iq  J=L  1/P*  ]•  (1.4-1) 

1  j  i  J 


g.  Aircraft  time  domain  performance  specifications 


are  converted  into  equivalent  frequency  domain 
speci f i ca  t i ons  . 

h.  QFT  converts  the  multiple  input-multiple 

output(MIMO)  synthesis  problem  into  equivalent 
SISO  synthesis  problems.  These  are  solved, 
resulting  in  loop  compensation  for  the  system, 
g.  The  design  is  simulated  and  results  analyzed 
over  the  design  range  of  flight  conditions  and 
failures.  Control  surface  deflections  and  rates 
are  examined  to  ensure  acceptaDle  limits. 

I . 3  Tnesis  Presentation 

Tnis  tnesis  is  organized  into  six  chapters.  Chapter  I  is 
the  introduction.  Cnapter  II  develops  the  control  structure 
and  develops  the  equivalent  plant  matrices  from  the  aircraft 
equations  of  motion.  In  Chapter  III,  QFT  for  SISO  systems  is 
explained,  as  well  as  the  development  of  the  equivalent  SISO 
systems  from  the  original  MIMO  system  and  reconfiguration 
theory.  Chapter  IV  applies  QFT  to  the  problem  at  hand, 
developing  compensators  and  prefilters  for  the  control 
system.  The  system  is  simulated  in  Chapter  V  to  determine  now 
well  the  design  satisfies  the  desired  performance 
specifications.  Chapter  VI  contains  the  results,  conclusions, 
and  recommendations. 


Derivation  of  Aircraft  Transfer  Functions  and 
Control  Structure  from  Aerodynamic  Data 


1 1 . 1  Introduction 

This  chapter  develops  the  aircraft  transfer  functions 
needed  in  Q F T  design.  The  aerodynamic  data  is  used  to  develop 
aircraft  equations  of  motion  around  an  equilibrium  point. 
These  equations  of  motion  are  then  arranged  in  state  equation 
format,  and  are  then  transformed  to  the  complex  frequency 
domain  to  form  transfer  functions  relating  each  input  to  each 
output.  Tne  individual  transfer  functions  are  then  grouped  to 
form  an  equivalent  two  input-two  output  plant  matrix.  In  the 
process  weighting  factors  are  introduced.  These  are  chosen  to 
maximize  tne  ratios  of  tne  diagonal  elements  of  the  plant 
matrix  to  its  off-diagonal  el ements ( ensure  diagonal 
dominance) . 

1 1 . 2  Ai rcraf t 

The  aircraft  that  forms  the  basis  for  this  thesis  effort 
is  the  Short-Takeoff-and-Landi ng( STOL  )  F - 1 5  currently  under 
development  by  McDonnell  Douglas[6j.  Figure  II. 2-1  is  a 
diagram  of  the  proposed  aircraft.  This  aircraft,  derived  from 
the  F - 1 5  Eagle  air  superiority  fighter,  incorporates  two 
dimensional  thrust  vectoring  vanes,  thrust  reversing  vanes, 
and  a  forward  canard.  Tne  original  purpose  for  the  extra 
control  surfaces  is  to  provide  extra  force  and  moment 


Fig. 1 1.2-1:  F - 1 5  STOL  Ai  rcraft 


production  for  STOL  operation.  However,  these  same  extra 
control  surfaces  promise  to  provide  redundant  control  under 
failure  conditions.  Basic  specifications  for  the  STOL  F - 1 5 
aircraft  are  in  Table  II. 2-1.  Appendix  A  contains  the 
a*rodynamic  data  us,ed  to  develooe  the  state  space  models. 


TABLE  II. 2-1 


STOL  F - 1 5  Aircraft  Data 


Aircraft  Parameters 

Wing  Mean 

Aerodynami c 

Cord  c 

15.8  ft 

Wing  Reference  Area 

S 

608  sq-ft 

Wing  Span 

b 

42.7  ft 

Weight 

w 

33576  lbs 

I  XX 

23634  slug-sq-ft 

iyy 

181837  slug-sq-ft 

Izz 

199674  si ug-sq-ft 

I  xz 

-3086  slug-sq-ft 

I  I  .  3  Aerodynamic  Model 

The  equation  of  motion  for  the  aircraft  are  written  in 
state  matix  format,  with  the  individual  elements  derived  from 
the  aerodynamic  data.  This  format  is: 

x.  =  A_x  +  Bu_  (2.3-la) 

^  =  £x  +  D£  (2.  3- lb) 


where 


x_  is  the  state  vector 

A  is  the  plant  matrix 

B_  is  the  forcing  function  matrix 

^  is  the  output  vector 

u  is  the  forcing  function  vector 


£  is  the  output  matrix 
D  is  the  feedforward  matrix 


The  elements  of  the  A  and  £  matrices  are  derived  from 
the  aerodynamic  data  in  Appendix  A.  The  £  matrix  is  chosen  to 
output  the  desired  quantities  while  D  is  equal  to  0,  the  null 
matrix.  The  data  provided  by  McDonnell  Douglas,  which  is 
defined  in  terms  of  the  body  axis,  is  converted  to  the 
s t a o i 1 i t y  axis.  Tnis  conversion  is  done  to  facilitate 
analysis  of  tne  state  space  model  by  other  thesis  students  using 
tnis  aircaft  ooth  in  the  longitudnal  and  lateral  directions. 
Since  this  thesis  concerns  only  the  longitudnal  equations  QFT 
will  give  a  stable  design  in  either  axis  system. 

This  thesis  investigates  only  the  longitudinal  equations 
of  motion  for  the  aircraft.  The  model  used  consists  of  four 
states:  the  f orwa  rd  velocity,  v,  the  pitch  rate,  q,  the  angle 
of  attack,  <x  ,  and  the  pitch  angle, 0.  Initially  the  outputs 
chosen  were  the  forward  velocity,  and  the  flight  path  angle, Y. 
There  are  two  longitudinal  variables  one  would  like  to 
control,  especially  while  landing.  The  state  and  output 
equations  are  then: 


•  t  - 

- 

V 

_  _ 

V 

q 

q 

= 

A 

<* 

i  J 

0 

- 

© 

B 

i  J 


6.  (2-3-2) 


1.0 

0.0 

o 

o 

o 

o 

O 

o 

0.0 

-1.0 

1.0 

(2.3-3) 


The  components  of  the  A  and  j3  matrices  are  stability 
derivatives  generated  by  linearizing  the  actual  non-linear 
aircraft  equations  of  motion  about  an  equilibrium  point. 
Entries  in  C  are  chosen  to  give  the  desired  outputs. 

Of  the  various  control  surfaces  on  the  STOL  F - 1 5  the 
following  are  used  in  this  thesis: 

a.  Canard 

b .  Stabi 1 ator 

c .  Ailerons 

d.  Top  reversing  vanes 

e.  Bottom  reversing  vanes 

The  flaps  are  not  used  since  they  only  deflect  downward.  If 
they  are  used  it  is  very  likely  that  upward  flap  deflection 
would  at  times  be  demanded.  Since  this  is  impossible,  a 
nonlinearity  would  result  compounding  the  control  problem.  In 
this  initial  investigation  it  is  desired  to  avoid  such 
problems,  so  the  flaps  are  excluded  from  the  reconf i gurabl e 
controller.  However,  in  the  actual  aircraft  design  the  flaps 
should  be  included  as  another  possible  useful  surface 
in  case  of  failure.  The  2-0  nozzles  are  not  included  since 
the  stability  derivatives  required  in  the  matrix  are  not 
Known  at  the  time  of  this  work.  Appendix  B  details  the  state 
space  model  and  the  data  involved. 

Using  the  Compu te r- A i ded-Des i gn ( CAD )  package  T0TAL[7j, 
the  transfer  functions  relating  each  input  to  each  output  are 
calculated.  These  are  the  original  plant  transfer  functions 
that  are  combined  to  form  the  equivalent  symmetrical  plant 


ma tr i x ( i npu ts  =  outputs)  required  for  QFT.  For  the  design  in 
this  thesis  the  plant  matrix  is  2  X  2.  The  signal  flow  graph 
relating  the  surface  deflections  to  the  aircraft  responses  is 
shown  in  Figure  II. 3-1.  Assuming  that  the  transfer  functions 
do  exist  from  each  control  surface  to  each  output,  the 
transfer  function  matrix  representation  of  this  relationship  is: 


a 


a 


a 


tv 


is  the  canard  deflection 
is  the  stabilator  deflection 
is  the  aileron  deflection 
is  the  top  vane  deflection 
is  the  bottom  vane  deflection 
is  output  1,  forward  velocity 
is  output  2,  flight  path  angle 
As  is  seen  later  on  the  output  vector  had  to  be  changed  to 
forward  velocity  and  angle  of  attack,  but  still  with  the  same 
matri x  structure. 


bv 


v 

y 


The  plant  transfer  function,  P  ,  relates  the  J_t h  output 

i  j 

to  the  jth  input.  It  is  assumed  that  each  control  surface  is 


driven  by  a  servomechanism  with  a  transfer  function  of  M 


such  that 


P  i  = 
i  j  i 


A  A 
M  P  6 

i  i  j  i 


(2.3-5) 


where 


6  is  the  command  input  to  the  servo.  The  input- 
i 

output  relationship  now  becomes 
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(2-3-6) 


Since  two  outputs  are  being  controlled  by  two  inputs  through 
five  separate  control  surfaces,  there  are  many  possible  ways 
available  for  dividing  the  control  effort  between  the 
surfaces.  From  preliminary  discussions  with  personnel  from 
the  Air  Force  Flight  Dynamics  Laboratory[8]  it  was  decided 


that  the  control  of  the  flight  path  angle  would  come 

primarily  from  the  canard,  stabilator,  and  ailerons  while  the 

reversing  vanes(top  and  bottom)  would  provide  the  f orwa rd 

velocity  control.  A  variable  A  is  used  to  determine  how  much 

i 

control  authority  is  granted  to  the  j_th  surface  to  control 
its  respective  output.  This  term  does  not  have  to  be  a 
constant,  it  could  contain  terms  to  compensate  for  the 
frequency  response  of  the  surface  it  weights.  Next,  the  a 
coefficients  are  another  variable  transfer  function  in  the 
forward  path  that  combine  with  the  ft  terms  to  provide  the 
cross  coupling  required  for  robust  control  under  surface 
failures  and  disturoance  rejection  for  the  aircraft  under  no¬ 
fail  conditions.  These  terms  can  also  be  functions  of 
frequency  to  correct  for  differences  between  control 
surfaces.  Figure  II. 3-2  is  the  configuration  of  the  plant  and 
the  control  structure  previously  assumed.  The  equations 
relating  the  equivalent  2  X  2  plant  to  the  individual  control 
surface  transfer  functions  are 

P'  =  a  {  A  P  +  A  P  +  A  P  +  A  P  +  A  P  ) 

11  1  11  11  12  12  13  13  14  14  15  15 

+  ft  (AP  +  A  P  +  A  P  +  A  P  +  A  P  )  (2 . 3-7a ) 

21  21  11  22  12  23  13  24  14  25  15 

P*  =  a  (  A  P  +  A  P  +  A  P  +  A  P  +  A  P  ) 

12  2  21  11  22  12  23  13  24  14  25  15 

+  M  (AP  +  A  P  +  A  P  +  A  P  +  A  P  )  (2.3-70) 

12  11  11  12  12  13  13  14  14  15  15 

P*=  <j(AP  +  A  P  +  A  P  +  A  P  +  A  P  ) 

21  1  11  21  12  22  13  23  14  24  15  25 


%T  ^ 

t  -• 


U 


V/- 


+  H  {A?  +  A  P  +  A  P  +  A  P  +  A  P  )  ( 2 .  3-  7c ) 

21  21  21  22  22  23  23  24  24  25  25 


F i q . 1 1 . 3-2 ;  General  Plant  Structure 


P*  =  <7  (  A  P  +  A  P  +  A  P  +  A  P  +  A  P  ) 

22  2  21  21  22  22  23  23  24  24  25  25 

+  //  (AP  +  A  P  +AP  +  A  P  +  A  P  )  (2.3-7d) 

12  11  21  12  22  13  23  14  24  15  25 

This  is  quite  a  complicated  system,  especially  if  A  ,  a  , 

ij  1 

and  n  are  actually  A  ( j  w ) ,  a  ( j  w  )  ,  and  H  ( j  w ) .  The  best 
ij  ij  i  ij 

choice  for  the  A'  s ,  a  1 s ,  and  ^'s  constitute  an 

optimization  problem  beyond  the  scope  of  this  masters  thesis 

which  is  primarily  dedicated  to  the  reconfigurability 

synthesis  problem.  In  order  to  simplify  these  relations  to 

provide  a  better  "feel"  to  the  designer  on  how  best  to  pick 

and  choose  the  respective  weighting  terms,  the  a  terms  are 

set  equal  to  unity,  the  A's  are  reduced  to  five  in  all  - 

set  in  the  main  forward  path  for  the  dominant  variable  each 

set  of  surfaces  is  primarily  controlling.  A  ’  A  ’ 


and  A  are  in  the  forward  path  for  the  flight  path  angle 
3 

channel  while  A  ,and  A  are  in  the  velocity  channel.  This 

4  5 

simplified  control  scheme  is  shown  in  Figure  II. 3-3.  Some 
freedom  is  of  course  lost  by  doing  this,  but  the  design 
process  and  assignment  of  the  intermediate  plant  transfer 
functi ons (  see  below)  are  facilitated. 

The  equivalent  plant  transfer  functions  now  become 


P'  =  AP  +  AP  +  AP  +  AP  +  AP  )  ( 2 . 3-8a ) 

11  1  11  2  12  3  13  21  4  14  5  15 

P'=  AP  +  A  P  +  //  (  A  P  +  A  P  +  A  P  )  (  2 . 3  -  8b  ) 

12  4  14  5  15  12  1  11  2  12  3  13 

P'  =  AP  +  A  P  +  AP  +  if  (  A  P  +  AP  )  (2.3  —  8c) 

21  1  21  2  22  3  23  21  4  24  5  25 

P'  =  AP  +  A  p  +  A  P  +  A  P  +  A  P  )  ( 2 . 3-8d ) 

22  4  24  5  25  12  1  21  2  22  3  23 

Intermediate  plant  transfer  functions  in  the  parenthesis 
are  defined  as  follows: 


AP  +  A  P  +  A  P 
1  11  2  12  3  13 


( 2 . 3-9a  ) 


P  =  A  P  +  A  P 
12  4  14  5  15 


( 2 . 3-9b  ) 


P  =  AP  +  A  P  +  A  P 
21  1  21  2  22  3  23 


(2.3-9c) 


P  =  AP  +  AP 
22  4  24  5  25 


(  2  .  3  -  9  d  ) 


These  equations  are  substituted  into  equations  2.3-8a-d  to 


yield: 


P  +  n  P 
11  21  12 


(2.3-10a) 


P '  =  P  +  fi  P 

12  12  12  11 


(2.3-10b) 


P  +  n  P 


(2.3-lOc) 


F i g . 1 1 . 3-3 :  Simplified  Plant  Structure 


P'  =  P  +  H  P  (2.3-10d) 

22  22  12  21 


The  entire  control  structure  with  the  feedback  loops  in  place 

is  in  Figure  II. 3-4.  The  f  and  g  elements  are  p re  filters 

i  i  i 

and  loop  compensators  respectively  to  be  designed  by  QFT. 

The  QFT  design  process  is  discussed  in  Chapter  3. 


III.  QFT  Theory 
III.l  I ntroducti on 

Most  control  system  design  techniques  currently  used 
design  for  a  specific  set  of  outputs  given  a  particular 
system  model.  The  designer  hopes  that  in  the  process  enough 
robustness  has  been  built  in  to  the  design  that  it  will 
tolerate  some  plant  variations,  failures,  noise,  and  other 
types  of  uncertainty.  If  the  first  design  is  not 
sati sfactory,  the  designer  goes  back  and  increases  the  gain, 
moves  a  pole,  changes  the  weighting  matrix,  or  possibly  a 
combination  of  these  things  in  order  to  arrive  at  a 
satisfactory  design.  Quantitative  Feedback  Theory(QFT), 
developed  by  Dr.  Isaac  Horowitz,  allows  the  designer  to  put 
in  the  uncertainties  £  priori  ,  and  results  in  a  design 
guaranteed  to  meet  the  system  specifications.  QFT  has  many 
advantages,  including: 

a.  Uncertainties  are  included  at  the  outset,  resulting 
in  a  robust  design. 

b.  Frequency  domain  design  techniques,  which  are  well 
developed  and  understood,  are  used. 

c.  Complicated  multiple-input  multiple-output 
systems(MIMO)  are  reduced  to  a  set  of  single-input 
single-output  sys terns ( S I  SO )  which  simplifies  the 
design  process. 

d.  The  design  process  is  transparent,  that  is,  it 


allows  the  designer  to  see  exactly  what  "the  cost  of 


feedback" ( requ i red  loo?  transmission  bandwidth)  is 
for  the  range  of  uncertainties  chosen.  A  problem 
which  requires  unrealistic  compensation  will  be 
apparent  from  the  outset  of  the  design. 

This  chapter  is  broken  into  four  separate  parts.  The 
introduction,  QFT  SISO  system  design  techniqe,  MIMO  reduction 
into  equivalent  SISO  systems,  and  the  application  of 
reconf i guraDl i ty .  For  a  much  more  indeptn  treatment  the 
reader  is  urged  to  study  the  references!.  9,  10,  11,  12  ,13,114  J. 


III. 2  QFT  SISO  Design  Technique 
a .  Overv i ew 

In  order  to  develop  QFT,  a  specific  SISO  structure  used 

in  this  thesis  is  introduced.  Next  the  system  time  domain 

specifications  are  translated  into  frequency  domain 

specifications.  The  plant  model  and  system  specifications 

define  the  bounds  of  the  system  at  every  specific  frequency 

value  of  interest.  Uncertainties  in  the  plant  P  at  each 

frequency  of  interest  are  represented  by  an  area  on  the 

Nichol's  chart  known  as  the  plant  template.  Usiny  the  plant 

templates  and  system  specifications,  bounds  on  the  nominal 

loop  transmission  L  are  derived.  An  L  is  then  designed 

0  0 

that  satisfies  the  bounds.  From  the  L  a  loop  compensator  G 

0 

is  determi ned  ,  i . e . 


G 


t  /  P 


(  3 . 2  .  a  - 1 ) 


where  P  is  the  nominal  plant  transfer  function  chosen, 
o 

A  prefiTter,  F,  is  then  designed  to  give  the  required 
tracking  specifications  over  the  range  of  plant 
uncertainties. 


b .  S I S 0  System  Definition 

The  system  type  that  this  thesis  is  concerned  with  is 
the  two-degree  of  freedom  structure  in  Figure  III. 2-1.  In 
this  thesis  tne  disturbance  is  assumed  to  appear  at  the  input 
of  the  plant  rather  than  at  the  output.  It  is  also  assumed 
tnat  the  plant  P  has  a  known  range  of  uncertainties,  that  the 
tracking  input  r(t),  disturbance  input  d(t),  and  the  output 
y(t)  are  given.  Here  y(t)  is  actually  a  member  of  the  set  of 
acceptable  output  responces  Y ( t ) .  The  signal  x(t)  is  the 
plant  input  such  that  y(t)  =  x  ( t )  *  p  ( t ) . 

It  is  assumed  that  r(t)  and  y(t)  are  measureable,  the 
output  y(t)  is  fed  back,  and  that  Laplace  transforms  exist 
for  every  signal  and  system  element  present.  For  the  two- 
degree-o f- freedom  structure  the  designer  has  to  design  two 
compensation  elements,  the  prefilter  F,  and  the  loop 
compensator  G.  In  terms  of  the  loop  transmission  L  =  P  G  , 
four  separate  transfer  functions  can  be  defined  to  describe 
the  behavior  of  the  system.  The  responses  due  to  the  command 


input  and  disturbance  input  are: 


y  =  ^F  G  P  /  ( 1  +  GP)  Jr  =  j^FL/(l  +  L)J  r 
y,  =  Fp/(1  +  GP) 1  d  =  fp/(i  +  L )  1  d 


(3.2-1) 


(3.2-2) 


F i g . 1 1 1 . 2 -  1 :  Two  Oegree  of  Freedom  Structure 

The  transfer  functions  relating  the  plant  input  to  command 
and  disturbance  inputs  are  respectively: 

I  =  X(s)/R(s)  =  F  G / ( 1  +  L)  (3.2-3) 

r 

I  =  X ( s )/D( s )  =  1/(1  +  L)  (3.2-4) 

d 

This  thesis  uses  only  the  output  equations  to  develop  the 
loop  compensations.  The  internal  variables  are  ignored  in 
this  respect.  This  is  not  to  say  that  the  internal  variables 
are  not  important.  Some,  such  as  control  surface  deflections 
or  rates,  can  be  critical.  Quantitative  methods  exist  to  deal 
with  the  internal  variables,  but  they  are  beyond  the  scope  of 
tnis  thesis[  15  J.  Internal  va ri abl es ( such  as  the  control 
surface  deflections)  are  checked  at  the  end  of  the  design  to 


make  sure  that  limits  have  not  been  exceeded.  If  they  have, 
the  designer  must  revaluate  the  problem,  make  appropriate 
changes,  and  begin  again. 


c .  System  Specifications 

System  specifications  can  be  given  in  the  time  domain, 

such  as  rise  time,  settling  time,  oversnoot,  peak  time,  and 

final  value,  or  they  can  be  expressed  in  the  frequency 

domain,  such  as  bandwidth,  OC  gain,  resonant  frequency,  phase 

margin,  and  gain  margin.  Which  domain  they  are  given  in  does 

not  matter  since  they  can  always  be  transposed  to  the  other. 

Typical  time  response  curves  are  shown  in  Figure  III. 2-2.  The 

faster  rising,  underdamped  response  has  a  higher  oand width 

and  is  known  as  the  upper  bound  T  .  The  overdamped  response 

U 

has  a  lower  bandwidth  so  it  is  known  as  the  lower  bound  T 

L 

The  disturbance,  d  ,  is  shown  at  its  maximum  allowable  value. 
These  same  specifications,  now  expressed  in  the  frequency 
domain,  are  in  Figure  III. 2-3  .  The  upper  and  lower  bounds 
are  evident  as  is  the  "typical"  disturbance  rejection  shape 
of  d(ju>)  .  However  the  specifications  are  displayed,  the 
range  between  the  bounds  reflects  the  uncertainty  in  the 
output  which  leads  to  an  acceptable  response. 

The  time  response  of  a  desired  output  variable  of  a 
control  system  is  very  rarely  a  single  function,  ratner  a  set 
of  functions  is  specified  tnat  meet  tne  given  response 
specifications. 


Tracking  and  Disturbance  Bounds 
in  the  Frequency  Domain 


d.  Derivation  of  Plant  Templates 


In  QFT  contraints  on  the  free  functions  F,  6  are  found 

so  that  the  specifications  are  satisfied.  The  great  advantage 

of  frequency  response  is  poi ntwi se  synthesi s-the  constraints 

are  found  pointwise  at  each  u  value.  The  set  of  all 

possible  plants  P^(jw)  at  a  specified  frequency  u  form  an  area 

i 

on  the  log-magnitude/  phase  angle  chart,  the  Nichol's 
Chart(N.C.),  Figure  III. 2-4.  In  this  respect  it  matters  not 
whether  the  plant  has  one  state  or  a  thousand.  This  means 
that  for  the  frequency  in  question  all  possible  plant  phasors 
lie  within  this  area.  The  size  and  shape  of  the  area,  known 
as  the  plant  template, is  a  function  of  frequency.  It  is 
assumed  that  the  plant  does  not  have  any  RHP  zeros(non- 
minimum  phase  zeros).  RHP  zeros  limit  the  maximum  loop 
oandwidth  which  can  be  used  in  design,  and  the  bandwidth  is 
the  factor  which  QFT  manipulates  to  yield  a  design  that  meets 
system  performance  tolerances.  Normally,  the  plant  templates 
start  out  as  a  vertical  line  at  DC,  fatten  out  as  the 
frequency  increases,  then  converge  to  a  vertical  line  as  w—  oot 
The  height(in  decibels)  of  this  line  is  given  by  the  high 
frequency  gain  uncertainty. 

Figure  III. 2-5  shows  various  plant  templates  as  a 

function  of  frequency  for  some  simple  plant.  The  shape  is  not 

typical,  it  could  have  been  an  amorphous  blob,  or  even 

several  disjoint  areas.  However  for  design  purposes  any 

separate  areas  should  be  combined  into  a  single  area.  In 

order  to  draw  the  plant  template  at  any  (j  simply  connect  the 

i 

periphery  plant  points  to  create  an  area  on  the  Nichol's 


disturbance  rejection  and  tracking  responses.  From  these 
Dounds  a  nominal  loop  transfer  function  is  derived. 


having  magnitude  and  phase  markings  outside  the  chart.  Thus 
w'  from  the  intersection  of  L(j  u)  with  the  M  contours  on  the 

N.C.  yield  data  points  to  plot  the  closed  loop  response  as  a 
function  of  frequency. 

f .  Nominal  Plant 

In  order  to  arrive  at  compensation,  a  nominal  plant  P 

o 

is  chosen  from  all  the  possible  plants.  This  plant  is  used  to 

develop  the  loop  compensation  G  from  the  shaped  loop  transfer 

function  L  ,  where  G  =  L  /P  .  This  plant  should  be  stable  (if 
o  oo 

possible)  and  should  lie  on  the  lower  left  hand  side  of  the 

plant  template.  This  location  results  in  loop  bounds  closer 

to  the  center  of  the  Nichol's  chart[i6].  Other  nominal  plants 

result  in  loop  compensators  that  work  as  well,  but  the  actual 

compensator  design  process  is  more  difficult  because  the 

bounds  are  located  farther  up  on  the  Nichol’s  Chart.  The 

nominal  loop  transmission  function  l  is  determined  by 

o 

L  =  G  P  ,  or  the  compensator  using  this  P  is  G  =  L  /  P  . 
oo  ooo 

g .  Derivation  of  Open  Loop  Bounds 

For  both  the  disturbance  rejection  and  tracking  parts  of 
a  control  problem  the  loop  transmission  must  satisfy  certain 
bounds  such  that  the  specifications  are  met  over  the  range  of 
uncertainty.  QFT  represents  these  bounds  as  lines  traversing, 
or  closed  contours,  on  the  Nichol’s  Chart.  The  loop 
transmission  must  stay  above  the  lines,  or  stay  outside  the 
closed  contours,  at  each  frequency.  Two  bounds  exist  for  each 
frequency,  tracking  and  disturbance  rejection,  with  the  loop 


■  *-  •- 


•»" 
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To  find  the  bounds  it  is  necessary  to  ensure  that  the 

output  of  tne  system  y  satisfies  both  the  disturbance 

rejection  and  tracking  specifications  of  Figure  III. 2-3  . 

First  of  all  tne  maximum  allowable  value  of  the  closed  loop 

system  frequency  response  M  over  its  oandwidth  is 

m 

determined.  Tiiis  nu moer  corresponds  to  an  M  contour  on  tne 

N.C.  .  Tne  closea-loop  era n sm i s s i on  ,  M  =  L / ( 1  +  L)  <  M  , 

L  m 

required  to  stay  outside  the  particular  M  contour 

corresponding  to  tne  cnosen  M  .  Tne  area  enclosed  by  tne  A 

m 

contour  is  referred  to  as  the  "Forbidden  Region".  As  an 

exam ole,  in  Figure  III. 2- 5  the  I  M  I  is  5  dS,  and  the 

'  m ' 

"Forbidden  Region"  is  enclosed  by  the  I  M  I  =  5  dB  contour. 

‘  tn  ‘ 

The  uncertainty  in  the  closed  loop  response  is  entirely 
due  to  tne  plant  since  F  and  G  are  known.  This  implies 
that  tne  uncertainty  in  L  is  the  same  as  P  .  The  uncertainty 
in  tne  Magnitude  of  tne  output  is: 


l.:i  UAV  v  j  w  i  g  =  Lm  [AT  (  j  oj  )  j  (3.2-5) 


on  e  re 

Lm,AT^co;_  =  L-i l  T  \JtjoiJ  -  i-.H  T  (joojj 

U  L 

=  AL..U  L/il  *  L )  j  (3.2-6) 


L  must  oe  found  suen  tnat  T  ( j  co )  remains  .vitnin  the 
acceptable  range  of  Figure  II I. 2-3.  To  ensure  this  the  plant 
temolate  at  each  frequency  is  first  placed  on  the  Nichol's 
Cha rt ( remember  tnat  the  olant  temolate  represents  the  range 
of  uncertainty  of  the  olants  at  any  given  frequency).  Thus 


>  .»  ~  P  J"  »-  I  .  JT  .  fc.  •  ■  - 


since  Lm[  L  ]  3  Lm[  P  ]  +  Lm[  6  ],  Ang[  L  ]  =  Ang[  G  ]  + 

Ang[  P  ],  and  G  is  LTI,  then  the  template  can  be  translated 

vertically  and  horizontally  to  determine  the  bounds  that 

ensure  that  the  performance  specifications  are  met.  The 

template, cannot  be  rotated,  only  translated.  With  the 

template  at  an  angle  marking  near  tne  left  hand  side  of  the 

N.C.,  vertically  translate  tne  template  up  or  down  until  tne 

T  and  T  values,  tne  largest  and  smallest  magnitude  M 
>n  a  x  m  i  n 

contours  respectively  that  the  template  is  touchi ngj sati sfy 
the  requirement: 


max 


m  i  n 


U 


(3.2-7) 


When  this  requirement  is  satisfied  mark  the  nominal  plant 
point  on  the  N.C.,  then  move  over  ten  degrees  and  reoeat  the 
process.  If  in  this  process  the  template  penetrates  the  M 

m 

contour,  the  template  must  'oe  translated  on  tne  M.C. 
until  it  is  out  of  the  contour  even  if  the  tracking  oounos 


are  not  met.  Since  tne  gap  oetween  T  (j«)  and  T  (Jcj)  oecomes 

U  L 

wioer  as  c*j— *«e»  tnis  translation  is  vertical.  At  loner 


frequencies  cue  template  can  literally  "riae"  tne  upper 


portion  of  tne  .4  contour  wnen  tnis  nappe  ns.  At  nigner 

m 

frequencies  tne  difference  in  nagnituue  can  result  in  the 


oound  completely  enclosing  tne  .4  contour.  From  this  steo, 


the  lines  representing  tne  track  i  n  g  oounos  3  (jw)  on  L  are 

r 

derived. 


To  derive  the  disturoance  bounds  the  output  under 
disturbance  is  considered.  Assuming  a  disturbance  inDut  of  d, 
the  closed-loop  output  of  the  system  is: 


3 


y  =  d  P/(i  +  L) 
d 


(3.2-3) 


This  can  oe  rearranged  until  it  is  of  the  form 


1 1  +  L 1 

>  d  P/b 

d 

where  b  is  the  upper  bound  of  the  disturbance  response, 
d 

At  this  point  the  modified  Nichol's  Chart(MNC)  is  used 
to  finish  the  bound  determination.  Flip  the  normal  Nichol's 
Chart  upside  down,  change  the  magnitude  signs  of  the  open- 

loop  grids(  -24  dB  becomes  24  dB)  and  reverse  the  open-loop 

O  0  0  o 

angle  markings(-60  becomes  -300,  -170  becomes  -190,  etc.). 

Now  it  is  possible  to  design  L  directly  on  the  MNC  with  the 

superimposed  lines  now  representing  the  magnitude  and  phase 

of  l/[  1+  L  ].  For  disturbance  rejection  the  phase  of  the 

signal  does  not  matter,  just  the  magnitude.  This  assumes  that 

the  disturbance  has  the  maximum  effect  at  every  frequency, 

i.e.,  worst  case  scenario.  Take  the  plant  template  and  place 

it  on  the  MNC.  Now  translate  the  template  until  its  lowest 

point  rests  on  the  contour  of  constant  Lm[  1  +  L  J  equal  to 

-Lm[  D(jw)  ]  at  that  frequency.  Shift  the  nominal  point 

along  the  chart  as  before,  recording  the  bound  for  the  nomial 

plant.  Repeat  this  for  each  template  in  the  plant  set.  Invert 

the  chart  and  translate  the  B  (jw)  bounds  to  the  normal 

d 

Nichol ' s  Chart. 

g .  Derivation  of  Loop  Transmission 

With  the  bounds  in  place  on  the  Nichol's  chart, the 
design  of  L  itself  is  started.  The  initial  choice  of  where  to 
start  depends  on  a  number  of  factors.  If  the  open-loop  plant 
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contains  an  unstable  pole  this  must  be  included  in  L  .  For 
steady-state  error  reduction  a  Type  1,  2,  3,  or  more  system 
could  be  desired.  If  so,  the  proper  number  of  poles  at  the 
origin  must  be  included  in  L  at  the  onset  of  design.  The 
resulting  L  should  barely  meet  the  bound  at  each  frequency  to 
give  minimum  bandwidth  compensation.  If  such  an  L  is 
achieved,  then  the  system  is  optimally  compensated( wi th 
respect  to  the  performance  specifications)  A  typical  L 
satisfying  a  set  of  typical  bounds  is  shown  in  Figure  III. 2-7 
The  actual  "shaping"  of  L  is  a  cut-and-try  process  where 
poles  and  zeros  are  added  and  shifted  to  conform  to  the 

bounds  at  each  frequency  .  Sometimes  a  complicated  bound 

i 

structure  can  result  in  very  compl i ca ted ( h i gh  order) 
compensation.  In  this  case  it  may  be  better  to  use  a  smaller 
order  compensator,  trading  simplicity  for  a  wider  bandwidth. 

h .  Loop  Compensation 

The  compensator  G  for  the  nomi na 1  plant  is  determined 
as  follows: 

G  =  L  /  P  (3.2-10) 

o  o 


i .  Prefilter  Design 

Once  the  loop  compensation  has  been  found  the  tracking 


response  must  be  shifted  to  give  the  specified  output  to  a 
command  input.  The  prefilter,  F,  shifts  the  frequency 


response  of  the  closed-loop  system(  L / ( 1  +  L)  )  until  it  lies 
within  the  tracking  response  bounds. 


in«< 

InYll 


Fig. III. 2-7:  Typical  Loop  Transmission 


Place  the  plant  template  for  a  given  frequency  u  on  L 

i 

with  the  nominal  point  resting  directly  over  the  top  of  the 

corresponding  frequency  point  on  L  .  Record  the  maximum  and 

o 

minimum  closed-loop  magnitudes  that  the  template  touches, 

Lm  [  T  ( j  w  )  ]  and  Lm  [  T  (jw  )].  Next  obtain  the  upper  and 
max  min 

lower  tracking  bounds  for  the  same  frequency,  Lm [  T  ( j u>  ) J 

U 

and  Lm  [  T  (jw  )]  as  in  Figure  III.  2-3.  Now  defining  two 
L 

quantities: 


A  =  Lm  [  T  (jw)]  -  Lm  [  T  (jw  )] 


(3.2-11 ) 


(3.2-12) 


and 


B  =Lm[T(jw)J-Lm[T  (ju  )] 
U  min 


The  response  of  the  prefilter  must  be  such  that 


A(jw)  <  F(joi) 


<  B  ( ju>) 


(3.2-13) 


A  good  way  to  do  this  is  to  make  an  Lm ( •  )  verses  <j)  plot  for 
botn  A  and  B,  and  then  chose  an  Lm  £ F  J  bounded  by  A  and  B 
(  Lm  [f]  plot  lies  between  A  and  B).  Figure  III. 2-8  shows 
one  possible  design  of  F  given  the  indicated  bounds.  This  F 
is  not  unique,  since  any  function  that  does  not  violate  those 
bounds  satisfies  them. 


j .  Summary 

This  completes  the  description  of  the  Q F T  design  method 
for  compensating  SISO  plants  utilizing  output  feedback  with 
wide  plant  uncertainty.  Using  the  design  technique  is 
straightforward  without  any  guesswork.  If  the  design  is 
followed  to  the  letter,  a  design  guaranteed  to  meet  the 
specifications  over  the  range  of  uncertainty  results.  Also, 
QFT  is  transparent,  allowing  the  designer  to  see  the  trade¬ 
offs  required  for  meeting  the  specifications  as  the  design 
process  procedes,  which  can  eliminate  time  wasted  on  tries  at 
tne  compensation  of  "u ncompe nsa ta b 1 e "  plants. 
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Reduction  of  MIMO  systems  into  Equivalent  SISO 
Systems 


1 1  I  .  3 


In  order  to  use  the  QFT  design  technique  for  MIMO 
systems  they  must  first  be  decomposed  into  a  matrix  of  SISO 
systems.  Dr.  Horowitz  has  developed  a  method  to  do  this 
using  fixed  point  theory.  This  section  describes  the  steps 
taken,  but  not  the  functional  analysis  behind  it,  leaving 
that  to  the  reader  from  references[  17,  18,  19,  20  ] , 

As  a  further  simplification  only  a  2  X  2  system  is  examined^ 
since  that  is  the  type  of  system  used  in  this  thesis. 


a.  Description  of  M I M 0  Plant 


The  general  description  of  the  MIMO  input-output  relationship 
i  s : 

l  =  £  u  (3.3-1) 

where 

£  =  vector  of  plant  outputs 
£  =  matrix  of  plant  transfer  functions 
u_  =  vector  of  plant  inputs 

£  is  derived  from  eitner  the  state-space  matrices  or  the 
system  linear  differential  equations.  The  plant  has  m  inputs 
and  m  outputs,  i.e.,  ^  and  u.  are  m  X  1  vectors.  The  plant 
uncertainty  is  bounded  with  £  being  a  member  of  all  possible 
plantsj£}.  Notice  that  the  input  and  output  dimensions  are 
the  same.  With  m  inputs  at  most  m  outputs  are  independently 
control  1 abl e .  If  more  than  m  inputs  exist,  the  system  is 
modified  until  it  is  m  X  in. 

For  this  2X2  system  the  MIMO  compensation  structure  is 
shown  in  Figure  III. 3-1  .  It  consists  of  the  plant  matrix  £  , 
a  diagonal  loop  compensation  matrix  G  ,  and  a  prefilter 
matrix  F  .  The  matrices  are  of  the  form: 


11 

P 

12 

21 

P 

22 

1 

0 

(3.3-2) 


(3.3-3) 


Fig. III. 3-1:  MIMO  Compensation  Block  Structure 


f 

f 

11 

12 

f 

f 

21 

22 

(3.3-4) 


Figure  III. 3 - 2  is  a  more  detailed  look  at  the  2X2  MIMO 

compensated  system.  In  general  £  can  be  fully  populated; 

however,  in  this  thesis  the  cross  coupling  terms  f  and  f 

12  2 

are  set  equal  to  zero.  This  is  done'  since  in  this  design  the 
velocity  input  commands  only  the  velocity  output,  with  the 
same  thing  true  of  the  flignt  patn  angle.  No  cross  coupling 
is  desired. 

The  plant  matrix  must  satisfy  a  few  constraints: 
a.  >  must  be  invertible.  In  other  words  £  must  not  be 
singular. 


F  i .  1 1  I  .  3  - 2 :  MIPIU  Compensation  Signal  Flow  Grapn 

b.  The  main  diagonal  terms  must  dominate  as  s — ►  «»  for 
al!  P  €  {  P]  ,  i.e.. 


L  i  m  [ 

1  p  p 

P  P  | 

]  >  o 

$— »Of> 

1  11  22 

12  21 

This  is  referred  to  as  the  diagonal  domi nance 
c  o  n  o  i  t  i  o  n  . 


c  .  Gener a t i on  of  Ejui valent  SI  SO  too  os 

Jit.n  £  ue  r'iiieu  and  nee  ting  t;ie  conditions  a  cove,  tne 
inverse  of  P  ,  P  ,  exists,  w n e r e 


-  i 
P 


A  new  matrix,  Q 


P*  P* 

il  U 

P*  P* 

21  22 


is  now  introduced  with  the  individual 


inverse  of  t in e  individual  elements  of  ? 


elements  being  tine 


The  MIMO  problem  begins  loop  shaping  by  designing 


first  the  loop  that  has  the  strictest  tolerances. 

These  can  be  due  to  tracking,  disturbance  rejection,  or  plant 

uncertainty.  Doing  so  results  in  the  smallest  loop 

transmission  band widths.  Once  the  decision  is  made  as  to 

which  loop  to  shape  first,  each  output  of  each  transfer 

function  of  that  row  is  examined  to  find  out  which  one  places 

tne  highest  oouiids  on  L.  From  these  bounds  L  is  shaped.  Tnis 

L  is  then  usee  to  define  a  modified  q  for  the  remaining 

i  i 

luop.  fn  i  s  last  loop  design  is  an  exact  solution,  minimizing 
tne  o  a  n  a  w i a  t  n  . 


f .  Suinina  ry 

Designing  compensation  for  a  MIMO  sytem  using  QFT 
requires  that  the  MIMO  system  pe  broken  down  into  a  set  of 
equivalent  SISO  plants,  with  the  interaction  in  the  MIMO 
system  represented  as  disturbance  inputs  in  equivalent  SISO 
systems .  Once  tne  decomposition  is  performed,  QFT  is  used  as 
j  e  s  c  r  i e  a  e  a  r  1  i  e  r . 


III. 4 


Introduction  to  Reconfigurability  Theory 


a .  Introduction 

As  is  seen  in  Chapter  IV,  QFT  can  be  used  to  design 
fixed  compensators  to  meet  tolerances  in  the  face  of  wide 
plant  parameter  uncertainties.  In  this  thesis  the  plant 
uncertainties  are  due  to  the  range  of  flight  conditions  and 
of  control  surface  failures.  The  advantage  of  fixed 
compensators  over  ones  that  are  scheduled  with  the  failure  is 
that  they  don't  require  failure  identification  with  its 
associated  time  delay.  If  acceptable  performance  is  possible, 
then  there  is  time  for  reliable  identification  after  which  an 
"optimum"  compensator  can  be  switched  in.  This  chapter 
develops  the  QFT  method  for  dealing  with  the  recon fi gu rab 1 e 
aspects  of  the  control  problem.  The  available  free  n  terms 
of  Chapter  II  are  chosen  so  that  the  healthy  aircraft  is 
basically  non-interacting(BNIA)  ,  i.e.,  the  off  diagonal  terms 
are  much  less  in  magnitude  than  the  diagonal  terms.  From  the 
theory  an  equivalent  MIMO  plant  matrix  is  developed  that  is 
used  in  Chapter  IV. 


d .  Reconf i gurabl e  System  Model 

The  feedback  control  structure  is  shown  in  Figure 

III. 4-1  .  The  plantiT  are  the  constituent  surface 

i  j 


Fig. III. 4-1:  Reconf i gurabl e  Plant  Structure 

deflection  vs.  output  transfer  functions  from  Chapter  II.  The 

objective  is  to  design  the  system  such  that  there  is  as 

little  interaction  as  possible  between  the  loops  for  the 

healthy(no  failure)  aircraft,  but  also  allowing  cross  coupling 

during  failures  which  results  in  a  redistribution  of  the 

remaining  control  authority.  The  p  are, chosen  to  give  this 

i’j 

type  of  interaction.  From  Figure  III. 4-1,  the  equivalent 
plant  matrix  is: 


p'  = 

P  + 

M  P 

{ 3 . 4—2  b ) 

12 

12 

12  11 

P‘  = 

?  + 

M  P 

( 3 . 4-2c  ) 

21 

21 

21  22 

P 1  = 

P  + 

M  *P 

( 3 . 4- 2d ) 

22 

22 

12  21 

c.  Choice  of  Weighting 

Factors 

For  an  ideal  non-interacting  system,  it  is  desired  that 


requiring 


P'  =  p'  >  0 
21  12 


p.  f  p. 

12  11 


p.  /  p. 

21  22 


(3.4-3) 


( 3 . 4-4a  ) 
( 3 . 4- 4b  ) 


These  m's  in  general  are  functions  of  frequency  since  in  the 

majority  of  systems  each  P'  has  its  own  set  of  poles,  or  at 

i  j 

least,  zeros.  An  exact  ^  may  not  always  be  the  best 

i  j 

solution,  especially  when  the  denominator  terms  in  equations 
(3.4-4a,b)  contain  non-minimum  phase(right  half  plane)  zeros. 
In  this  case  the  appropriate  m  contains  an  unstable  pole, 
and  requires  an  exact  cancellation  to  keep  the  system 
response  stable.  Since  exact  cancellation  is  impossiole  in 
reality,  this  choice  of  m  can  cause  problems.  In  this  case 
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a  minimum  phase  n  that  will  diminish  the  cross  terms  as 
much  as  possible  is  required  -  not  perfect,  out  practical. 
Note  that  under  failure  conditions  equations  3.4-4a,b  are  no 


longer  valid,  i.e.,  since  the  n  terms  are  derived  for  a 
particular  no-failure  plant,  failures  will  change  those 
plants,  eliminating  the  cancellation.  This  causes  the  plant 
matrix  to  contain  non-zero  off  diagonal  terms.  These  off 
diagonal  terms  lead  to  loop  interaction  which  allows  control 
surfaces  in  the  "healthy"  loop  to  assist  in  control  of  the 
loop  containing  the  failed  surface(s). 


1 1 1 . 5  Summary 

Chapter  III  presents  an  overview  of  QFT  as  used  in  this 
thesis.  Rigorous  mathematics  are  not  used,  rather  a  more 
informal  approach,  explanation  with  words,  is  used  to  give 
the  reader  some  idea  of  the  power  and  the  procedures  of  QFT. 
In  Chapter  IV  more  of  the  exact  equations  are  introduced  when 
used.  If  the  reader  wishes  to  delve  deeper  into  the 
background  behind  QFT  he  is  encouraged  to  read  the  references 
listed  in  the  bibliography. 


I  V .  Compensator  Design 

I  V  .  1  I  ntroducti on 

This  cnanter  uses  QFT  to  design  tne  necessary 

compensation,  g  ,  a  ,  f  ,  and  f  ,  such  that  the  design 
1  2  11  22 

soec i f i cati ons  are  satisfied  at  tne  three  flight  conditions  for  a 
set  of  possible  failures.  Tne  equivalent  plant  matricies  are 
ueveiopec  first.  Froi:  these  matrices  tne  £  matrices  are 
f  o  u  1 1  c* .  nsxt  tne  mOm  control  cnannel  is  developed  using  SISC 
aesiud  techniques.  T.ien  using  information  found  during  tne  design 
of  tne  hjh  channel,  tne  second  cnannel  design  is  completed. 

To  -jive  trie  proper  tracking  response,  tne  prefilters  are  tnen 
developed  u  s  i  ri  a  frequency  response  data  of  each  loop 
transmission  function.  During  tne  design  process  it  is  noted 
that  proolems  occur  due  to  the  aircraft's  conf i gu ra t i on .  The 
problems  stem  from  tne  non-ninimun  phase  characteristics  of 
the  flignt  patn  angle  responses.  These  problems  lead  to  a  re- 
evaluation  of  tne  control  system,  and  a  suggestion  to  ease 
tne  application  of  QFT  to  MIMO  systems  with  non-minimum  phase 
terms  and/or  unstable  open-loop  poles.  From  the  steps  taken  a 
plant  emerges  that  is  amenable  to  QFT  design  techniques.  It 
is  snown  tnat  certain  failures  can  be  tolerated  while  others 
leua  to  plant  matricies  for  wnicn  the  QFT  technique  cannot  be 
applied.  For  such  plants  a  modification  of  the  QFT  technique 
has  ueen  developed  to  acnieve  an  acceptable  system  design  for 
such  plants,  but  doing  this  is  beyond  the  scope  of  tnis 
thes  i  sill]. 


IV.  2  Plant  Equations 


The  first  step  to  design  compensators  for  the  control 
problem  at  hand  is  to  construct  the  equivalent  plant 
matricies,  £'  ,  from  which  the  £  matricies  are  obtained. 

These  plant  matrices  are  developed  from  the  state  equation 
models  in  Appendix  B.  Initially  the  transfer  functions 
relating  surface  deflections  to  outputs  are  used  to  construct 
plant  matrices  as  is  done  in  Chapter  II.  From  these  transfer 
functions,  proper  /x  to  satisfy  3 I A  system  requirements  are 

i  j 

founj.  Using  tne  /a.  ana  tne  plant  matrices,  anotner  set  of 

i  j 

equivalent  matrices  with  reconf i gurabi 1 i ty  "built  in"  are 
designed,  tfext  various  surfaces  are  failed  and  new  plant 
matrices  are  developed  at  each  failure  condition.  Each  matrix 
requires  examination  to  determine  whether  it  can  be  used  in 
the  QFT  design  process,  or  whether  it  contains  non-minimum 
phase  terms  which  cause  that  particular  flight  condition/failure 
not  to  be  included  in  the  plant  set.  The  equivalent  plant 
matrices  are  inverted  and  the  equivalent  SISO  systems  are 
developed.  For  tnese  equivalent  SISO  systems,  frequency 
response  oounds  at  each  frequency  are  asea  to  descrioe  tne 
uncertainty  in  the  plant.  Tnis  uncertainty  results  in  areas 
on  tne  Mi  cools  cnart  within  wnich  tne  plant  lias  at  eacn 
frequency.  Using  tnese  areas  and  oounos  constructed  wi  :n 
eicner  tracxinq  or  aiscuroance  rejection  criteria  (.vine never 
piaces  stricter  requ i remen ts  on  tne  loop  transmission,  i.e., 

.  "iitar  ounawiatn  ,  which  is  a  "cost  of  feedback”.  A  loop 

a  j  c '  .  *  *  a  s  t  '  -  ‘  ;  .  i  ■  '■  ■  -■  '  '  "  ■  t 


minimizes  loop  bandwidth  is  then  designed.  From  this  loop 

transmission,  the  required  compensation  g  is  constructed. 

i 

Prefilters  f  are  built  such  that  the  tracking  response 
i  i 

meets  the  performance  specifications  in  the  frequency  domain, 
and  thus  the  time  domain  specifications  are  also  satisfied. 

Two  condidtions  must  be  satisfied  such  that  the  design 
specifications  are  guaranteed  to  be  met. 

a.  The  inverse  of  the  equivalent  plant  P*  must  exist. 
This  is  a  controllability  condition  and  must  be 
satisfied  for  each  flight  condition  and  failure 
case.  Thus  the  determinant  of  the  plant  cannot 

be  allowed  to  become  zero. 

b.  The  second  condition,  the  diagonal  dominance 
condition,  requires  that 


11m  Q  Q  /  Q  Q 

s— »oo  LI  11  22  12  21 


(4.2-1) 


This  must  be  satisfied  for  each  £  £{£_[  and  all 

frequencies  greater  than  some  frequency  u>  , 

1  I  I  b 

wherever  the  quantity  |  1  +  L  >  1  ,  where  L  is  the 

loop  transmission  function. 


From  the  state  space  models  in  Appendix  B  transfer 
functions  relating  each  control  surface  input  to  each  output 
are  generated  using  the  Computer  Aided  Design(CAD)  package 
T0TAL[7  ].  The  resulting  transfer  functions  are  in  Appendix  C 
The  relationship  between  the  flight  path  angle,  pitch  angle, 
and  angle  of  attack  to  the  jth  control  surface  input  is 


r( S )/  6  (S) 


e( s)/  6  (s)  -  a(s)/  6  (s)  (4.2-2) 


The  actual  flight  path  equations  are  developed  directly  from 
the  state  space  representation  by  modification  of  the  C 
matrix.  The  aircraft  has  three  stable  and  one  unstable  root 
for  the  three  flight  conditions.  A  listing  of  these  roots  is 
given  in  Table  IV. 2-1.  Note  that  the  airplane  becomes  more 
unstable  as  forward  velocity  increases,  with  one  short  period 
root  moving  out  along  the  real  axis  and  tne  pnugoid  roots 

/ 

uproacning  tne  imaginary  axis.  Open -loop  instability  is  not  a 
problem  for  qFT  since  a  staole  closed-loop  system  is 
o e v e i o p e a  i n n e r e n t 1 /  o y  z n e  j e s i g n  technique. 

almost  imme a i a te 1 j  a  problem  arose  in  tne  assign 
process.  Tne  surface  to  flignt  pat.i  angle  transfer  functions 
are  all  non-minimum  pnase,  having  at  least  two  zeros  in  the 
3HP.  These  transfer  functions  always  resulted  in  a  non- 
minimum  piiase(nrnp)  equivalent  o  1  a n t  matrix.  The  sensitivity 
reduction  properties  of  a  nno  plant  is  inherently  limited,  so 
it  may  be  imoossiole  to  achieve  the  desired  performance.  One 
could  orocaea  using  the  "Singular  G"  cec unique  for  nmp 
plants [11],  put  tne  result  may  be  that  a  fixed  L T I ( linear, 
time  invariant]  compensator  cannot  handle  cne  uncertainty 
range.  Another  way  of  loo  King  at  tnis  is  tnat  tne  gain  cannot 
oe  increased  too  muc.n  since  an  open- loop  pole  will  migrate  to 
tne  nmp  zero  causing  instability.  From  discussions  witn  Dr. 
riorowicz,  tne  decision  was  maue  to  continue  tne  analysis 
using  a  different  uutpot  variable  otner  than  flignt  path 
angle  to  control,  wnicn  woulo  nopefully  provide  a  n.np  plant 


determinant.  It  is  determined  that  the  angle  of  attacx  always 
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has  minimum  pnase  transfer  functions.  Thus  the  output 
variables  chosen  are  the  velocity  and  the  angle  of  attack 
rather  tnan  the  velocity  and  flight  path  angle.  This  is  not 
ideal  for  aircraft  landing  control;  however,  it  does  have  the 
advantage  of  guaranteeing  beforehand  that  a  design  does  exist 
which  satisfied  the  specifications  over  the  range  of 
uncertainty.  No  sucn  guarantees  can  be  made  apriori  for  nmp 
pi  ants  . 

TA3LE  IV.  2-1 


Eigenvalues  of  the  Open  Loop  System 
for  tne  Tnree  FI  1  ant  Conditions 


Flig.ot  Condition  1 :  t  Sea  Level  100  kts) 

s  =  0.2646 

S  =  -0.01697  +  jG.2500 
s  =  -1.545 

Pliant  C  o  n  d  i  t  i  o  n  2 : ( S  e  a  Level  120  kts) 

s  =  0.3351 

s  =  -0.01104  +  j 0.2031 
s  =  -1.979 

Flight  Condition  3: (Sea  Level  130  kts) 

c  =  0.1603 

s  =  - .  0032  Jo  +  .’0.1090 
3 


i  na  next  step  is  to  cnose  wnic.o  control  surraces  are 
to  execute  tne  primary  control  over  wnicn  variaole.  Tnis  is 
ao  tie  o  j  examining  tne  ua  _;ni  tjue  of  tne  frequency  response  for 
eac.i  transfer  function,  ana  cnoosiny  tne  dominant  control 


surface  of  3a  cn  output  as  t  n  a  t  control  .v  n  i  c  n  has  tne  n  o  s ' 


effect  upon  that  variable.  The  stabilators  and  canards 
influence  angle  of  attack  about  equaly,  much  greater  than 
any  of  the  other  surfaces.  Likewise,  the  top  and  bottom  vanes 
have  the  greatest  effect  on  velocity.  The  ailerons  influence 
velocity  slightly  more  than  they  influence  the  angle  of 
attack;  however,  they  are  added  to  the  angle  of  attack 
control  to  give  that  channel  greater  roDustness  under 
failure.  This  thesis  assigns  canards,  stabilators,  and 
ailerons  as  the  primary  angle  of  attack  control,  while  the 
reversing  vanes  provide  primary  velocity  control.  This  means 
that  the  "crosstalk  terms"  are  the  velocity  change  due  to 
canard,  stabilator,  or  aileron  deflection,  and  the  change  in 
angle  of  attack  due  to  reversing  vane  deflection. 

Once  the  primary  surfaces  for  each  variable  are 
defined  the  amount  of  control  authority  each  surface  has  in 
controlling  the  assigned  variables  must  be  calculated.  The 
method  used  in  this  thesis  first  assigns  a  percentage  of 
control  effoct  to  each  control  surface  involved  in  the 
primary  control  of  a  variable.  The  percentages  for  the  angle 
of  attack(AOA)  are: 


Canards 

40% 

S  t  a  d i 1 ators 

50% 

Ai 1 erons 

10% 

These  percentages  were  suggested  by  Capt.  Greg  Mandt[s]and 
assume  the  STOL  aircraft  has  the  ability  to  symmetrically 
deflect  the  ailerons.  As  mentioned  earlier,  the  stabilators 
and  canards  influence  angle  of  attack  about  equally,  with  the 
stabilator  being  slightly  stronger.  The  small  percentage 


chosen  for  the  aileron  reflects  the  fact  that  this  surface 
must  also  be  used  to  control  lateral  variables  as  well. 
Choosirtgla”  large  percentage  could  cause  saturation, 
especially  if  large  lateral  manuevers  are  called  for.  For 
velocity  control  the  percentages  picked  are: 


Top  vanes 

50% 

Bottom  vanes 

50% 

In  order  to  achive  these  percentages  the  amount  of  control 

exerted  by  a  particular  surface  over  a  variable  is  assumed  to 

be  represented  oy  the  DC  gain  of  its  transfer  function.  The 

weignting  terms  4  from  Chapter  II  are  chosen  such  that  the 

i 

new  DC  gains  reflect  the  percentage  values  picked.  These  DC 
gain  terms  are  in  Table  IV. 2-2  for  the  three  flight 
conditions.  The  signs  of  these  A  terms  are  chosen  such 

that  all  the  control  surfaces  work  together.  This  is  done  by 
referencing  the  Bode  plot  for  each  transfer  function  and 
choosing  the  sign  that  results  in  similar  frequency  response. 
It  turns  out  that  simple  sign  changes  on  the  A  terms  ensure 
that  this  occurs.  As  an  example,  for  flight  condition  one  the 
DC  gains  for  the  angle  of  attack  are: 

Canards  1.349 

Stabilators  2.835 

Ailerons  0.3961 

As  a  first  try  4  is  set  to  unity.  This  implies  that  40%  of 

1 

tne  control  authority  is  equal  to  1.349.  The  negative  sign  is 
bropped(but  is  dealt  with  later).  This  implies  that  the  total 
effort  has  a  magnitude  of  3.373.  In  order  to  achieve  the 
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the  terms  are: 

Top  vanes  1.000 

Bottom  vanes  1.662 

Now,  to  make  sure  all  the  surfaces  are  working  together, 
i.e.,  are  in  phase,  signs  on  4's  2,  3,  and  4  must  be 
changed.  For  the  first  flight  condidtion  the  A  vector  is: 

T 

A  =  [  1.000  -0.5948  -0.8514  -1.000  1.662] 

The  A  vectors  for  all  flight  conditions  are  in  Table  IV. 2-3 
At  this  point  the  reader  might  ask  if  choosing  a  A_  vector 
at  each  flight  condition  is  not  actually  scheduling  A  » 
implying  that  the  design  is  not  really  robust.  The  intent  is 
examine  the  plants  resulting  from  each  of  the  vectors,  and 
select  the  A  vector  which  results  in  the  greatest  number  of 
minimum  phase  P'  under  failures  for  all  three  flight 
conditions. 

TABLE  IV. 2-3 

Weighting  Factors  for  the  Three 
Flight  Conditions 


FI i ght 

condi ti on 

1 

2 

3 

Deltas 

A 

1 

A 

1.000 

1.000 

1.000 

-0.5948 

-0.6357 

-0.7191 

2 

A 

-0.8514 

-0.5681 

-0.2303 

Using  these  A  values,  equivalent  plant  matricies  £  are 
developed  for  each  flight  condition.  Here  a  problem  occurs 
because  of  the  structure  of  the  velocity  to  top  and  bottom 
vanes  transfer  functions.  For  instance,  at  flight  condition 
one  the  velocity  to  top  vane  transfer  function  is: 


v ( s  )  =  -9.496(s  +  0. 2 3 28  + j 0.4 762  )  (s  +  4.979)  (4.2-3  ) 

JJJJ  (  s +  d.  01 89 7  + j 0.2 508) (s  +  1. 5 46  ) (s+0. 2846) 


while  the  velocity  to  bottom  vane  transfer  function  is 


v  (  s )  =  -9.496(s  +  0. 5715  +  jO. 1485) (s-2.415 

atTT  (  s  +  0".  01897+ jO .  2  508 7  ( s  +  1  .  546  )  ( s - 0 . 2&46) 


(4.2-4) 


Using  the  proper  A  vector  from  Table  IV. 2-3,  the  two 

transfer  functions  are  added  together  to  form  P  ,  i.e., 

22 


22 


v ( s )  =  -3.143(s-12.29)(s+0.4202+j0.4092)  (4.2-5) 

TTsT  ( s  +  0. 0189  7  +  j 0.2 5  08 )( s  +  1 .546  )  (s-C.  2846) 


which  is  non-minimum  phase.  Using  this  P  transfer  function 

22 

results  in  //-  values  with  RHP  poles.  For  this  transfer 


function  to  be  minimum  phase  the  condition 


45<  *4 


(4.2-6) 


must  be  satisfied.  However  A  must  not  be  too  small.  This  „ 

5 

places  much  more  control  emphasis  on  the  top  vanes,  possibly 

leading  to  their  saturation  during  a  maneuver  or  after  a 

failure.  This,  in  turn,  can  lead  to  the  loss  of  aircraft 

control  from  effects  of  nonlinear  el ements ( such  as  limit 

cycles).  As  a  compromise,  A  is  chosen  as  0.926.  This 

5 

delta  value  results  in  a  minimum  phase  system  and  also  leads 


The  servo  transfer  function  for  the  top  and  bottom  vanes  is 
i  (s)  =  89 

TJsfT  r^T-59  (4.2-9) 

With  the  servo  transfer  functions,  the  *£  matrices  for  the 
three  flight  conditions  are  calculated  using  the  modified 
vectors  and  are  in  Tables  IV. 2-5  through  IV. 2-7.  Now  that  a 
satisfactory  plant  model  is  designed,  QFT  reconfigurability 
theory  can  be  applied  as  discussed  in  Chapter  III. 


I V . 3  Tracking  Specifications 

In  a  control  problem  one  is  seldom  confronted  with  the 
situation  where  a  specific  response  to  an  input  is  required 
Usually  a  range  of  acceptable  responses  is  specified. 

QFT  provides  the  minimum  bandwidth  compensators  which 
guarantee  that  the  closed-loop  response  remains  within  the 
assigned  tolerances  over  the  given  range  of  the  plant 
parameter  uncertainty.  The  acceptable  responses  are  often 
defined  in  either  the  time  domain  or  the  frequency  domain, 
they  are  in  the  former  they  can  be  converted  to  the  latter. 
This  step  is  essential  since  QFT  design  Is  performed  in  the 
frequency  domain. 

From  data  provided  by  McDonnell  Douglas[6  ],  ranges  of 
frequency  response  for  both  velocity  and  flight  path  angle 
are  defined.  For  the  velocity,  the  desired  response  is  a 
single  order  system 

T ( s )  =  5a  (4.3-1) 

s  +  a 

where  'a'  can  range  Detween  0.6  and  1.0  .  The  lower  bound 


94692 (s+34. 51 )  (  s  +  0 . 0 1 85Q±j 0 . 2563)/ A 

1 

• 

-6. 668(s  +  16.87)(s  +  0. 01834 ±jO. 2576)  /  * 

2 

-6438000 (s+26. 41) (s+0. 3369±j 0.4395)/  A 

1 

62.54(s+88. 7)  (s  +  0. 366 l±j 0.4445  )/  a 

2 

(s+30.62) ( s+133.6*j235. 1)( s -0.2846) (s+1 .546) 
(s+0. 01897±j0. 2508) 

(  s  + 1  .  546  )(s-0. 2846 )( s  +  0. 01897±j0. 2508 )(s+89) 


9 502  5  ( s  +  46 . 83 ) ( s  +  0 . 01 510± j 0 . 2 1 78 )  /  4 

3 


-5.269(s+19.39)(s+0.01480±j0.2182)/  4 


-7008144(s+26. 33) (s+0. 6516* j 0.2793)/  4 

3 


44 . 64( s  +  91 . 19 (s  +  0. 697  3. j0±2  348)/  4 


(  s  +  30 . 62  ) ( s+1 38±  j235.1)(s+1.979)(s-0.3851) 

( s+0. 01104±j0. 2081 ) 


(  s  +  1  .  9  79  )  (s-0. 3851  ) (s  +  0. 01 104±j0. 2081  ) ( s+89 ) 


1227  1  l(s  +  79. 60) { s  +  0 . 0l05  0±j0 . 1459)/  4 

5 

-3.505(s+27.24){s+0. 01 060±j 0 . 1 4  72 ) /  4 

6 

-6146000(s+32.98)(s+0.3911±j 1.084)/  4 

5 

20. 06(s  +  92.29)(s  +  0.4517±j 1.150)/  4 

6 

(s+30.62) (s+138±  j  235  .  1  ) (s -0.4608) (s  +  2.686  ) 
*(s+0.003206±j 0.1090) 

(s  +  89)  (s-0. 4608)  (s  +  2.686)  (  s  +  0 . 003206dkj  0 . 1090) 


transfer  function  is 


—  I  (s)  =  3.0  (4.3-2) 

u  s  +0.6 

while  the  upper  tracking  bound  is 

l  (s)  «  5.0  (4.3-3) 

U  s  +1.0 

In  the  time  domain  T  (s)  has  the  step  response 

L. 

specifications: 

Rise  time(t  )  ■  3.662  s 

r 

Settling  time(t  )  *  6.520  s 

s 

Peak  value(M  )  «  5.000  ft/sec 

P 

Final  value(F  )  *  5.000  ft/sec 

v 

Likewise,  the  upper  bound  specifications  are: 

Risetime  *  2.197s 

Settling  time  *  3.912  s 

Peak  value  »  5.000  ft/sec 

Final  value  *  5.000  ft/sec 

Figure  IV. 3-1  shows  the  time  responses  for  both  the  upper  and 
lower  velocity  response  bounds. 

It  is  essential  in  any  practical  design  that  at  large 
the  permissable  range  of  |  T(jw)|  exceed  the  range  of 
|  P  ( J  w )  |  .'This  is  always  guaranteed  if  an  infinite  difference 
between  the  upper  and  lower  bound  is  specified  asio— « 

Stated  more  succinctly: 

1  i  m  [  Lm  [t  (s)  -  T  ( s  )]  ]  =  00  (4.3-4) 

s-oo  L  U  L  J 

One  way  to  do  this  that  works  well  is  to  add  a  zero  to  T 


.  This  gives  an  order  of  two  difference 


and  a  pole  to  Tw 

between  the  bounds  and  ensures  an  infinite  spread  aso>-*ao  . 

The  required  poles  and  zeros  for  this  purpose  are  added  as 

close  as  possible,  subject  to  the  time  response  being 

negligibly  affected  -  this  is  easily  done  on  the  computer. 

The  first  choice  can  be  very  far  off.  If  no  great  deviation 

in  time  response  is  noticed,  the  pole  and/or  zero  can  be 

brought  in  toward  the  origin.  This  is  continued  until  an 

excessive  deviation  in  time  response  results,  after  which  it 

is  backed  out  a  bit.  This  thesis  uses  a  5  percent  deviation 

in  tne  settling  time  and/or  rise  time  as  the  limit  of 

pole/zero  movement.  For  the  velocity  reponse  upper  bound  a 

zero  cannot  be  added  realistically  since  it  is  a  first  order 

transfer  function,  instead  two  poles  are  added  to  T  to  give 

L 

the  same  effect.  These  modified  upper  and  lower  bound 

velocity  transfer  functions  are: 

T  (s)  =  5.0  (4.3-5) 

v  (s  +  1) 

U  _ 

and 

T  (s)  =  2100  (4.3-6) 

v  (s  +  0.6)  (s  +  20)  (.s  +  3-5) 

L 

Tne  resulting  frequency  response  bounds  for  velocity  are  in 

Figure  IV. 3-2  .  The  range  between  T  and  T  gives  the 

WU  L 

set  of  allowaDle  frequency  responses  that  result  in  the 
desired  time  responses. 


F(T)  RESPONSE 

,_p.00  1.25  2.50  3.75  5.00 

i.  i  i  i  i  i.  ■  i  »  I  ■  ■  i  i  ■  »  ■  ■  »  »  »  ■  ■  ■  i  i  *  .  .  I  ..lll.i  11  1  .  j _ l 


TRACKING  RESPONSE  FOR  VELOCITY  CHANNEL 


AOA  bounds  are,  for  the  lower  bound: 


Rise  time 

=  0.5859 

Settling  time 

=  1.062 

Peak  value 

=  3.500° 

Final  value 

=  3.500° 

and  for  the  upper  bound 

Rise  time 

=  0.1389 

Peak  time 

=  0.3354  s 

Settling  time 

*  1.098  s 

Peak  value 

0 

=  4.582 

Final  value 

=  3.500° 

The  plots  of  the  time  response  bounds  for  AOA  are  in 
Figure  IV . 3-3  . 

As  is  done  earlier  with  the  velocity  response  the 

transfer  functions  are  modified  to  provide  an  infinite 

separation  of  T  and  T  as  ui-xja.  A  zero  is  added  to  T  and 

°L  “U  U 

a  pole  to  T  .  The  resulting  bounds  are: 

<* , 


T  (s)  =  7 ( s+50 )  (4.3-10) 

U  s  2  +  7  s  +100 

T  (s)  =  3920  (4.3-11) 

a L  ( s  +  4) ( s  +  14) ( s  +  20 ) 


Figure  IV. 3-4  shows  the  resulting  frequency  bounds. Again  the 
bounded  area  represents  the  range  of  frequency  responses 
which  satisfy  the  given  specifications.  The  next  step  is  to 
determine  the  bounds  on  disturbance  rejection. 


TIME  RESPONSE  BOUNDS  FOR  flOfl  CHANNEL 


Since  the  type  of  aircraft  this  thesis  uses  as  a  basis 
for  reconf i gurabl e  controller  design  is  just  in  the 
preliminary  design  stages  itself,  the  basic  control  laws  are 
still  in  development.  The  desired  disturbance  rejection  has 
not  been  established  at  the  time  of  this  writing,  but 
from  discussions  with  FDL  personnel  (.22]  the  level  has  been 
set  at  t  =  It  <  - 1 0 d B  for  all  flight  conditions  for 

I  12  I  21 

the  healthy  aircraft.  Under  failure  the  ideal  situation  is  to 
meet  the  same  bounds;  however,  the  remaining  surface 
deflections,  and  the  required  compensator  bandwidth,  may 
become  too  large  for  a  realistic  design.  In  this  case  the 
disturbance  rejection  bounds  are  sacrificed  for  stability  in 


unr 


Def  i 


.  b  Aj 
W 

1  C  J  1 

e  Q  i 
SO  s 


nd i  t 


Configuring  the  plant  in  this  manner  leads  to  smaller  loop 
transmission  bandwidths.  A  method  using  ideal  crossfeed 


terms  n  and  //  is  developed  in  Chapter  III.  These  ideal 
21  12 

H  terms  are  defined  in  equations  3.4-1  through  3.4-7  . 
i  j 

Using  the  equivalent  plant  matrices  in  Tables  IV. 2-3  through 

IV. 2-7  the  individual  //  terms  are  calculated.  As  a  sample 

i  j 

calculation,  at  flight  condition  one  the  n  term  is 

12 

-5 

U  =  7.042(10  )(s+16.87)(s+30.62)(s+138.6+j235.1)  (4.5-1) 

12  ( s  +  34 . 5 1 ) ( s  +  89  ) 


Likewise  the  fi  term  is: 

21 

U  =  102942 (s+26.41)  (4.5-2) 

21  (  s  +  3  CL  6  2  ")  (  s  + 1 3  8 . 5  +  j  2  3  5  . 1) 

The  rest  of  the  fi  ' s  are  developed  in  the  same  manner  and 


are  in  Tabl  e 

IV. 5-1 

.  With  the  ii 

values 

the  plant 

matrix 

becomes : 

* 

P  + 

/i  P 

P  + 

(i  P 

11 

21  12 

12 

12  11 

P  '  = 

**S|» 

(4.5 

P  + 

//  p 

P  + 

1/  p 

21 

21  22 

22 

12  21 

Using  Eq(4.5-3)  the  equivalent  r ec onf i gu ra b 1 e  plant  matrices 

composing  _P'  for  the  first  flight  condiditon  are: 

P1  =  1.639(s-0.4618)(s+5.772)(s+134.2+j221.3) 

11  (s+30.62)(s+138.6+j235.1)(s+89) 

( s  +  23.45)  (s  +  7.257)  (4.5-4) 

''  (s  +  1".  5'4"6)  (s-0.284"6)  (s  +  0. 01897  + jO.  2508) 


Three  Flight  Conditions 


405.9(s+0.2524+j0.4044)(s+4.912 


(s  +  140.3+j240.2H s+1.642 


Note  that  these  terms,  even  though  they  are  minimum  phase, 
have  resulted  in  a  non-minimum  phase  plant  transfer 
functions.  It  turns  out  that  the  plant  determinant  is  also 
non-mi  nimuin  phase,  but  this  stems  from  the  choice  of  A 
rather  than  n  terms(since  these  can  be  factored  out  of  the 
determinant  expression).  This  also  occurs  using  the  other 
vectors  as  shown  in  Table  IV. 5-2.  Non-minimum  phase 
determi nants  .  1 imi t  the  loop  transmission  bandwidth  and  can 
cause  instability  using  the  MIMO  QFT  design  process.  For  either 
problem  re-examination  is  required,  or  an  extension  of  QFT 
for  non-minimum  phase  plants,  the  "Singular  G"  method  should 
be  used[23]«  Singular  G  is  beyond  the  scope  of  this  thesis, 
thus  the  A  vector  is  examined  to  see  if  it  can  be  used  to 
drive  the  system  minimum  phase. 

Upon  a  close  examination  of  the  aircraft  and  the  related 
equations  of  motion  and  resulting  matrices,  the  main  cause  of 
problems  seem  to  be  the  close  coupling  of  v  and  a  within, 
the  plant  and  their  respective  signs.  With  the  current  scheme 
of  choosing  the  A  vector,  as  explained  above,  the 
equivalent  plant  before  the  reconf i gurabl e  y  terms  are 


added  is: 


Flight  Condition  One: 


P'  =  1.64(s-0.462)(s+5.77)(s+134±j221)(s+23.5)(s+7.26)/  A 
11  1 

p '  =  P '  =0 
12  21 

P'  =  405(s  +  0.252tj0.404)(s  +  4.91)( s  + 140± j 240 ) ( s  +  1 . 64 ) /  4 
22  1 


Flight  Condition  Two: 

P*  =  2.42(s  +  0.892)(s  +  0.049±j0.407)(s  +  32.5)(s  +  1331,j247  )/  a 
11  2 

p  •  =  p  •  =0 
12  21 

P'  =  128(s-0.484)(s+2.33)(s+0.385)(  s  +  1 352“j  245  )  (  s  +  35 . 4 )  /  a 
22  2 


Flight  Condition  Three: 

P'  =  3.15(s+0.0052±j0.183)(s+42.9)(s-0.688)( s  +  1 34±j  222  ) /  a 

11  3 

P  1  =  P '  =0 

12  21 

P'  =  164(s-0.476)(s+0.89)(s+2.53)(s+42.5)(s+134±j237)/  4 
22  3 


A  =  ( s+30 . 6 ) ( s+138±j  235 ) (s+89)(s  +  1.55) 

1  *  (s-0.285)(s+0. 0189 ±j 0.251) 

A  =  (s  +  30.6) (s  +  138  ± j  235) (s  +  89) ( si .98) 

2  *  (s-0.385)(s  +  0. 01 l±j 0.208) 

A  =  (s  +  30.6)(s  +  13  8±  j  235) ( s  +  89)  (  s  +  2 .69) 

3  *  (s-0.461)(s+0. 00321 ±j 0.109) 


where  the  IP  (operation  causes  the  leading  gain  value  in 
•  ijl 

the  transfer  function  to  be  greater  than  zero(an  'absolute 


value'  type 

function).  As 

stated 

earlier 

,  the  ii 

values 

are : 

II 

C\i 

tH 

^4. 

'  -\~PJ 

lpi.l  = 

1-  J 1 

P 

11 

, 

( 4 . 5-8a  ) 

^21 

■  •  \\x 

P  1  = 
221 

“Of 
ro 
►— » 

|P 

1  22 

( 4 . 5-8b  ) 

Using  these,  the  reconf i gurabl e  equivalent  plant  is: 


Note  that  to  zero  out  the  off  diagonal  terms  of  £'  a 

subtraction  -ef  positive,  minimum  phase  transfer  functions, 

occurs  to  develop  each  P*  .  The  different  dynamics  of  the 

i  i 

velocity  control  channel  compared  with  tjie  AOA  channel 
seem  to  result  in  the  non-minimum  phase  plants  when  the 
addition  using  /j  '  s  is  accomplished.  In  order  to  circumvent 
this  problem  various  combinations  of  different  4  vectors 
that,  see  from  the  surface,  might  give  minimum  phase  plants 
are  chosen.  The  det[  £'  ]  is  calculated  for  various  failure 
conditions  and  4  vectors,  again  without  finding  a  minimum 


phase  £' .  In  each  case  at  least  one  zero  is  in  the  RHP. 
At  this  point  it  becomes  evident  that  some 


problem modification  is  required  to  provide  a  solution  to 
this  problem  using  QFT  without  using  the  "Singular  G" 
method.  As  a  simplification  measure,  the  second  order  servo 
poles  are  deleted.  The  assumption  is  that  the  final  loop 
transmissions  will  have  a  small  enough  bandwidth  not  to 
excite  the  neglected  servo  poles,  or  have  the  additional 
phase  lag  drive  the  system  unstable.  To  guarantee  that 
combinations  of  control  surfaces  will  result  in  minimum  phase 
plants,  the  plant  determinant  must  not  have  RHP  zeros.  The 
determinant  of  £'  is  given  by 

P'  P'  -  P'  P*  (4.5-10) 

11  22  12  21 


Factoring  out  the  reconf i gurabl e  terms,  the  servo  transfer 
functions,  and  the  square  of  plant  poles,  results  in  the 
determinant  being 

det  [  P  *  ]  =  K(  s  +  30 . 62  )  (  s  +  89  )  (1-  //  /i  ) 

(plant  poles)2  12  21 

(4.5-11) 

[  N ( P  )N(P  )  -  N ( P  ) N ( P  )  ] 

11  22  12  21 


Where  N(.)  stands  for  the  numerator  of  the  transfer  function 

in  question.  Only  the  term  [  N(P  )N(P  )  -  N(P  )N(P  )]  can 

1  1  22  12  21 

contain  non-minimum  phase  terms(assuming  1-  fi  fj  is  mi  nimum 

12  21 

pnase).  This  quantity  expanded  out  with  the  A_  vector  terms 
arid  separate  surface  transfer  functions  (Appendix  I)  is: 


(  P  P  -  P  P  )  +  tff(  P  P  -  P  P  ) 


11 

24 

14 

21 

11 

25 

15 

21 

+  A  / 

4  [ 

P  P 

- 

P  P  +  \p[ 

P  F 

> 

-  P  P 

2 

1 

12 

24 

14  22 

12 

25 

15  22 

+  A  / 

4  [ 

P  P 

- 

P  P  +  lp{ 

P  F 

> 

■  P  P 

3 

1 

13 

24 

14  23 

13 

25 

15  23 
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(4.5-12) 


where  if/  =  4  /  4  • 

5  4 


The  separate  difference  terms  can  be  combined  to  simplify  the 
equation: 

A  +  if/  B  +4/4  (A  +^B) 

1  12  12  2 


+  4/4  (A  +  ^B)  (4.5-13) 

3  13  3 

There  are  three  unknowns  in  this  expression:  if/ ,  A  /A  »  and 

2  1 

4/4  .  The  nominal  plant  used  to  develop  4  is 

3  1 

FC2:No  Failures.  The  /  value  is  chosen  to  make  as  many  of 

the  individual  A  ,  B  pairs  minimum  phase  as  possible.  Then 

k  k 

the  remaining  A  /  A  anc^  4/4  terms  are  picked  to  de- 

2  1  3  1 

emphasise  any  remaining  non-minimum  phase  tendencies.  The 
vector  chosen  is: 


4  =  [  -1  1  1  1  -0.6  ] 

This  vector  also  results  in  minimum  phase  plants  at  the  other 

two  flight  cond i t i ons ( no  failures).  Variations  of  individual 

4  terms  ±20%  still  resulted  in  minimum  phase  plants.  In 
i 

addition  the  plants  remained  minimum  phase  for  the  following 
failures  at  each  flight  condition: 

a.  Canards 


b .  Ailerons 


c .  Stabi 1 ators 


d.  Stabilators  and  Ailerons 


e.  Stabilators  and  Canards 


Also  at  FC2  failure  of  the  bottom  vanes  results  in  minimum 


phase  P'.  Top  reversing  vane  failure  never  results  in  a 
minimum  phase  plant.  Other  4  vectors  have  been  tried. 


h ever _th.e  above  4  vector  gave  the  most  minimum  phase 

plants  with  failures  out  of  all  A  tried.  In  other  word, 

that  particular  4  results  in  the  'robustest'  behavior.  Even 

with  the  plant  being  minimum  phase  the  diagonal  elements  were 

not  dominant,  signified  by  the  determinant  being  negative. 

For  diagonal  dominance  the  columns  are  swithched,  and  the 

signs  of  the  individual  4  terms  changed  to  flip  the  sign  of 

i 

the  determinant.  Physically  this  means  that  the  role  of 
primary  control  of  the  variables  has  changed.  The  aerodynamic 
surfaces  are  now  controlling  velocity  while  the  reversing 
vanes  are  controlling  AOA.  This  result  is  interesting  since 
it  is  opposite  of  what  one  would  think  is  logical.  The  switch 
seems  to  stem  from  forcing  P'  into  minimum  phase  behavior  to 
use  QFT(but  not  caused  by  QFT  directly).  Why  this  happens  is 
not  clear;  however,  for  small  signals  the  approximation 
should  still  be  valid. 

The  A  vector  is  now: 


A  =  [  1.0  -1.0  -1.0  -1.0  0.6] 

The  P'  matrices  used  in  loop  design  are  in  Appendix  D  along 
with  the  equivalent  Q  matrices  derived  from  them. 

A  moral  to  this  story  is  that  when  dealing  with  a  MIM0 
QFT  design  problem  the  very  first  step  taken  should  be  to 
find  the  range  of  weighting  functions(  in  this  thesis) 
that  results  in  minimum  phase  plant  matrices  before  even 
thinking  about  division  of  control  authority(or  anything 
else).  Without  this  characteristic  the  problem  becomes  much 
more  difficult.  If  this  had  been  done  in  this  chapter  several 


*  *  a  -  v  ■  Oj.  \ 


weeks  of  work,  as  well  as  ten  pages  of  writing,  could  have 
been  spared.  Had  this  fact  been  brought  out  at  the  beginning 
of  this  effort,  possibly  a  A  vector  could  have  been  found 
to  enable  control  of  the  flight  path  angle. 


IV. 6  Loop  Bound  Design  L 


Since  it  is  desired  that  only  AOA  commands  should 

command  AOA,  and  likewise  with  velocity,  the  equivalent 

closed-loop  SISO  structure  is  in  Figure  IV. 6-1  .  Notice  that 

f  and  f  have  been  set  to  zero  to  reflect  this.  To  begin 
12  21 

the  design  of  the  loop  transmissions,  the  designer  must  first 

decide  which  loop  to  design  first.  Loop  1,  the  AOA  channel, 

has  less  uncertainty  and  smaller  associated  disturbance 

rejection  bounds(as  determined  by  looking  at  both  the  size  of 

the  plant  templates  and  the  area  between  tracking  response 

frequency  bounds),  thus  it  is  designed  first.  This  is  an 

empirical  decision  since  the  same  steps  are  used  to  design 

both  loops;  however,  beginning  with  the  loop  having  the 

stricter  tolerances  leads  to  a  design  with  smaller  loop 

bandwidth.  Plant  templates  for  both  Q  and  Q  are 

11  22 

drawn,  and  from  this  it  is  seen  that  Q  has  less  area 

11 

enclosed  by  the  contour  for  every  frequency,  thus  this  loop 
is  designed  first. 

The  loop  transmission  L  must  satisfy  all  its 

1 

frequency  bounds(deri ved  from  the  plant  uncertainty  and 
performance  tolerances).  Examining  each  FC:Failure  case  it  is 


79 


found  that  the  transfer  functions  have  the  form: 


K  (LHP  Zeros) 

!HP  Pole)(LHP  Poles ) 


(4.6-1) 


where  K  is  a  positive  constant.  The  plant  template  for 

Eq(4.6-1)  is  a  line  at  low  frequencies,  becoming  wider  at 

middle  frequencies,  and  narrows  to  a  line  as  cj  tends  toward 
infinity. 


F i g  .  I V . 6- 1 :  Equivalent  SISO  Plants 


In  order  to  derive  the  frequency  bounds  for  the 

disturbance  rejection,  B  (jcu),  it  is  noted  that  the  output 

d 

y  is  due  to  disturbance  only  since  f  is  set  to  zero(a 
21  12 
change  in  AOA  for  a  velocity  input  is  not  desired).  The 


output  in  this  case  is: 

y  =  d  Q  /(l  +  L  ) 

12  11  1 

=  t  Q  /[Q  (1  +  L  )] 


(4.6-za) 


( 4 . 6-2b  ) 


Eli] 


Assuming  the  worst  case,  i.e.,set  t  =  b  .and  also  assume 

22  22 

that  the  response  is  less  than  or  equal  to  b  ,  the 

12 

disturbance  rejection  bound.  Then 


y  =  I  b  Q  /[Q  (1  +  L  )]  <  b  (4.6- 3a) 

12  I  22  11  12  1  12 

rearrangi ng 


1/(1  +  L  )  <  b  Q  /[b  Q  ]  ( 4 . 6  -  3b) 

1  12  12  22  11 


Tnis  must  be  investigated  for  each  flight  condition  and 

failure  case.  Q  and  Q  at  each  FC:Failure  case  are  used  to 
11  12 

find  the  most  demanding  bounds.  The  FC2 : Stabi 1 ator  Failed 

case  is  chosen  as  the  nominal  plant.  At  low  frequencies  it  is 

normally  located  on  the  left  side  of  the  template,  close  to, 

if  not  at,  the  bottom  point  of  the  template.  Such  a  choice 

makes  it  easier,  from  the  author's  experience,  for  design  of 

L  by  moving  the  loop  bounds  lower  on  the  Nichol's  chart. 

1 

Using  this  equation  and  theinverted  Nichol's  chart,  bounds 

for  l  can  be  found.  However,  if  the  frequency  response  data 
1 

shows  that  |  Ll  »  1,  then  the  above-  equation  is  simplified 

i 

to : 

L  >  bQ/ibQj  (4 .6-4) 

1  22  11  12  12 


or 

l  >  o  Q  /Lb  Q  j  (4.6-5) 

■  1  22  llo  12  IZo 

wnere  l)  and  l)  correspond  to  the  FC:Failure  that  causes 
llo  12o 

tne  worst  bound.  As  long  as  I  l  |»  1  ,  the  bounds  for  L 


are 


plotted  directly  on  the  Nichol's  chart.  When  |  L  j=  1  the 

1 

inverted  Nichol's  chart  is  required  .  With  the  log-magnitude 

measure(  Lm(x)  =  20*1og(x)  ),  L  must  satisfy 

1 

Lm(L  )  =  Lm ( b  )  +  Lm(Q  )  -  Lm(b  )  -  Lm(Q  )  (4.6-6) 

1  22  11  12  12 

From  the  speci f ications ,  b  =  -10  dB,  while  b  is  found  in 

12  22 

Figure  IV. 3-2  .  The  calculations  show  that  ll  +  L  |  is  very 

1 

small  for  even  low  frequencies.  This  implies  that  L  is  very 

1 

close  to  the  -1  +  jO  point  at  tnose  frequencies,  well  within 

the  +4  dB  contour  on  the  Nichols  chart  chosen  as  the 

forD idden  region.  The  reason  that  I  1  +  L  Jis  so  small  is 

1 

that  I  Q  |»|q  I,  which  makes  sense  sincelP'  }>>  IP'  J.  Thus  the 

1  12  1  if  22  12 

bounds  on  the  first  loop  are  determined  by  the  output  y 

11 

This  output  is  composed  of  two  parts:  tracking,  y  ,  and 

lit 

disturbance,  y  where 
lid 


y 

f  Q 

g  /d  + 

L  ) 

=  f  L  /(I  +  L  ) 

(4.6-7) 

lit 

11  11 

l 

l 

11  1  1 

y 

d  "  q 

/(I  +  L 

)=  b 

Q  /[Q  (1  +  L 

)J  ( 4 . 6-8a  ) 

lid 

li  li 

1 

21  11  12  1 

y 

y  + 

y 

( 4 . 6-8b ) 

li 

lit 

lid 

' 

Tne 

m  a  x  i  in  u  m 

value 

for  the  d  i  s 

turpance  is 

1 

v  M 

O  Q 

/  [  Q  L 

)|-|» 

/(Q  y  )  1 

(4.6-9) 

1 

llo1  1 

21  11 

12  1 

21  12  1  1 

for 

cases  where! 

L,  l»l* 

With 

1  Q  |  »  1  Q  1 , 

b  >  1 ,  and 

1  1 2 •  1  11* 

21 

L  | 

larger 

than  ) 

q  |at 

every 

frequency ( wh i ch 

should 

l1 

. 

1 

ll' 

be 

so  sinde 

l  = 

g  Q 

,  and 

realistic  loops 

that  provide 

1 

l  li 

sensitivity  reduction  tend  to  be  larger  than  the  plants  they 


contain),  the  disturbance  output  is  very  small  when 

compared  to  the  command  output  within  the  passband  of  f  and 

11 

can  be  ignored.  At  higher  frequencies  the  disturbance  may  be 

larger  than  the  tracking  since  y  may  be  smaller  than  y 

lit  lid 

To  meet  specifications  in  this  case  L  must  stay  out  of 

1 

the  forbidden  region.  This  is  done  by  translating  the 

templates  arounG  trie  forbidaen  region  and  noting  the  trace  of 

toe  nominal  plant.  This  defines  a  zone  tnat  L  cannot 

1 

be  witnin  at  that  frequency.  In  tne  tracking  case  with  the 


disturbance  small,  the  bound  for 
equation: 


L 

1 


satisfies  the  following 


T 

r 

i  C  1  b  -  a 

(4.6-10) 

max 

m  i  n 

1  “  «  11  11 

where  T  and 

r 

are  the  maximum 

and  minimum  closed  loop 

max 

magnitude  line 

m  i  n 

touched 

by  the  template 

respecti vel y.  Table 

IV. 6-1  contains  the  tracking  bounds  for  loop  one.  The  lower 


bound  is  dropped  in  magnitude  to  ease  the  bounds  on  the  first 

loop;  however,  any  simplifications  here  are  made  up  by 

wider  Dandwidths  of  f  and  the  second  loop's  compensation. 

11 

Tne  plant  templates  for  Q  are  in  Figure  IV. 6-2  .  The 

11 

nominal  plant  is  FC2:Cunard  Failed.  From  the  Nicnol's  cnart. 
tney  are  transferred  to  a  clear  sneet  of  plastic  and  using 

tne  previously  defined  tracking  oounds,  tne  bounds  on  L  are 

1 

drawn  and  lie  in  Figure  IV. 6-3  .  At  low  frequencies  the  bound 
for  co  =  i ,  dominates  since  it  is  well  above  the  bound  for 
w  =  2  and  cj=  4.  For  high  frequencies  the  bounds  wrap  around 


the  forbidden  region,  satisfying  both  the  tracking  and 
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Fia.IV.6-2:  Q  Plant  Templates 

.  ll 

disturbance  rejection  cases,  finally  resulting  in  an  UHFB 

24  d  3  hign.  In  designing  the  actual  L  a  choice  has  to  be 

made  as  to  the  snaps  and  magnitude  of  L  .  Since  the  open  ltjop 

1 

transfer  function  contains  an  unstaole  pole  the  Nyquist 
5 1  a  o  i  1  i  t y  Criterion  s  n  o  a  1  a  a  a  used  to  e  s  t  a  n  1  i  s  n  ,i n  a  t 

frequency  cna  rac  to  r  i  s  t  i  c  s  tne  loop  transmission  snould  have 

for  stability  unuer  feeubaCK.  A  polar  plot  of  an  open-  loop 
transmission  containing  an  unstaole  pole  wiiicn  is  stable  for 
unity  feeuoacK  is  snown  in  Figure  IV.  0-4  . 


Fig.  I  V . 6 - 3  :  Loop  One  Bounds 

For  a  stable  system  Nyquist  states  that  the  number  of  RHP 
zeros  must  be  equal  to  zero[l5  144]  .  The  equation  which 
relates  the  number  of  these  zeros  to  the  number  of  ooen  loop 
RHP  poles  and  the  number  of  encirclements  of  the  -1  point  is: 

l  =  N  +  P  (4.6-11) 

Z  =  Numaer  of  close-loop  RHP  zeros 


d  5 


wnere 


Fig. IV. 6-4;  Polar  Plot  of  Loop  Transmission  with 
an  Unstable  Pole  that  is  Stable  for  Unity  Feedback 


N  *  Numoer  of  clockwise  encirclements  of  the 
-1  +  jO  point. 

P  *  Numoer  of  RHP  poles. 


Since  P  *  1  the  -1  point  must  be  circled  -1  times 

(counterclockwise)  in  order  for  Z  =  0.  Since  g  should 

1 

contain  minimum  phase  terms  it  is  of  the  form: 

g  =  K  (  LHP  zeros)' 

I  {  LHP  poles)  (4.6-12) 


If  K  >0  then  only  the  minimum  phase  transfer  functions  can 

be  used  for  a  plant  with  one  unstable  pole  and  no  RHP  zeros. 

If  K<  0  then  the  plants  containing  one  RHP  zero  may  be 

staole,  'tuit  this  implies  that  Q  has  a  RHP  zero  which  cannot 

11 

be  used  with  tne  form  of  QFT  exploited  in  this  thesis.  This 


can  be  seen  if  a  pole-zero  plot  of  tne  system  is  made(Figure 


■  %  *■  •-  *  i  ■ 


IV. 6-5).  Using  unity  feedback  and  minimum  pnase  compensation 

in  the  forward  loop  does  not  change  the  fact  that  a  branch  of 

the  root  locus  starts  at  the  open-loop  RHP  pole  and  migrates 

to  the  RHP  zero  as  the  gain  is  increased.  The  system  is 

unstable  for  all  positive{  and  most  negative  )  gain.  One  of 

the  possible  stable  loop  transmi ssons (-Fi gure  IV. 6-4)  has 

a  starting  point  of  -130  degrees  at  cj  =0  .  Since  one  pole  is 

in  the  RHP,  tnis  means  tnat  this  pole  is  included  in  L  .thus 

1 

jiving  trie  -ldO  degrees  without  any  free  integrations.  Tnis  implies 

tnat  L  is  Type  U,  ana  tnat  some  steady-state  error  to 
1 

command  inputs  nas  to  be  tolerated.  Tne  amount  of  steady 

state  error  is  dependent  on  the  DC  gain,  the  higher  tne 

oetter.  But  a  hign  DC  gain  also  means  tnat  several  poles 

might  have  to  De  added  to  drop  l  fast  enough  to  just  meet  the 

1 

bound  at  u  *  1  to  minimize  the  bandwidth.  Tnis  results  in 
tne  forced  addition  of  several  zeros  just  to  bring  L 

1 

to  the  right  of  the  forbidden  region.  An  L  with  this  shaoe 

1 

is  of  fairly  high  order  in  both  the  numerator  and 

•denominator,  leading  to  an  even  more  complex  g  .  Also  the 

1 

system  is  conditionally  staole  at  oest,  ano  a  large  gain 
Change  .nay  drive  the  system  unstaole.  On  the  other  hand,an 
./ic.i  a  1  j.i3f  DC  gain  can  be  used  to  reduce  the  complexity  of 
u  silica  no  great  drop  nas  to  occur  to  meet  tne  oouna  at 

i 

cm  *  i,  out  tne  sceaay  state  error  increases.  Tne  smaller  JC 
•  a  in  is  on o sen  si  nee  tne  amount  of  steamy 


state  error  winc.i 


L 


l  1  C 


(4.6-14) 


wnere 


C(S)  =  (s  +  U  .  3d5 1  )  / ( s  -  U.3851)  14.6-15) 

C(s)  contains  tne  unstaole  pole.  The  zero  could  be  placed 

anywnere,  Dut  placing  it  at  the  'mirror'  image  of  the 

unstable  pole  is  analytically  convenient.  This  can  be  seen 

from  a  frequency  response  plot  of  the  above  transfer  function 

C(jw)  shows  a  straight  line  undergoing  a  180  degree  phase 

shift  as  w — ►  «  from  -180°  to  0°.  This  shifts  L  '  from  0° 

Q  1 
to  -180  at  U)  -  0  to  give  the  shape  of  L  required  by  the 

1  o 

Nyquist  stability  criterion.  L  '  is  designed  to  start  at  0 

1 

and  satisfy  the  requirement  that  when  it  is  multiplied  by 
C  ( s )  an  L  is  formed  that  satisfies  all  bounds.  For  the 


oounas  in  Figure  I V  -  6  -  3  ,  an  L  tnat  satisfies  tne  oounjs  ;5 

la 

6 

L  *  _ 1.313(10  )  (  s  +  -0 . 3o5 1  ) _  (4.5-16) 

lo  is  -  U.3o31)(s  ♦  0.2)(s  ♦  200)(s  +  600) 


One  can  see  ino w  well  L  .neets  its  pounds  o y  examining  Figure 

lo 

IV. 6-6  .  L  is  close  to  optimum  at  1  rad/s  and  again  at 
lo 

20  rad/s,  but  far  from  the  bounds  at  2,  4,  and  3  ra d/s.  To 

bring  it  in  any  closer  at  tnese  points  requires  a  much  more 

complicated  compensator  which  might  not  reduce  the  3  dB 

bandwidth  much  from  its  current  13  rad/s  position,  and  only 

gives  slightly  less  high  frequency  response  along  with  a 

higher  resonant  peak  around  2  r ad/s.  Figure  IV. 6-7  is  a  Bode 

plot  of  L  (phase  suppressed).  For  this  L  the  required 
lo  lo 

compensator  g  is 
1 

5 

g  =  2. 37  (Id  )ts  +  0.3d51Hs»0.6067  +  jJ.  35  76) 

i  (s+0.2)is+2dd)is+6d0'(s+0.6124) 

(s+11.36)(s+3j.3d)(s*1.979)  (4.6-17) 

*  l$  +  i.iJj)lS  +  /'./23)ls  +  22.25y 


Figure  IV. 6-3  contains  the  frequency  response  of  g  ,  and 

1 

snows  that  it  looks  line  a  lead-lag  filter.  Looking  a:  a  root 
locus  of  the  system,  the  for  high  compensator  gain  values  tr.e 
unstable  open-loop  pole  migrates  across  the  j  axis  to  an 
open-loop  zero.  This  again  illustrates  why  a  non-minimum 
phase  zero  cannot  be  used  here.  This  zero  is  unaffected  by 
feedback,  and  under  high  gain  conditions  an  open-loop  pole 
-oi  grates  to  that  zero,  making  the  system  unstable. 


over  7  t  n  order.  From  Figure  IV. 5-3  it  seems  as  if 


NOMINAL  LOOP  ONE  T RfiNSM I S 5  I  ON 


Fig. I  V . 6- 7 :  Log  Magnitude  of  L  Verses  Frequency 

lo 


tne  compensator  can  be  simplified.  A  transfer  function  which 

mimics  the  snaoe  of  g  is: 

1 

6 

g  =  1.9 2  u  C  1 J  j  y  s  +  6 )  (4.6-19) 

1  ( s  +  100)  ( s  +  oUd) 


Inis  coitipe  n  sa  tor  leads  to  a  neater  bandwidth  since  it  tends 
tu  nave  a  greater  gain  at  each  frequency  tnan  the  exact 
compensator.  A  comparison  of  the  two  AOA  channel  compensators 
is  in  Figure  IV. 6-10  .  ^hen  the  simplified  compensator  is 
multiplied  by  the  nominal  plant  and  the  loop  closed,  the 
resulting  bandwidth  is  13  rad/s,  higher  as  classical  control 
theory  predicts.  The  simpler  compensation  is  used  in  the 
simulation  to  keep  the  order  of  the  whole  system  at  a  minimum 
to  reduce  simulation  costs  and  computer  loading. 
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Firj.  IV.6-10:  Comparison  of  Full  and  Reduced  Order 
Loop  One  Compensators 


I V .  7  Design  jf  Second  Loop  Compensation 

Once  tne  first  loop  is  desiyned,  the  knowledge  of  tnat 

first  loop  compensation  is  used  to  give  more  exact(more 

optimal  in  tne  sense  of  minimizing  oanowidtn)  compensa ti on 

f or  tne  seconc  loopLl8j.  Tne  exact  equations  for  t  and  t ^ 

are  suosituted  into  tne  equations  for  t  and  t  to  give 

£1  22 


where 


t  »  f  L  /(I  +  L  ) 
22  22  2 e  2e 


L  =  g  Q 
2e  2  22e 


(4.7-1) 


0  (1  +  L  )/(l  +  L  -  y  ) 


(4.7-2) 


*r  s  Q  Q  /( Q  Q  ) 
11  22  12  21 


(4.7-3) 


The  Q  calculated  using  the  first  loop  compensation  are 
22e 

in  Tables  IV. 7-1  tnrough  IV. 7-3  ,  differing  by  flight 


condition.  At  low  frequencies  where 
approximated  Dy 


this  is 


t  =  f  (4.7-4) 

22  22 


When  L  I  «:  1  at  higher  frequencies  the  equation  can  be 
I  2el 

represented  by 

t  -  f  L  (4.7-5) 

22  22  2e 


At  frequencies  in  between,  the  simplifications  are  not  valid 

and  the  true  relationships  are  used.  The  disturbance  output 

t  is  now  given  by: 

21 

t  =  f  L  /lQ  (1  +  L  )(1  +  L  -  7  )J  (4.7-6) 

21  11  1  21  2e  1 


Again  at  <o  sucn  tnat  J  L  |  »  1,  the  response  is 

t  =f  /lQ  (1+L  )  (4.7-7) 

21  11  21  2e 


And  at  hign  frequencies  wnere 


t  =  f  L  / CQ  (1  ♦  L  )(1  -  Y  )] 


(4.7-8) 


TABLE  IV.  7-1 

Equivalent  Q  Plant  Transfer  Functions  for 


Fliqnt  C o n d i t i o n  One 


Ho  Failures: 

l*i.jlS+J.ju97±jU.4Doi)is+13.24)($+loU)(s+59.67 


Stabilator  Failed: 

21«).7(s  +  0.337ij0.4401)(s  +  53.77)(sn3.24)(sna0)(s  +  60.93) 


S  +  l  .  34 j 


(,  s  +  0  .  Jld3± j U.257)(s  +  io.6±jo.57)(s  +  12o±j95.1; (s-0 . i u46 ) { s  + 1  .  5  46  } 


vS  +  U.U14o/+jU.21ai)^s  +  ll3.d+jDO.JdHS  +  1.9/9){s-J.3d51) 
Cdridra  F d  i  lea: 


(s+0.0147+j0 

.2ly)(s+lll+j209)(! 

mtp.uEBiiti:niRHR»nfa 

Ailerons  Fai 

224 . 3 ( $  +  0 . 6 

led: 

696-i-jO.  2279)  (  s  +  130 

)  ( s  +  23  .  7  7  ) ( s  +  45 . 1 1 ) ( ■ 

s+13.24 

(s+0.3147+jO 

. 2 1 3 ) ( s  +  20 . 5+ j  2.33 

)Ts+i^r+ji20) 

(T-0.351 

rnr+rr 

Stadilator  Failed: 

134. 5( s+0. 6124) (s  +  1 . 153)(s  +  7.  : 

723  )  (s  +  13.24) 

■  ( 5  +  130) 

( s  +  53. 3 

Staoildtor  and  Ailerons  Failed: 

73*3^3  (  s  +  O.  3  J34+ j  0 . 3  0  C5  3  )  (  s  + 1  b  .  3  9 + j  3  . 1 3 1 )  (  s  + 1 3  0  )  (  s  +  4  6 . 3  9 


s+ 1 1 5+ 


TABLE  IV. 7-1 


Equations  4.7-1  through  4.7-8  are  used  to  find  both  the 

tracking  and  disturbance  bounds  for  the  second  loop.  Notice 

that  t  contains  only  one  term  due  to  tracking  response.  The 
22 

disturbance  rejection  has  been  "built  in"  by  the  time  the 
second  loop  is  designed. 

The  dominant  bounds  for  the  second  loop  are  due  to 

disturbance  rejection  at  the  low  frequencies  and  tracking 

response  at  hign  frequencies.  This  differs  from  the  first 

loop  wnere  tne  bounds  are  derived  from  tracking  requirements 

only.  Again  this  is  aue  to  tne  same  reason,  tne  derivation  of 

Q  from  P_‘  that  madejQ  |>|  Q  |also  madejQ  |>  |Q  J.  If 

I  L  >>  1  then  t  is  given  by 
1  1  1  21 


t  =  f 

/[Q 

(1  - 

►  L  ) 

(4.7-9) 

21 

11 

21 

2e 

which 

can  be  rearranged  to 

the  form 

1/(1  + 

L  ) 

< 

b  Q  /f 

(4.7-10) 

2e 

21  21  11 

which 

for  1  L  1 

>>  1 

is  approximated 

as 

'  2e 

L  1  > 

1  f 

/(b  i 

3  )  1 
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is  simply 

referred 

to  as  L 

2e 

2 

Tne  plant  templates  for  Q  are  in  Figure  IV. 7-1  .  Notice 

22e 

that  tne  uncertainty  is  increased  from  the  first  loop(Figure 


IV. 6-5).  As  was  done  on  the  first  loop,  the  lower  tracking 


Fit].  IV.  7-1:  Q  Templates 
2  2eq 

bound  is  decreased  at  nijner  frequencies  to  provide 

bandwidth  reduction.  Tnis  is  seen  as  an  increase  in  A  T  (  j  w  ) 

wnicn  t.ne  te up  lutes  lost  not  violate.  Tne  tracKinj  response 

cun  oe  restored  oy  usin.j  t.ie  actual  oounu  values  when 

desitjninj  f  .  Tne  original  ana  modified  A  T  ( j  u»  )  bounds 
■  .  ^2 

are  listed  in  Table  IV. 7-4. 
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Modified 
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Figure  IV. 7-2  shows  the  bounds  for  both  tracking  and 

disturbance  rejection  for  loop  two.  Up  to  20  rad/s  disturbance 

rejection  dominates,  with  the  high  frequency  bounds  determined 

both  by  tracking  and  disturbance  considerations.  The  UHFB  is 

a  26  dB  vertical  line.  As  before  L  contains  an  unstable  root 

2 

and  must  start  off  at  -130  degrees  for  stability.  With  the 

nominal  plant  again  FC2:Canard  Failed,  an  L  that  satisfies 

2o 

the  loop  bounds  witn  minimum  bandwidth  is: 


2  ,  732  (  10  ) ( s  +  0. 3s51  )  ( s+25) 
s  +  1)  (s-0.3351  J  (  s  +  300)  ( s+900 . 


(4.7-12) 


L  is  shown  in  both  the  Nicnol's  chart  in  Figure  IV. 7-3  and 
2o 

the  frequency  response  in  Figure  IV. 7-4  .  Again  some 

overdesign  of  L  is  evident  since  the  bounds  are  not  "just 
2o 

met"  at  each  specific  frequency.  As  before,  this  is  done  to 


Fi  j  .  I V . 7-2 :  3ounds  on  L 

2 

so:newridt  greater  tnan  L  at  y  =  26  rad/s.  i  t  n  tne  nominal 

1  o  o 

condition  FC2 : Staoi 1  a  tor  Failed,  tne  required  loop 
c o.npen sa  t  i  on  is: 

j  =  333ouis+j1jis+2o./1)(s+17.26)is+1.973j 

2  U  +  l)ls  +  d;(s  +  doJHs  +  ^jj)is  +  o./l^4  +  jo.2Ulyj 

(s  +  d.3d5)(s  +  25)(s  +  lll+j20oHs  +  d.ul4o  +  ju.2l3)  (4.7-13) 

*  ( s  +  2 9. 9  + j 1 9. 03 )( s  + 2  32. 6  )( s  + 13. 24  )(s  +  l&u) 

Tne  frequency  response  of  g  is  in  Figure  IV. 7-5  .  Notice 

2 


J 


20.72(s+a.2)(s+0.4) 


(4.7-14) 
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design  tecnn i Hues .  To  i  i  :.i  i  t  effects  from  saturations  of 
internal  elements,  feedback  patrts  around  tne  saturating 
elements  are  required.  QFT  has  been  extended  to  nandle 
feedback  of  internal  variables,  out  that  is  beyond  tne  thrust 
of  this  thesis[24].  After  the  system  is  designed,  it  is 
noticed  that  if  L  is  shaped  a  bit  differently  then  the  gain 
can  be  saved  at  both  low  and  high  frequencies.  This  in  turn 
reduces  the  closed-loop  bandwidth.  Unfortunately,  this  cannot 
be  included  in  this  chapter  due  to  printing  deadlines; 
however,  the  improved  loops  and  their  compensation  are  in 
Appendix  H  . 

Tne  original  plan  to  control  the  flignt  patn  angle  and 
velocity  in  tne  SToL  landing  scenario  nad  to  be  modified  to 
controlling  the  angle  of  attack  ana  velocity  in  order  to  get 
iiiini.nu. ii  pnase  plants.  QFT  can  be  used  in  tne  control  of  non- 
minimum  pnase  plants,  out  restrictions  placed  on  the 
bandwidth  of  the  loop  transmission  make  it  impractical  to 
control  the  type  of  plant  in  this  thesis  containing  both  open 
loop  RHP  poles  and  zeros.  Just  controlling  AOA  implies  that 
landing  control  is  lost  since  the  exact  pitch  angle  is 
unknown  .  Thus,  this  control  scheme  cannot  be  used  to 
effectively  control  landing  since  direct  control  of  the 
decent  is  lost;  however,  it  can  be  used  to  insure  stability 
in  tne  event  of  control  surface  failures.  Tnere  is  a  chance, 
not  investigated,  tnat  oy  expanding  uet|£  }  vectors  can  oe 
founu  tnat  force  tne  plant  into  minimum  phase  form.  Tnis  is 
stated  a^ain  in  Cnapter  o. 

before  designing  a  ■■ll.'IO  control  scneme  using  QFT  one 


snould  first  form  trie  jetermi nant  arid  jeter,  nine  tne  range  of 
weightings  which  result  in  mi  ni.nu.,i  pnase  plant  Matrices.  Only 
then  may  QFT  continue.  The  reversal  of  diagonal  dominance 
seen  in  the  equivalent  plants  could  be  caused  by  forcing  the 
plant  to  be  minimum  phase,  and  does  not  seem  to  result  from 
QFT  directly. 

Restating,  using  QFT,  compensators  are  found  wnich 

enables  the  system  to  meet  the  given  performance  tolerances. 

Tnese  higher  order  compensators  are  approximated  by  lower 

order  compensators  without  significant  error.  Loop 

transmission  oandwidths  are  IZ  and  do  rad/sec  for  the  AOA  and 

velocity  cnannel,  respectively.  Since  development  of  these 

compensators ,  reshaping  of  tne  loop  transmissions  nave 

significantly  reuucea  tne  .magnitude  of  g  and  g  as  a 

1  Z 

function  of  frequency.  These  improved  compensators  are  in 
Appendix  H. 


V .  Simulation  Results 


V . 1  Introduction 

This  chapter  presents  the  results  of  computer 
simulations  of  the  compensated  aircraft  under  normal  and 
failure  conditions.  Nineteen  separate  cases  over  three  flight 
conditions  are  examined.  The  commanded  outputs  are  evaluated 
to  see  if  performance  tolerances  are  met.  The  control 
surfaces  are  examined  to  see  if  any  physical  limits  not 
included  in  the  controller  synthesis  have  been  exceeded,  or 
saturated.  Using  these  results,  the  practical  application  of 
tne  controller  to  an  aircraft  is  examined,  and  any  possible 
corrections  or  improvements  are  suggested. 

V . 2  Computer  Model 

The  aircraft  is  simulated  using  the  CAD  package 
MATRIX  X  C  25 ]—  This  program  allows  the  system  as  a  whole  to 
be  simulated  in  block  form  using  transfer  functions  rather 
than  transforming  the  entire  system  into  state  space. 
Individual  control  surface  failures  are  simulated  by  removing 
connections  between  the  specific  blocks.  Appendix  G  goes  into 
further  detail  on  the  actual  mechanics  of  the  computer 
simulation. 

Before  simulation  can  occur  the  system  must  be  modified 
to  reflect  tne  changes  done  during  loop  compensation  design. 
Tne  original  signal  flow  cnart  for  the  control  system  is  in 
Figure  V.2-1  .  Since  the  columns  in  _P'  are  switched  this 


grapn  is  no  longer  valid.  The  new  system  in  Figure  V.2-2 
reflects  the  changes  accomplished  to  compensate  the  system. 
The  simulation  uses  the  lower  bandwidth  compensators 
developed  in  Appendix  H  rather  than  the  ones  developed  in 
Chapter  IV.  The  first  computer  runs  are  made  without  taking 
into  account  surface  rates  or  deflections.  This  is  to  verify 


Fig.V.2-1:  Signal  Flow  Graph  of  Original  Plant 

that  the  compensators  designed  using  QFT  did  indeed  force  the 

plant  to  meet  the  performance  tolerances.  After  this  MATRIX  X 

allows  the  designer  to  include  saturations  directly  in  tne 

plant  model.  Botn  surface  rate  and  deflection  limits  are 
* 

introduced.  The  aerodynamic  surfaces  are  limited  to  maximum 
aeflections  of  ±20  ueyrees  at  rates  up  to  30  deg/sec.  Of 


course  eucn  type  of  surface  nas  its  own  individual  limits, 
out  just  one  set  of  limits  for  all  surfaces  is  chosen  to 
simplify  the  simulation.  Tne  limits  on  tne  vanes  have  a 
faster  deflection  rate  of  30  deg/sec,  out  since  the  vanes  are 
very  non-linear  over  large  deflections  within  the  range  of 
their  stops,  their  deflection  is  limited  to  ±8  degrees  even 
though  they  can  actually  deflect  much  farther  . 


V . 3  Simul ati on  Resul ts 

Figures  V.3-la  tnrough  V . 3- 1 d  show  the  simulation 
results  for  all  nineteen  flight  condi ti on/f a i 1 ure 
combinations.  Step  inputs  for  both  velocity  and  angle  of 
attack  are  used.  Figure  V.3-la  is  the  AOA  tracking  response 
to  an  inch  deflection  stick  input.  Without  error,  an  inch  of 
stick  deflection  should  result  in  a  change  of  AOA  by  3.5 
degrees.  QFT  guarantees  that  the  response  is  within  the  given 
response  tolerances.  The  figure  shows  that  all  the  transient 
parts  of  the  responses  except  one  lie  within  the  bounds. 
Possible  cause  for  the  response  lying  just  outside  is  the 
simplified  compensators  used.  That  response  is  underdamped 
and  corresponds  to  a  plant  on  the  upper  right  hand  side  of 
tne  plant  template.  The  simplified  compensation  raises  the 
frequency  response  of  tne  tracking  output  enougn  to  slightly 
exceed  tne  upper  tracking  Dound.  The  increased  overshoot  is 
not  very  great,  on  the  order  of  2  percent,  whicn  should  be 
toleraole.  Since  the  system  is  Type  0,  some  steady  state 
error  is  expected.  The  worst  steady  state  tracking  error  is 
about  7  percent.  The  Figures  of  Merit  for  the  AOA  command 
input  are  contained  in  Table  V . 3 - 1  and  lie  within  the 
specified  ranges(see  pages  61+  and  68  f°r  The  ranges  of 
the  Figures  of  Merit  for  the  two  channels)  except  for 
FClrAilerons  Failed  case  which  shows  the  high  overshoot  noted 
previously.  The  disturbance  output  on  the  velocity  channel 
for  the  step  AOA  command  is  in  Figure  V . 3 - 1 b .  The  largest 
deviation  is  -1  ft/sec  with  a  steady  state  value  of  -0.9  . 
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This  corresponds  to  a  DC  gain  of  12  dB  down  from  the  tracking 

response,  within  the  -10  dB  limit.  This  response  is  for  the 

nominal  plant  P  which  is  on  the  bottom  of  the  plant  template 
o 

for  most  frequencies,  and  for  which  the  DC  gain  and  open  loop 

bandwidth  is  at  a  minimum.  Most  of  the  plants  exhibit  much 

smaller  variations  which  indicates  greater  loop  gain  and 

Da nd width  leading  to  better  disturoance  rejection.. 

Tne  tracking  response  for  a  velocity  input  is  shown  in 

Figure  V . 3-lc.  A  +5  ft/sec  velocity  change  is  commanded.  All 

of  tne  responses  are  initially  within  the  previously 

described  bounds,  but  at  some  time  deviate  since  the  steady 

state  error  is  non-zero.  The  worst  error  is  again  for  the 

nominal  plant,  FC2:  Canard  Failed  with  the  lowest  dc  open- 

loop  gain.  The  responses  show  the  first  order  characteristics 

desired  in  the  original  specifications.  Table  V. 3-2  contains 

the  velocity  response  Figures  of  Merit,  showing  that  the 

bounds  are  met  for  all  cases.  The  AOA  disturbance  for  the 

step  ve 1 oc i ty- comma nd  is  in  Figure  V  .  3 -  1  d .  The  'worst'  curve 

has  a  steady  state  value  of  17  dB  below  the  command  input  , 

well  witnin  tne  tolerances.  This  indicates  that  the  loop 

transmission  is  overaesigned  and  can  be  decreased  without 

violating  tne  disturbance  tolerances. 

As  a  cnee*  of  tne  QFT  design  proceedure,  tne  system  is 

simulated  for  FC2:Canard,  Ailerons  Failed  and  FC3:Canard, 

Ailerons  Failed.  Both  of  these  plants  nave  non-minimum  phase 
✓ 

Q  and  should  be  unstable  under  high  forward  gain  and  negative 
1 J 

feedback.  Figures  V.3-2a-d  do  indeed  show  the  projected 
instability  for  both  AOA  and  velocity  step  command  inputs. 


TABLE  V.3-2 


o~  a 


The  system  has  acted  as  predicted  by  QFT(and  classical 
control  theory). 


In  summary,  tne  linear  simulation  snows  that  tne  control 
system  satisfies  tne  specifications  over  tne  range  of 
uncertainty  for  wnicn  tne  design  is  made(is  rooust  for  tne 
aiven  control  surface  failures  at  tne  tnree  flight 
conditions;.  This  is  guaranteed  wnen  QFT  is  applied 
correctly.  What  is  not  shown  are  the  control  surface 
deflections  and  rates.  For  a  more  realistic  simulation  these 
must  be  included.  This  is  done  in  the  next  section. 


V . 4  System  Simulation  With  Non-Linearities  Included 

The  previous  section  shows  that  the  system  meets  the 
desired  performance  tolerances  for  a  wide  variety  of  failure 
conditions  a  summing  a  linear  plant.  What  is  not  shown  are  the. 
internal  variaoles  of  surface  deflection  and  deflection  rates 
wnicn  nave  definite  limits.  These  internal  variaoles  must  oe 
exduiineu  to  determine  whetner  or  not  tne  controller  works  'in 
tne  real  world'.  To  do  tnis  the  limits  mentioned  in  Section 
V.3  are  included  directly  in  the  simulation.  The  figures  in 
this  section  contain  plots  of  control  surface  movement.  The 
two  surfaces  shown  are  the  canard  and  top  reversing  vane. 
Given  these  responses,  the  response  for  the  rest  of  the 
surfaces  can  be  calculated  using  the  A  vector.  The  other 
surfaces  deflect  the  amount  of  degrees  indicated  on  the  plot 
times  the  A  that  premultiplies  the  particular  surface.  For 
instance,  the  stabilator  deflection  is  -1  times  the  canard 

113 


deflection.  The  same  procedure  applies  to  calculate  the 
bottom  vane  deflection  from  the  top  vane  plot,  the  top  vane 
deflection  is  multiplied  by  -0.6  .  This  is  a  bit  misleading 
since  not  all  surfaces  saturate  at  the  same  deflection.  Thus 
it  gives  the  impression  that  all  surfaces  have  saturated  when 
tnere  is  actually  some  control  authority  left.  The 
approximation  is  made  to  simplify  the  computer  simulation. 

Figures  V . 4- la-d  are  graphs  of  the  system  response  to  a 
3.5  degree  AOA  command.  Figure  V.4-la  is  the  AOA  output  and 
snows  that  most  of  the  responses  do  approach  3.5  degrees,  but 
a  few,  connected  with  single  failure  cases  at  each  flight 
condition,  have  significant  steady  state  error.  The  case 
FC2:Canard  failed  has  the  worst  steady  state  error,  over  20 
percent.  The  reason  these  responses  fail  to  meet  the 
specifications  could  stem  from  the  control  surfaces 
saturating,  failing  to  gi.ve  the  extra  amount  of  control 
requires  to  reach  the  desired  value.  Possibly  simplifications 
made  adding  tne  saturation  blocks,  or  in  other  places  in  the 
simulation,  uroppeu  tne  loop  transmission  enougn  tnat  with 
tnese  severe  failures  tne  dc  gain  dropped  enough  to  increase 
steady  state  error.  The  system  is  unstable  for  all  double  . 
failure  casesvwnicn  makes  sense  since  when  all  surfaces 
saturate  tne  system  becomes  essentially  'open-loop',  and  the 
plant  is  open-loop  unstable).  This  instaoility  can  be  seen  in 
Figure  V-.4-2  which  plots  the  AOA  output  for  the  unstable 
cases.  For  the  rest  of  the  cases  the  output  is  acceptable, 
with  the  Figures  of  Merit  in  Table  V . 4 - 1 .  Single  failure 
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F i g . V . 4 - 2 :  Unstable  AOA  Output  for  Step  AOA  input  with 
Plants  Containing  Double  Surface  Failures 

cases  also  increase  the  disturbance  level.  Note  that  in 

Figure  V .  4  - 1  c  the  velocity  drops  as  .nuch  as  3.5  ft/sec,  the 

sa;ne  amount  as  commanaed(AOA)  ,  or  0  dB  disturbance  rejection 

This  can  oe  disasterous,  especially  at  V  ,  since 

ini  n 

dece 1 1 e ra t i on  at  tnat  flignt  condition  causes  a  stall.  The 
rest  of  tne  cases  experience  cnanges  in  velocity  within  tne 
10  aB  tolerance.  Figures  V .  4 - 1  b  and  V.4-ld  are  the 
deflections  of  the  aerodynamic  surfaces  and  thrust  reversing 
vanes  respectively.  The  surfaces  'nit  the  s tops  1 ( def 1 ect i on 
saturation)  for  only  the  severe  douDle  failure  conditions, 
causing  instability;  however,  other  cases  undergo  rate 
saturation.  This  rate  saturation  does  not  cause  instability, 
possibly  due  to  the  small  magnitude  of  the  unstable  root 


ro  kO  ro 


(around  .3)  allowing  tiia  surface  to  deflect  to  the  proper 

position  for  countering  the  unstable  moments  before  they  'get 
out  of  hand'.  However,  even  though  the  surfaces  do  not 
saturate  for  single  failures  in  this  simulation,  the  trim 
position  of  the  surfaces  has  not  been  accounted  for.  Thus 
some  of  the  previously  non-saturated  cases  may  actually  be  in 
saturation  due  to  this.  The  reversing  vanes  show  the  same 
pattern  of  saturation,  except  tney  also  saturate  for 
FC2:Bottom  Vanes  failed. 

The  response  of  the  system  to  a  5  ft/sec  step  velocity 
command  is  shown  in  Figures  V.4-3a-d.  The  velocity  response 
response  in  Figure  V.4-3a  shows  that  the  velocity  follows  the 
input  with  only  a  slight  increase  in  steady-state  error 
except  for  the  multiple  failure  cases.  The  Figures  of  Merit 
for  the  velocity  command  are  in  Table  V. 4-2.  Multiple 
failures  again  cause  instability.  For  some  single  failures 
the  final  velocity  value  is  significantly  less  than 
desired.  The  loss  of  the  control  surfaces  decreased  the 
control  authority  to  the  point  where  the  desired  output  is 
pnysically  imposs i o 1 e ( for  that  particular  configuration).  In 
Figure  V.4-3c  tne  aisturbance  rejection  also  stays  below  the 
10  dB  level(less  than  1.6  degrees  in  magnitude)  for  most 
single  failure  cases  and  douole  failures(for  which  it  is 
unstaole).  In  the  worse  staole  case,  FC 3 : S ta b i 1  a  to r  Canard 
Failed,  the  0  dB  level  is  approacned  with  the  angle  of  attack; 
decreasing  5  degrees  for  a  5  ft/sec  velocity  command  input. 
From  Figures  V.4-3c,d  the  surface  deflections  saturate 


•  4  -  3  a  :  Velocity  Response  Fig,V.4-3b:  Aerodynamic  Surface 
or  Velocity  Step  Input  Def 1 ecti on  for  Velocity  Step 
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totally  only  for  double  failures.  The  vanes  saturate  for  all 
double  surface  failure  conditions. 

The  simulations  including  saturation  effects  show  that 
the  system  remains  stable  and  within  performance  tolerances 
for  most  single  failure  cases,  but  double  failures  lead  to 
severe  saturation  and  instability  that  certainly  keeps  the 
system  from  meeting  the  specifications. 


V  .  b  Sum.iia  r.y 

For  a  linear  system,  the  compensators  developed  using 
gFT  meet  the  perf ormance  tolerances  for  tne  flight  conditions 
and  surface  failures  which  they  are  designed  for.  This  is  not: 
surprising  since  QFT  guarantees  the  desired  performance  in 
the  linear  ( nonsaturated)  case.  When  saturations  of  control 
surfaces  are  included,  the  tolerances  are  met  for  the  healthy 
aircraft  and  most  single  control  surface  failures.  Double 
failures  lead  to  instability  due  to  the  lack  of  available 
control  authority.  QFT  has  been  extended  to  include,  and 
compensate  for,  saturated  elements  in  a  SISO  system,  but  has 


not  yet  been  extended  to  the  MIMO  case. 


In  other  words, 


QFT  provides  a  robust  controller  for  this  aircraft  subjected 
to  surface  failures  and  changing  flight  conditions. 


V I .  Conclusions  and  Recommendations 

V 1 . 1  Discussion 

This  thesis  has  demonstrated  the  application  of 
Quantitative  Feedback  Theory  in  designing  control  laws  for  a 
reconf igurable  flight  control  system  aircraft.  Fixed 
compensation  provides  robust  control  for  three  different 
flight  conditions  with  numerous  control  surface  failures.  The 
use  of  fixed  compensation  is  important  since  it  implies  that 
identification  of  failures  and/or  scheduling  for  change  of 
flight  condition  can  be  minimized,  reducing  i dent i f i ca t i on 
failures,  false  alarms,  required  memory,  and  time  delay.  This 
compensation  results  from  application  of  QFT(minimum  phase 
technique)  to  the  linearized,  small  perturbation  equations  of 
motion  for  an  open-loop  unstable  a i rcra f t ( STOL  F - 1 5  ) .  The  use 
of  QFT  is  simple  and  strai ght f orwa rd  once  the  plants  are 
described  prpperly,  i.e.  non-minimum  phase  terms  are 
eliminated.  In  order  to  eliminate  non-minimum  phase  terms 
proper  control  surface  weightings,  derived  from  an  equation 
for  the  determinant  of  the  plant,  are  used.  QFT  affords  the 
designer  insights  on  how  wide  the  loop  bandwidth  should  be’ 
for  various  failures  and  disturoance  rejection  levels.  QFT 
enables  the  designer  to  eliminate  those  conditions  which 
cause  an  unrealistically  nigh  loop  bandwidth.  Finally,  QFT 
does  the^above  from  initial  stages  of  the  design 
ef f ort ( ra ther  than  finding  out  it  "doesn't  work"  at  the  end 
of  the  design  process). 


VI . 2  Conclusions 

1.  Quantitative  Feedback  Theory  provides  robust  control 
for  the  complicated  aircraft  in  this  thesis.  Fixed  control 
laws  for  both  healthy  and  damaged  aircraft  are  designed  such 
that  performance  is  within  the  established  tolerances  using 
linear  model  simulation. 

2.  The  original  idea  of  controlling  the  flight  path 
angle  and  velocity  is  modified  to  controlling  the  angle  of 
attacK  and  the  velocity.  Tnis  is  due  to  the  non-minimum  phase  . 
plants  associated  witn  tne  flight  path  angle  which  limit  the 
usefulness  of  QFT  in  this  instance. 

3.  Control  surfaces  must  be  combined  in  such  a  fashion 
that  the  resulting  plant  matrices  are  minimum  phase.  This 
should  be  the  first  design  step  for  problems  having  MIMO 
plants  that  are  open-loop  unstable,  or  exhibit  non-minimum 
phase  terms.  Doing  so  ensures  that  the  problem  formulation  is 
one  that  QFT  can  'handle'  without  using  the  more  involved 
"Singular  G"  -method.  Choosing  the  proper  weightings  requires 
examination  of  the  influence  of  the  individual  surface 
weightings  on  the  plant  determinant. 

4.  CornD  i  nation  of  surfaces  using  weignting  factors 
designed  to  drive  the  plant  mini  mum  pnase  tended  to  reduce 
tne  control  surfdce  influence  on  some  output  variaoles  and 
increase  tne  effectiveness  of  others.  In  fact  the 
effectiveness  is  reduced  so  much  that  the  diagonal  dominance 
reverses.  The  necessary  process  of  eliminating  RHP  zeros 
appears  to  cause  the  reversal  effect  for  this  particular  plant. 
Previous  theses  have  not  noticed  this  effect,  possibly  for  two 


reasons:  one,  the  problems  may  be  related  to  tne  particular 
plant,  or  two,  the  plants  used  in  QFT  theses  up  to  this  time 
have  been  open-loop  stable  (except  for  one  channel  of  Arnold 
and  Walkes  use  of  Singular  G  on  the  X-29[-J])  and  the  effect 
only  shows  up  for  open-loop  unstable  plants.  The  opinion  of 
this  thesis  is  that  the  particular  plant  data  is  to  blame, 
not  the  technique. 

5.  When  control  surface  rates  and  deflection  limits  are 
added  to  tne  simulation  the  surfaces  saturate  only  for  the 
double  failure  cases.  Saturation  in  these  cases  leads  to  loss 
of  aircraft  control.  This  does  not  take  into  account  any  trim 
position  tne  surface  may  oe  at  before  the  command,  thus 
saturation  could  occur  for  even  lesser  failure  cases.  In  all 
single  failure  cases  the  system  remains  stable;  however,  for 
some  instances  the  steady  state  tracking  and  disturbance 
rejection  tolerances  are  not  reached.  One  possible 
explaination  of  these  effects  could  be  the  reduction  of 
control  surface  effectiveness  resulting  from  the  desire  to 
maintain  minimum  phase  plants.  Enough  effectiveness  seems  to 
be  reduced  so  that  under  double  surface  failures  the  aircraft 
does  not  have  enough  control  authority  left  to  remain  stable 
while  single  surface  failure  reduces  tracking  response  and 
increases  disturbance  output. 

o.  Reversing  vane  failure  can  only  oe  tolerated  at  one 
flignt  condition,  and  that  is  only  for  tne  oottoin  vanes.  For 
otner  conditions  non-minimum  phase  plants  resulted.  This 
indicates  that  failure  of  these  surfaces  could  cause 


instability.  Anotner  failure  c a se ( Ca na rds  ,  A i 1 e rons  Failed) 
that  QFT  indicated  would  be  unstable  is  checked  and  found 
unstable.  This  provides  a  check  of  QFT's  ability  to  forcast 
cases  it  cannot  compensate  for. 

7.  Digitizing  effects  are  not  included  in  the  loop 
transmission  design.  As  with  all  future  aircraft  the 
STOL  F - 1 5  aircraft  will  have  a  digital  flight  control 
system.  The  digitizing  effects  introduce  effective  lags(as 
seen  from  the  analog  design),  finite  word  lengtns,  round-off 
errors,  and  maximum  sampling  rates.  When  using  analog  design 
proceedures  for  developing  an  eventual  digital  flight  control 
system  tne  designer  might  want  to  consider  adding  extra  lead 
to  tne  system  to  ensure  stability  wnen  impl imented.  Dr. 
Horowitz  has  recently  extended  QFT  into  the  discrete 
domain [26  j. 

3.  Loop  compensation  for  both  channels  is  high  order, 
being  bth  over  7th  order  in  one  case  and  10th  over  11th  in 
the  other.  Approximate  compensators  with  reduced  order  were 
found  that  mimiced  the  frequency  response  of  the  full  order 
compensators.  The  full  and  reduced  order  compensators 
designed  included  considerable  overdesign.  Compensator 
designs  with  less  bandwidth  are  included  in  Appendix  H,  and 
are  used  in  the  simulation. 

9.  The  design  of  compensators  using  QFT  would  be  much 
less  time  consuming  if  a  computer  aided  design(CAD)  package 
is  implemented.  Witn  sucn  a  program  tne  design  of 
compensati on  for  nigner  oroer  pi  ants ( fourth  on  up)  will  be 
possible  witnin  tne  time  allotted  for  a  Masters  Thesis.  Doin' 


so  allows  more  "realistic"  design  problems  to  be  considered. 
For  example,  the  STOL  landing  proolem(as  with  any  aircraft) 
is  actually  a  six  degree-of-freedom  problem.  Using  a  CAD 
package  such  thesis  topics  may  become  practical  while  using 
QFT.  Current  AFIT  efforts  include  such  a  package,  and  these 
should  be  supported[  27] . 

V  I . 3  Recommendations 

1.  The  use  of  QFT  to  design  robust  reconf i gurabl e  flight 
control  systems  should  be  continued.  Future  proolems  should 
attempt  the  design  of  compensa t i on  for  systems  of  third  and 
riigner  orders,  or  as  in  this  thesis,  include  more  non- 
conventional  control  surfaces.  In  particular,  redesign  of  tne 
STOL  F - 1 5  aircraft  for  control  of  flight  path  angle,  side 
velocity,  forward  velocity,  and  roll  angle  to  provide 
increased  landing  control  would  be  an  useful  extension  of  QFT 
design  effort.  Another  extension  of  this  thesis  would  be  to 
retain  the  2X2  equivalent  plant  matrix,  but  include  control 
devices  on  the  aircraft  but  omitted  here  such  as  flaps, 
thrust  vectoring  vanes,  and  engine  throttle.  This  plant  might 
also  have  more  general  frequency  dependent  A  terms  for 
each  surface,  causing  the  surfaces  to  work  together  better 
tnan  simple  gains  terms  can  accomplish.  Including  more 
control  surfaces  has  tne  possibility  of  reducing  or 
eliminating  tne  proolems  wi to  non-minimum  phase  plant 
matrices  encountered  in  tnis  thesis.  This  in  turn  could  allow 


the  application  of  reconf  i  gu rao i 1 i ty  tneory  in  a  greater 
extent  than  done  here.  A  signal  flow  chart  of  this  system  is 
in  Figure  VI. 3-1  .  The  use  of  internal  loops  to  decrease  the 

effects  of  saturation  elements  could  also  be  investigated. 


Non-linear  simulation  of  future  designs  -using  STOL  F-15 
models  is  available  and  should  be  used  to  validate  the  linear 
designs  in  the  vicinity  of  the  appropriate  flignt  condition's 
equilibrium  point.  This  will  provide  a  much  greater  check  of 
the  control  system's  robustness. 

2.  When  using  QFT  to  design  compensators  for  MIMO  systems 
containing  unstable  poles  and/or  non- mini  mum  phase  zeros  the 
first  design  step  should  oe  to  expand  out  the  expression  for 
the  plant  determinant  including  the  weighting  terms.  Then  the 


range  of  weighting  terms  leading  to  minimum  phase  plants  must 
be  establisned.  Then  working  within  this  range  concepts  such 
as  “control  authority"  and  "effective  use  of  surfaces"  can  be 
explored.  Following  this  procedure  could  save  the  designer 
time  and  effort  wasted  on  impractical  plant  configurations. 

3.  Digitizing  effects  should  be  included  in  future  QFT 
flight  control  system  designs  using  analog  techniques.  These 
could  be  as  simple  as  first  order  lags.  Current  analog 
designs  from  QFT  must  be  converted  into  difference  equations 
for  implementation  in  modern  flight  control  systems.  The 
sampling  and  processing  delays  within  the  flight  control 
system  snould  oe  included  in  the  plant  uncertainty  to 
guarantee  stability. 

4.  Future  control  system  designs  should  also  include  at 
least  a  simple  model  of  the  most  important  part  of  the  flight 
control  system,  the  numan  pilot.  Doing  so  will  partially 
account  for  another  source  of  time  lag  and  possible 
instability. 

5.  A  QFT  CAD  package  should  be  implemented  to  provide 
SISO  compensation  design  and  MI  MO  -system  reduction  into 
equivalent  SISO  plants.  Doing  so  will  result  in  a  substantjal 
decrease  in  time  required  for  compensation  design,  and  will 
allow  the  designer  t  contemplate  more  difficult  problems. 

b.  Extension  of  QFT  theory  to  the  discrete  domain  for 
direct  design  of  digital  compensators,  rattier  than  analog 
design  assuming  fast-enough -sampling-rates,  should  be 
investigated.  Future  FCS  will  be  digital,  thus  the  design 


STQL  Aircraft  State  Space  Equations 


The  following  taolcs  contain  the  aerodynamic  data  used 
to  ootain  the  state  space  equations  for  the  STOL  F - 1 5 . 

Trie  data  is  provided  Dy  McDonnell  Douglas  Corporation.  This 
data  is  preliminary,  cursory,  and  experimental  and  should  not 
he  looked  at  as  representing  any  finished  aircraft.  The 
stability  derivatives  are  originally  given  in  the  aircraft 
body  axis,  but  are  converted  to  the  stability  axis  for  FCS 
design.  Also  some  derivatives,  especially  those  related  to 
control  surface  deflections,  have  units  per  degree  and  are 
changed  to  a  per  radian  measure.  During  simulation  the  output 
is  converted  back  into  a  per  degree  measure. 


Non-Dimensional  Body  Axis  Stability  Der i va t i ves ( 1/deg ) 


CZA  = 

-0.0613756 

CZD3 

= 

-0.00036784 

CZQ  = 

0.0 

CZD4 

= 

-0.00036784 

CZH  = 

-0.000026124 

CZD5 

= 

0.0016781 

CZU  = 

0.011254 

CZD6 

-0.0016781 

CZD1  = 

-0.001452119 

CZD7 

= 

0.0016781 

CZD2  = 

-0.00549256 

CZD3 

= 

-0.0016781 

CMA  = 

0.012469 

CMD  3 

0.00096437 

CMQ  = 

-0 . 0985o 1 

CMD4 

= 

-0.0012491 

CMH  = 

-0.0000346518 

CMD  5 

= 

0.0014572 

CM  U  = 

-0.0115264 

CMD6 

= 

-0.0011263 

CMD 1  = 

0.0051492 

CMD7 

= 

0.0014572 

CM02- 

-0.0098354 

CMD8 

= 

-0.0011263 

CXA  = 

0 . 0025bl4 

CXD3 

= 

0.001395584 

CXQ  = 

0.0 

CXD4 

= 

0.001395584 

CXH  = 

0.00043821 

CXD5 

= 

-0.0045798 

CXU  = 

-0.173411 

CXD6 

= 

-0.0045798 

CXD1  = 

-0.00065298 

CXD7 

= 

-0.0045798 

CXD2  = 

-0.00099344 

CXD8 

= 

-0.0045798 

TABLE  A-2 


Aerodynamic  Data  for  Flight  Condition 


Non-Dimensional  Body  Axis 


CZA  =  -0.07062331 
CZQ  =  0.0 

CZH  =  -0.00001818665 
CZU  =  0.006534185 

C  ZD  1 =  -0.002571649 
CZD2  =  -0.009552322 

CM A  =  0.009313564 

CMQ  =  -0.169491 
CMH  =  -0.0000388546 
CMU  =  -0.0140683 
CMO 1 =  0.0052887 

C  M  D  2  =  -0.0107546 

CXA  =  0.00208763 

CXQ  =  0.0 

CXH  =  0.000631123 

CXU  =  -0.246589 
CXD1 =  -0.00153014 
CXD2=  -0.002016566 


Stability  Deri vati ves ( 1/deg) 


CZD3  =  -0.004384277 
CZD4  =  -0.004515593 
CZD5  =  0.00135028 

CZD6  =  -0.00135028 
CZD7  =  0.00135028 

CZD8  =  -0.00135028 

CMD3  =  0.00112899 

CM04  =  -0.00214211 
CMD5  =  0.00129075 

CMD6  =  -0.00137616 
CMD7  =  0.00129075 

CMD8  =  -0.00137616 

CXD3  =  0.001365414 
CXD4  =  0.000946531 
CXD5  =  -0.00340353 
CXD6  =  -0.00340353 
CXD7  =  -0.00340353 
CXD8  =  -0.00340353 


Non-Dimensional  Body  Axis  Stability  Deri  vati ves ( 1/deg ) 


CZA  = 

-0.07941805 

CZD3 

= 

-0.003408011 

CZQ  = 

0.0 

CZD4 

= 

-0.003625448 

CZH  = 

-0.00006055214 

CZD5 

= 

0.00196522 

CZU  = 

0.013001412 

CZD6 

= 

-0.00196522 

CZD1  = 

-0.003086522 

CZD7 

= 

0.00196522 

CZD2  = 

-0.010978221 

CZD8 

= 

-0.00196522 

CM  A  = 

0.004264454 

CMD3 

— 

-0.00362015 

CMQ  = 

-0.1564511 

CMD4 

= 

-0.00451965 

CMH  = 

-0.000027818 

CMD5 

0.00103652 

CMU  = 

-0.0049774 

CMD6 

= 

-0.00115411 

C  M  0 1  = 

0.0066521 

CMD7 

= 

0.00103652 

C  M  D  2  = 

-0. 0125407 

CMD8 

= 

-0.00115411 

CXA  = 

-0.00095604 

CXD3 

— 

0.000578922 

CXQ  = 

0.0 

CXD4 

= 

0.000354177 

CXH  = 

0.000984122 

CXD5 

= 

-0.00201548 

CXU  = 

-0. 193601 

CXD6 

= 

-0.00201548 

CXD1  = 

-0.000741196 

CXD7 

= 

-0.00201548 

CXD2  = 

-0.001530024 

CXD8 

= 

-0.00201548 

ho o end in  B 


State  Space  Models 


This  appendix  contain  the  state  space  models  used  to 
derive  the  plant  transfer  functions.  The  three  separate 
tables  represent  these  state  space  models  for  the  three 
flight  conditions.  The  matrices  are  generated  from 
the  aerodynamic  data  in  Appendix  A  by  the  program 
3T0LCAT. FOR ,  a  listing  of  which  is  in  Appendix  F.  The  state 


vector  x  is: 


x 


B-l) 


with  the  output  vector  y  set  equal  co  x.  The  state  equations 
are  of  the  form: 


x  =  A  >j  +  3  u  (B-2) 

Y  -  C  x  ( 3-3  > 

The  input  vector  u  represents  the  control  surface 

deflections  £  .  The  flight  path  angle  y  is: 

i 

y  =  $  -  a  =  x  x  (B-4> 


t  ' 


0. 0017287 

-4 1 

0.  001630.34 

-0. 

-0. 00147334 

0. 

0.  0 

. 1 037  2. 42547 

387155  1.40620 

373867  -0. 627534 

1.0  0.0 


-31. 


0.  0 


-0. 


0327601 


0.  0 


J 


—2. 43336 
6. 3 l 3132 

—  0  .  0  C-  -  0  /  *+ 
0.  0 


— 0 .  441 3 0. 7 
-  * . 433 1 4 
-0.  066476 
0.  0 


4 .  Z.J.  363 
“0.  W  i.  0  £  7 
-0. 036040 
0.  0 


—  10.  i  ■  3  0  o 
0. 75054 

0  .  0  J.  lX.  1 

0.  0 


'*  ±  *0  .  0  W  0  w 

-0. 362174 

.  4. 

0.  0 


01.  0 


0.0  0.0 

t'.0  0.0 


L. 


"s 


0.  0143226 
e.  00137631 
-0.  0006624c: 
0.  0 


— 20.  D  /  6c' 

la 

0. 997636 

1.0 


31.2520 
2. 45191 
-0. 632834 
0.  0 


—2c.  lc'58 

0.  0 

-0. 0072024 

0.  0 


Cl  m  ZiCiCiClC/ 

-0. 158766 

0093 

— 10. 87c 1 

-10. 37 

. 373&I 

-2. 967 1 7 

— v.  3660 1 

0.  75232 

-0. 27Z 

0. 028296 

-0. 035536 

— 0.  05S232 

•0.  13692 

—  -  1 
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This  appendix  contains  the  individual  control  surface 
deflection  to  response  output  transfer  functions  for  the 
state  space  models  in  Appendix  B.  The  transfer  functions  are 
calculated  using  the  CAD  package  T0TAL[7] .  The  transfer 
functions  are  arranged  according  to  the  output  variable,  with 
the  responses  for  a  single  output  to  all  of  the  surfaces  over 
the  three  flight  conditions  arranged  in  a  single  table.  The 
output  variables  are  the  angle  of  attack  a ,  pitch  angle  9 , 
and  the  forward  velocity  v.  These  transfer  functions  do 
not  contain  tne  servo  poles. 


*.  *.  \  A  V  \ 


TABLE  C— 2 


Pitch  Annie  Response  to  Control  Surface  Input  Transfer  Functions 


Flight  Condition  One: 


0/j  =  C0. 6347 (s+0. 1374) (s+0. 2527) 3/ ^ 


6 /  i  =  [  — 1 . 178  (s+0. 0B939)  (s+0. 4035) 3/ 


0 


1 


0/j  =  [-J. 1472(s+0. 1011) (5+0.4798) 1/  ^ 


0/  x  =  [0. 7946  (s+0. 1556)  (s+0. 3015) 1/ 


tv 


=  [—0.5014 (s+0. 03614) (s+0. 5029) 3/  o 
b  v  1 


0  =  (s+1.546) (s-0.2846) (s+0. 01897+J0. 2508) 
1 


Flight  Condition  Two: 

0/,  =  [0.8191 (s+0. 01344) (s+0. 5747) 3/ 


0. 


6/s  =  [-1.498(5+0.006778) (s+0. 6854) 3/  g 


0/j  =  [-0.3101 (s+0. 01372) (s+0. 771 1 ) 3/  £ 


d/x  =  [0.7506(5+0.03004)  (s+0. 6347.)  3/ 


tv 


0„ 


0/a  =  [-0.3822(5-0.02779) (s+0. 0115) 3/ 


bv 


0 , 


S  =  (s+l. 979) (s-0.3Sri) (s+0. 01 104+10.2081) 


Flight  Condition  Three 


tv 


bv 


[ 1.880 (s+0. 04249) (s+0. 7665) 3/  a 

*3 

[-2. 967(5+0.04023) (5+0.8865)1/  Q 
[-0.8661 (s+0. 04034) (s+0. 9704) 3/ q 

[0. 7524 (s+0. 04998) (s+0. 8427) 3/  a 

PZ. 

[-0.2738(5+0.01886)  <s+1.062>3/  ^ 

* 

(s+2.686) (s-0.4608) < s+0. 003206+ j0. 1090) 


Flight  Condition  One: 


f 

L'- 


«/ 


6  ~ 


a 

tv 

bv 


P 


1 


C-0.2740(s+0.01895+J0.2564) (s-21.41) 3/ 

C-0. 7660(s+0. 01834±j0. 2571) (s+15.49) 3/ 

C-0. 2744 (s+0. 01734* j0. 2572) (5+6. 045) 3/ 

[0. 3890 ( s+0. 01 799± j0. 2564  >  (s+20.80>3/ 

C -0.3890 (s+0. 01900±j0. 2583) (5+12.61)3/ 

(s+1.546) (s— 0. 2846) (s+0. 01897* j0. 2508) 


Fi 1 ght 


Condition  Two: 

a/0  =  C-0. 2107(s+0. 01501+J0. 2171 >  (s-37.26) 3/ 


c 


a 

tv 

bv 


0  = 


C-0. 6648(s+0.01504±j0. 2182)  (s+23.01) 3/ 

C-0.  3604 (s+0. 01434*j0. 2169)  (s+9.  444)  3/ 

C0. 3032(5+0. 01387±j0. 2177)  (s+25.68) 3/ 

C— 0. 3032 (s+0. 01719±J0. 2195) (s+12.80) 3/ 

(s+1.979) ( s-0. 3851 > (s+0. 01104+J0. 2081) 


Flight  Condition  Three: 


C-0.  2830 (s+0. 01037+J0. 1445)  (s-64.66)  1/ 

[—0.9554 (s+0. 01066+ j0.  1445)  <s+32.36>  3/ 

t— 0. 5823 (s+0. 01030+ j0. 1485) (s+16.23) 3/ 

[0. 1989(5+0.007381+ j0. 1466) <s+39.48) 3/ 

t-0. 1989 (s+0. 01938+ j0.  1485)  (s+14.73) 3/ 

(s+2.686) ( s—0. 4608) ( s+0. 003206+ j0. 1090) 


iacTi] 


»j 


TABLE  D-l 

Plant  Matrix  for  FCi  No  Failures 


28.87  (S+0. 01 826+30. 2567) 
(S+19.29HS+31.74) 


-16.20  (S+0. 01 827+30. 2568) 
(S+22.03) (S+35. 87) 


-216.6  (S+0. 3396 jO. 4492) 
(S+32 .02+3  3. 985) 


180.8  (S+0. 3432+jO. 4449) 
(S+32. 33) (S+53. 28) 


TABLE  D-2 

Plant  Matrix  for  FCI  Canard  Failed 


29.29  (S+0. 01823+j0.2568) 
(S+16.64)(S+32.64) 


-17.03  (S+0. 01817+j0.2569) 
(S+l 5 . 24) (S+43. 13) 


-193.9  (S+0. 3i9tt]0. 4590) 
(S+28. 12+J8. 148) 


134.4  (S+0. 3593±30. 4578) 
(S+38. 28ij 12. 30) 


TABLE  I)- 3 

Plant  Matrix  for  FCi  Ailerons  Failed 


28.45  (F-tO. 01828+.) 0,2566) 
(S+20.59)(S+:«.79) 


-15.36  (S+0. 01831110.2567) 
(F+?7.94±j4.67P) 


-194.5  (S+0. 3414110.4439) 
(S  +  7<*.  1 7±j  1 3.  55) 


134.5  (S+0. 3450lj0. 4371) 
(S+41 .24+120.57) 


TABLE  0-4 


Plant  Matrix  for  FC1  Stahllator  Failed 


F' 


»  1 

A, 


27.27  (S+0. 0 1 828+ j 0.2564) 
(S+20. 31) (S+28. 48) 


-13.89  (S+0. 01832+jO. 2567) 
(S+25.04+)7.554) 


-213.9  (S+0. 3192+0.4502) 
(S+14.46) (S+43. 13) 


174.3  (S+0. 31 14+) 0.4423) 
(S+14. 25) (S+60. 93) 


TABLE  D-5 

Plant  Matrix  for  FC1  Stahllator  and  Aileron  Failed 


P’ 


A, 


28.12  (S+0. 01825+J0. 2567) 
(S+15.10) fS+33.54) 


-14.69  (S+0. 01808+)0.2566) 
(S+10.29) (S+39. 16) 


-191.7  (S+0. 32051)0. 4419) 
(S+18.24) (S+35.08) 


129.9  (S+0. 3124+jO. 4250) 
(E+21.01) (S+47.60) 


TABLE  D-6 

Plant  Matrix  for  FC1  Stahllator  and  Canard  Fal 1 ed 


r 


=  l 


28.12  (S+0. 01822+JO. 2563) 
(S+ 16.87) (S+32. 18) 

-191.2  (S+C. 32611)0.4703) 
(S+l 1 . 58) (S  +  39, 59) 


-14.58  fS+0.01P21+J0.2569) 
(S+l 5.11)  (S+35.91) 

128.9  (S+0. 3192+JO. 4893) 
(S+8.665) (S+33.6R) 


D-  a 


TABLE  D-10 

Plant  Matrix  for  FC2  Stabilator  Failed 


■ 

44.75  (S+0. 01466+J0. 2181) 

-24.73  (S+0. 01471+10. 2182) 

(S+20. 34) (S+33. 09) 

(S+19. 60) (S+38.95) 

-374.3  (S+0. 6577+10.3211) 

174.3  (S+0. 6819+10. 2773) 

- 

(S+23. 16±)6.916) 

(S+29. 91+119. 08) 

TABLE  D-Il 

Plant  Matrix  for  FC2  Stabilator  and  Aileron  Fulled 


■ 

43.50  (S+0. 0I462±j0. 2181) 

-22.24  (S+0. 01457±JO. 2183) 

(S+lfi . 22) (S+32. 99) 

(S+l 3. 77) (S+37.44) 

-395.3  (S+0. 6072±j0. 3609) 

301.3  (S+0. 6597510. 3404) 

- 

(S+l 1 . 48) (S+35.54) 

(S+12.70XS+4B.89) 

TABLF.  D-12 

Plant  Matrix  for  FC2  Stabilator  and  Canard  Failed 


43.73  (£*0.01444+10.2180) 
(S  +  19.98MS  +  32.09) 

-380.3  (S+0. 6199+10. 3847) 
(S+10.06) (S+34.41) 


-37. 64  (S  +  0. r’.462+)0. 2183) 
(S+16.27) (P+35.49) 

128.9  (S+0. 6275+10. 3960) 
(S+9. 553) (S+43. 13) 


TABLE  D-13 

Plant  Matrix  for  FC2  Bottom  Vanes  Failed 


* 

28.23  (S+0. 0141 4+J0. 2178) 

-14.74  (S+0. 014491J0. 2179) 

(S+30.91+J3.114) 

(S+33.98) (S+39.86) 

-500.9  (S+0. 6073+JO. 3696) 

332.4  (S+0. 6471+j0. 3059) 

(S+12.75HS+30.66) 

(S+29.54) (S+32.57) 

TABLE  D-14 

Plant  Matrix  for  FC3  No  Failures 


16.08  (S+0. 00993+30. 1463) 

-10.92  (S+0. 01024+jO. 1463) 

(S+34.89+J9. 739) 

(S+40.60+J 12.82) 

-248.2  (S+0. 3492+jl. 1.02) 

206.3  (S+0.36465J1.105) 

(S+27.61+J4.906) 

(S+35.79) (S+46.07) 

TABLE  D-l 5 

Plant  Matrix  for  FC3  Canard  Failed 


TABLE  D-16 

Plant  Matrix  for  FC3  Stabllators  Failed 


* 

14.62  (S+0.00976+J0. 1461) 

-7.997  (S+0. 00999+jO. 1460) 

(S+32. 31+1 10.55) 

(S+37. 35+119. 29) 

-250.6  (S+0. 2675+11. 060) 

2 11.1  (S+0. 3092+11. 066) 

- 

(S+11.65) (S+40.53) 

S+14.82) (S+58.96) 

TABLE  D-17 

Plant  Matrix  for  FC3  Ailerons  Failed 


p 

-1 

15.19  (S+0. 00991+10. 1461) 

-9.139  (S+0. 01023±0. 1461) 

(S+33.85+1T2.59) 

(S+41 .84+120.87) 

-229.8  (S+0. 3482+1 1.087) 

169.4  (S+0. 3888+1 1.88) 

(S+25. 67*19.221) 

(S+38. 55+115. 33) 

TABLE  P-13 

Plant  Matrix  for  FC3  Stabllators,  Ailerons  Failed 


P '  -  2  j 

A 


13.73  <S+0.00971±0.  145$) 
(S+30. 99+113. 27) 

-232.3  (S+0. 2467+) 1.017) 
(S+10.90MS+37.22) 


-6.129  (S+0. 00993+jO. 1455) 
(S+32. 97+1  27 . 80) 

160.6  IS +0.2935+jl,011) 

(S+14.81XS+52.64) 


TABLE  D-19 


Plant  Matrix  for  FC3  Stabllator,  Canard* 


TABLE  D-23 

Matrix  for  FC1  Ailerons  Failed 


■ 

1 

1 

- 

135.4  (S+0.34S91J0.4371) 

15.36  (S+0.018311J0.2568) 

(S+41.24tj 20.52) 

<S+27.94+j4.878) 

1 

1 

TABLE  D-24 

0  Matrix  for  FC1  Stabllator  Failed 


1 

1 

13.89  (5+0.01832*30.2567) 

(S+14 . 25) (S+60.93) 

(S+25 .04IJ7.554) 

1 

1 

.9  (S+0.31921J 0.4502) 
(S+14.46)(S+43.13) 


27.72  (S+0.01828ij0.2567) 
(S+20.31)(S+28.48) 


A-  1674  (S+0. 3570*10. 4401 )(S+1 7. 89) (S+5B. 77) 
(S+1.546) (S-0.2846) 


TABLE  D-25 

0  Matrix  for  FC1  Stabllato  .  Ailerons  Failed 


129.9  (S+0.3124ij0.4250) 
(S+21.01) (S+47.60) 


14.69  (S *0.01 808 tjO. 2566) 
(5*10.29) ( S ♦ 39 . 16) 


'•  .VtV 


Matrix  for  FC2  Stab 11a tor.  Ailerons  Failed 


285  <S+0.5975+jO.3404) 
(S+12.7) (S+46.89) 


22.24  (S+0.01457±J0. 283) 
(S+13 . 77) (S+37 .44) 


395.3  (S+0. 6072tj0. 3609) 
(S+l 1 . 48) (S +35. 54) 


43.50  (S-fO.01462ijO.2181) 
(S+ 18 . 22) (S+32. 99) 


2646  (S+0.5834+<0.3083)(S+15.39ij3. 131) 
(S+l. 979) (S-0.385I) 


TABLE  D-31 


I 


Matrix' for  FC2  Stabllator,  Canard  Failed 


(S+0. 6275ij0. 3960) 
(S+9. 553) (S+43. 13) 


(S+0. 6199tj0. 3847) 
(S+10.06) (S+34.41) 


22.69  (S+0. 01462+j0. 2182) 
(S+18.27)(S+35.49) 


5.73  (S+0. 01494+j0. 2181) 
(S+19.98)(S+32.09) 


A”  1908  (S+0.6496+jP. 4263)  (S+8.661)  (S+23. 253 
(S+l .979^(S-0. 3851)  r~  ‘ 


TABLE  n-3? 


0  Matrix  fr.r  FC2  Bottom  "jnos  Failed 


332.4  (S+0. 6471+jO. 3059) 
(S+29.54)(S+32.57) 


500.9  (S+0. 6073+J0. 3696) 
(S+l 2. 75) (S+30. 66) 


14.74  (S+0. 01449+)0.2179) 
(S+35.98) (S+39.86) 


28.23  (S+0. 01414+10.2178) 
(S+  30. 91  +  j  3. 1 14) 


a  -  615  (S+0.6567+j0.2744)(S+33.53)(S+95.98) 
^  (S+l. 979)(S-0. 3851) 


TABLE  E-l 


)er  Frequency  Bound  -for  AOA  Channel  T 


(RAD/SEC) 

DECIBELS 

DEGREES 

10. 8854 

0. 

. 100000E-01 

10.8854 

- . 286667E— 0 1 

. 200000E-0 1 

10.8855 

-.573336E— 01 

.300000E-01 

10.8855 

- . 860009E— 0 1 

. 400000E-0 1 

10. 8855 

-. 114669 

. 500000E-0 1 

10.8856 

-. 143338 

. 600000E-01 

10. 8857 

172007 

. 700000E— 01 

10.8858 

-.200678 

. 800000E-01 

10. 8859 

-. 229351 

. 900000E-01 

10. 8860 

-.258025 

. 100000 

10.8861 

-.286700 

. 200000 

10.8881 

— . 573603 

. 300000 

10.8915 

-.860911 

. 400000 

10.8962 

-1. 14883 

.500000 

10.9023 

-1.43756 

. 600000 

10.9097 

-1. 72730 

.700000 

10.9184 

-2.01828 

. 800000 

10.9285 

-2.31068 

.900000 

10. 9400 

-2.60473 

1 . 00000 

10.9528 

-2. 90063 

2.00000 

11. 1556 

-6.01052 

3. 00000 

11.4950 

-9.56744 

4 . 00000 

1 1 . 9703 

-13.8706 

5. 00000 

12. 5721 

-19.3200 

6. 00000 

13. 2689 

-26.4513 

7. 00000 

13.9784 

-35.9109 

8. 00000 

14.5284 

-48.2084 

9.00000 

14.6567 

-63.0524 

10.0000 

14. 1497 

-78.7284 

20.0000 

1. 12672 

-133.185 

30.0000 

-6. 13492 

-134. 329 

40. 0000 

-10.6410 

-130. 767 

50. 0000 

-13. 8041 

-126. 703 

60. 0000 

-16. 1883 

-122.963 

70. 0000 

-18.0757 

-119.709  . 

80. 0000 

-19.6252 

-116.926 

90.0000 

-20. 9337 

-114.552 

1 00 . 000 

-22. 0633 

-112.521 

SHHE 


llEHUU 


W  CRAD/SEC) 


'  . 100000E-01 

. 200000E-01 
. 300000E— 01 
. 400000E-01 
. 500000E— 01 
.600000E-01 
. 700000E— 01 
. S00000E— 01 
. 900000E-01 

. 100000 
. 200000 
.300000 
' . 400000  ' 

. 500000 
. 600000 

e..  .700000 

.800000 
. 900000 

1 . 00000 
2. 00000 
3 . 00000 
4.00000 
5. 00000 
6. 00000 
7 . 00000 
8. 00000 
9.00000 

10. 0000 
20. 0000 
30. 0000 
40. 0000 
50. 0000 
60. 0000 
70. 0000 
60. 0000 
90. 0000 
100. 000 


DECIBELS 

DEGREES 

10.3814 

0. 

10. 8813 

-.212813 

10. 8812 

-.425623 

10. 8811 

-.638430 

10.8809 

-.851232 

10.8806 

-1.06403 

10.8803 

-1.27681 

10.8799 

-1.48958 

10.8794 

-1.70235 

10.878? 

-1.91509 

10. 8783 

-2. 12782 

10. 8692 

-  -4.25380 

10.8540 

—6. 37610 

10.8329 

-8.49293 

10.8058 

—10. 6025 

10.7729 

-12.7032 

10-7342 

-14.7932 

10.6899 

-16.8710 

10.6402 

-18.9352 

10.5851 

-20. 9843 

9.78131 

-40.4057 

8.65155 

-57.4954’ 

7 . 35993 

-72.2553 

6.01007 

-85. 0303 

4.65594 

-96.2078 

3.32246 

-106. 1 10 

2.01991 

-114.981 

.  751.756 

-123.000 

-.481476 

-130.301 

-11. 1086 

-178.698 

-19.2908 

-203. 698 

-25. 7721 

-218.434 

-31.0726 

-227.982 

j j ■ 5o04 

-234.617 

—39 . 3655 

-239.474 

-42. 7248 . 

-243. 175 

-45.7106 

-246.085 

-48.3959 

-248.430 

W  < RAD /SEC) 


DECIBELS 


DEGREES 


0. 

13. 9794 

0. 

.  100000E-01 

13.9790 

-.572939 

.  200000E-01 

13.9777 

-1. 14576 

. 300000E-0 1 

13.9755 

-1.71836 

. 400000E-01 

13.9725 

-2. 29061 

,  500000E-01 

13.9686 

-2.86241 

.  600000E-01 

13. 9638 

-3. 43363 

.  700000E-01 

13.9582 

-4.00417 

.  800000E-01 

13.9517 

-4.57392 

. 900000E— 01 

13.9444 

-5. 14276 

.  100000E+00 

13.9362 

-5.71059 

. 100000 

13.9362 

-5.71059 

. 200000 

13. 8091 

-11. 3099 

. 300000 

13.6051 

-16.6992 

. 400000 

13.3348 

-21 v 8014 

-500000 

13.0103 

-26.5651 

. 600000 

12.6440 

-30.9638 

. 700000 

12.2475 

-34. 9920 

. 800000 

11.8310 

-38.6593 

. 900000 

1 1 . 4026 

-41.9872 

1.00000 

10.9691 

-45.0000 

1.00000 

10.9691 

—45 . 0000 

2. 00000 

6.98970 

-63. 4349 

3. 00000 

3. 97940 

-71 . 5651 

4.00000 

1 . 67491 

-75.9638 

5.00000 

-. 170333 

-78. 6901 

6. 00000 

-1  .'70262 

-80. 5377 

7.00000 

—3^01030 

-81.8699 

8.00000 

-4. 14973 

-82. 8750 

9.00000 

-5. 15874 

-83.6598 

10. 0000 

— 6 . 06 -'8 1 

-84. 2894 

10. 0000 

-6. 06381 

-84.2894 

20. 0000 

-12. 0520 

-87. 1376 

30.0000 

-15.5678 

-88.0908 

40. 0000 

-18.0645 

-88.5679 

50. 0000 

-20.0017 

-88.8542 

60. 0000 

-21 . 5843 

-89. 0452 

70.0000 

-22 . 9234 

-89. 1815 

30 . 0000 

-24. 0831 

-89.2838 

90 . 0000 

-25 . i 060 

-89 . 3634 

1 00 . 000 

- 2 :j  .  021 0 

-89.4271 

1  < RAD/SEC > 

DECIBELS 

DEGREES 

13.9794 

0. 

. 100000E— 01 

13. 9702 

-1.02360 

. 200000E— 0 1 

13. 9746 

-2.04666 

. 300000E— 01 

13. 9685 

-3.06867 

. 400000E— 01 

13. 9601 

-4.08909 

. 500000E— 0 1 

13. 9492 

-5. 10741 

. 600000E— 0 1 

13. 9360 

-6. 12312 

. 700000E-01 

13. 9205 

-7. 13570 

. 800000E-01 

13. 9026 

-8. 14467 

. 900000E-0 1 

13. 8824 

-9. 14955 

•  1 00000E+00 

13.8600 

-10. 1499 

. 100000 

13.8600 

-10. 1499 

. 200000 

13.5200 

-19. 8099 

. 300000 

13.0062 

-28.6272 

. 400000 

12.3752 

-36. 4390 

. 500000 

11.6778 

-43.2409 

. 600000 

10.9529 

-49.  121 1 

. 700000 

10.2262 

-54.2050 

. 800000 

9.51361 

-58. 6207 

. 900000 

8.82412 

-62.4839 

1 . 00000 

8. 16209 

-65.8926 

1 . 00000 

8. 16209 

-65.8926 

2.00000 

2.97028 

-86.9013 

3.00000 

-.560205 

-98.8229 

4. 00000 

-3.26769 

-107.845 

5.  00000 

-5.51138 

-115.433 

6.00000 

-7. 46130 

-122.096 

7.00000 

-9.20749 

-128. 063 

8.  00000 

-10. 8020 

-133. 461 

9.  00000 

-12.2770 

-138. 377 

10.  0000 

-13.6538 

-142.376 

10. 0000 

-13. 6538 

-142. 876 

20. 0000 

-24. 1164 

-173.518 

30. 0000 

-31.3371 

-191. 3B3 

40. 0000 

-36. 9527 

-203.764 

50. 0000 

-41 . 5976 

-2 1 3 . 003 

60.0000 

-45.5769 

-220. 159 

70. 0000 

-49.0625 

-225.841 

80. 0000 

-52. 1633 

-230.440 

90.0000 

-54.9544 

-234.223 

100. 000 

-57. 4906 

-237. 381 

Appendix  F 


STOLCAT,  written  by  the  author  and  Captain  Greg 
Mandt(AFWAL/FIGX)  takes  body  axis  aerodynamic  data  and 
transforms  it  into  lateral  and  longitudinal  state  space 
aircraft  equations  of  motion  in  the  stability  axis.  The  form 
of  these  equations  is  in  Appendix  B.  This  thesis  uses  only 
tne  longitudinal  equations,  but  provides  the  lateral 
equations  for  other  Air  Force  Institute  of  Technology 
students  working  with  the  same  aircraft  in  their  thesis.  The 
stability  axis  is  chosen  as  the  output  axis  to  make  the  other 
student's  analysis  simpler.  The  choice  has  little  effect  on 
this  thesis  since  design  is  limited  to  the  longitudinal  mode. 
The  program  is  in  FORTRAN  77,  and  set  up  to  run  interactive, 


rather  than  batch. 


PROGRAM  STOLCAT 

C************************************************************************* 

c 

C  DECLARE  VARIABLE  TYPES 
C 

C******************************************************************* ****** 

REAL  ALPHA,  Q,  S,  C,  B,  U,  DTHETA,  W,  B IXX,  B I YY,  BI ZZ, 

1BXXZ  DALPHA  DPR  VT 

2CZA,  CZQ,  CZu!  CZD1,  CZD2,  CZD3,  CZD4,  CZD5,  CDZ6,  CDZ7,  CZD8, 

3CXA,  CXQ,  CXU,  CXD1, CXD2,  CXD3,  CXD4,  CXD5,  CXD6,  CXD7,  CXDB, 

4CMA, CMQ, CMU, CMD1 , CMD2, CMD3, CMD4, CMOS,  CMD6, CMD7, CMD8, 

5Z1,  ZA,  ZH,  ZQ,  ZU,  ZD1,  ZDS,  ZD3,  ZD4,  ZDS,  ZD6,  ZD7,  ZD8, 

6XA,  XH,  XO,  XU,  XD1,  XDS,  XD3,  XD4,  XDS,  XD6,  XD7,  XD8, 

7M1 ,  MA,  MH,  MQ,  MU,  MD1 ,  MDS,  MD3,  MD4,  MD5,  MD6,  MD7,  MD8 
REAL  CNB,  CYB,  CLB,  L,  N 
DIMENSION  AMAT (4, 4) , BMAT (4, 8) 

DIMENSION  DI RMAT (S, 5) , DI RBMAT  <5, 9) 

CHARACTER*3  KEY,  KEY1,  DATA1,  DATA?,  DATA3,  RUN 
CHARACTER* 1  STAB  1 , ST AB2 

C************************************************************************* 

C 

C  INITIAL  DATA  VALUES  FOR  PROGRAM  CHECK 
C 

C************************************************************************* 

DATA  □  /48.  1/,S  /608.  /,  C  /IS.  94/,  B  /42. 7/,  U  /20I./ 

DATA  DTHETA  /I 1.8030/,  DALPHA  /11.8030/,U  / 33576. 14/ 

DATA  BIXX  /23644.  /,  BIW  /181847.  /,  BIZZ  /199674. /,  BIXZ  /-3086.  / 

DATA  CZA  /-7. 84976E-2/,  CXA  /l. 5095276E-3/,  CMA  /9. 5741 18E-3/ 

DATA  CZQ  /0. / ,  CXQ  /0. /,  CMQ  /-. 16951603/ 

DATA  CZU  /-l. 06551597/,  CXU  /-6. 1932E-3/,  CMU  /6. 3942B9E-2/ 

DATA  CZH  /-1.676463E-4/,  CXH  /6. 662777E-4/,  CMH  /1.76622E-4/ 

DATA  CZD1  /-2.63634E-3/,  CXD1  /-l. 552420E-3/,  CMDI  /5. 57696E-3/ 

DATA  CZD2  /-8. 3151  IE-3/,  CXD2  /-2. 749671E-4/,  CMD2  /-l.  02066E-2/ 

DATA  CZD3  /-5. 59102E-3/,  CXD3  /l. 157373E-3/,  CMD3  /B.52107E-4/ 

DATA  CZD4  /-4. 50843E-3/,  CXD4  /9. 421 1093E-4/,  CMD4  /-2.  1 1 1 1BE-3/ 

DATA  CZD5  / 1. 896349E-3/,  CXD5  /-3. 120989E-3/,  CMD5  /2.  55459E-3/ 

DATA  CZD6  /-7. 422954E-4/,  CXD6  /-3. 595656E-3/,  CMD6  /-l.  30123E-3/ 
DATA  CZD7  /l. 896349E-3/,  CXD7  /-3. 120989E-3/,  CMD7  /2. 55459E-3/ 

DATA  CZD8  /-7. 422954E-4/,  CXD8  /-3. 595658E-3/,  CMD8  /-l.  30123E-3/ 
DATA  CLB  /-2. 973933E-3/,  CNB  /-5. 5065055E-4/,  CYB  /-l.  637941E-2/ 

DATA  CLP  /— 5. 740524E-3/,  CNP  /-2. 3099719E-3/,  CYP  /  0.000000000/ 

DATA  CLR  /  3. 902348E-3/,  CNR  /-9. 6998151E-3/,  CYR  /  0.008800000/ 

DATA  CLD1 / 1 . 0017E-4/,  CND1 /-l. 3256E-3/,  CYD1/3.  0606E-3/ 

DATA  CLD2/-1.  14999E-4/,  CND2/5. 1323E-4/,  CYD2/ 1 . 3139E-3/ 

DATA  CLD3/8. 5184E-4/,  CND3/4. 4837E-4/,  CYD3/-1. 0622E-3/ 

DATA  CLD4/7. 5284E-4/,  CND4/7. 6138E-5/,  CYD4/-1. 5235E-4/ 

DATA  CLD5/6. 9959E-4/,  CND5/0. 00/,  CYD5/®. 00/ 

DATA  CLD6/9. 6816E-5/,  CND6/ 1 . 5934E-4/,  CYD6/0.  0/ 

DATA  CLD7/-3. 7897E-5/,  CND7/ l . 8357E-4  ,  CYD7/0.0/ 

DATA  CLD8/-9. 6816E-5/,  CND8/-1. 5934E-4/,  CYD8/0.0/ 

DATA  CLD9/3.  7897E-5/,  CND9/-1. 8357E-4/,  CYD9/0. 8/ 

DPR  =  57.2957795 

C************************************************************************ 

C 

C  SCREEN  INFO  PACKAGE 

Z  ************************************************************************ 

WRITE!*, 5) 

5  FORMAT (IX, '****************************************************’ > 


non 


WRITEt*, 10) 

10  FORMAT < IX, ’ ***  STABILITY  DERIVATIVE  TRANSFORMATION  PROGRAM  ***’ ) 
WRITEt*, 20) 

£0  FORMAT (IX,1 ****************************************************’  ) 
WRITEt*, 100) 

100  FORMAT (IX, ’ENTER  BODY  AXIS  ( NON-D I MENS I ONAL I  ZED)  COEFFICIENTS  ’) 
WRITEt*, 101) 

101  FORMAT (IX, ’FOR  TRANSFORMATION  TO  DIMENSIONAL  I  ZED  BODY  AXIS’) 
WRITEt*, 102) 

102  FORMAT (IX, ’AND  TO  GENERATE  STATE  AND  INPUT  MATRICES.’) 

WRITEt*, 41) 

41  FORMAT (IX, ’NOTE:  ALL  COEFFICIENTS  ARE  REQUESTED  WHEN  COMPUTING’) 

103  CONTINUE 
WRITEt*, 30) 

30  FORMAT (IX,’ ****************************************************’ > 
WRITEt*, 106) 

106  FORMATtlX, ' TO  TRANSFORM  ONLY  LONGITUDINAL  DATA  -  TYPE  LONG’) 
WRITE(*, 107) 

107  FORMAT (IX, ’TO  TRANSFORM  ONLY  LATERAL-DIRECTIONAL  DATA  -  TYPE  LAT’ ) 
WRITEt*, 108) 

108  FORMATtlX, ’TO  TRANSFORM  BOTH  LONG  AND  LAT-DIR  DATA  -  TYPE  BOTH’) 
WRITEt*, 111) 

111  FORMAT (IX, ’KEYWORD  «  •) 

READ (*,109)  KEY 

109  FORMAT (A3) 

IF (KEY  . EQ.  ’ LAT’ )  GO  TO  104 
IF (KEY  .EQ.  ’LON’)  GO  TO  104 
IF (KEY  .EQ.  ’BOT’)  GO  TO  104 
IF (KEY  .EQ.  ’GAM’)  GO  TO  396 
GO  TO  103 

C********************************************************************** 
INPUT  DATA 

C********** ************************************************************ 

104  CONTINUE 
WRITEt*, 500) 

500  FORMAT (IX,’ ****************************************************’ ) 
WRITEt*, 510) 

510  FORMAT ( 1 X, ’ Q  (DYNAMIC  PRESSURE  -  LBS/FT**2)  *  ’) 

READ  ( *,  * )  0 
WRITE(*, 520) 

520  FORMATtlX, ’S  (WING  REFERENCE  AREA  -  FT**2>  -  ’) 

READt*,  *)  S 
WRITE ( *, 530) 

530  FORMATtlX, ’C  (WING  MEAN  AERODYNAMIC  CORD  -  FT)  =  ’> 

READ (*, *>  C 
WRITEt*, 540) 

540  FORMAT ( IX, ’B  (WING  SPAN  -  FT)  *  ’  > 

READt*,  *)  B 
WRITEt*, 350) 

550  FORMATtlX, ’ VT  (TRIM  VELOCITY  -  FT/SEC)  *  ’> 

READ  (*,  *>  U 
VT=U 

WR I TE (*, 560) 

560  FORMATtlX,  '  THETA  (PITCH  ANGLE  -  DEGS)  =■  ’) 

READt*,*)  DTHETA 
WRITEt*, 570) 

570  FORMATtlX, ’W  (WEIGHT  -  LBS)  »  ’) 

READ!*,*)  W 


PROGRAM  STOLCAT 

C*********************************** ************************************** 

c 

C  DECLARE  VARIABLE  TYPES 

C 

C************************************************************************* 
REAL  ALPHA,  Q,  S,  C,  B,  U,  DTHETA,  W,  BIXX,  B1YV,  BIZ Z, 

1BIXZ  DALPHA  DPR  VT 

2CZA,  CZO,  CZuJ  CZD1,  CZD2,  CZD3,  CZD4,  CZD5,  CDZ6,  CDZ7,  CZDQ, 

3CXA,  CXO,  CXU,  CXD1,  CXD2,  CXD3,  CXD4,  CXD5,  CXD6,  CXD7,  CXD8, 

4CMA,  CMQ,  CMU,  CMD1,  CMD2,  CMD3,  CMD4,  CMOS,  CMD6,  CMD7,  CMDS, 

5Z1,  ZA,  ZH,  ZO,  ZU,  ZD1,  ZDS,  ZD3,  ZD4,  ZDS,  ZD6,  ZD7,  ZD8, 

6XA,  XH,  XQ,  XU,  XD1,  XD2,  XD3,  XD4,  XDS,  XD6,  XD7,  XD8, 

7M1 ,  MA,  MH,  MO,  MU,  MD1 ,  MD2,  MD3,  MD4,  MD5,  MD6,  MD7,  MD8 
REAL  CNB,  CYB,  CLB,  L,  N 
DIMENSION  AMAT (4, 4) , BMAT <4, 8) 

DIMENSION  DIRMAT <3,  5) ,  DIRBMAT (5,  9 ) 

CHARACTER*3  KEY,  KEY1,  DATA1,  DATA2,  DAT  A3,  RUN 
CHARACTER* 1  STAB 1 , STAB2 

C***********************************************************************»* 

c 

c  INITIAL  DATA  VALUES  FOR  PROGRAM  CHECK 

C 

C************************************************************************* 

DATA  Q  /48.  I /,  S  /60S.  /,  C  /IS.  94/,  B  /42.7/,  U  /201./ 

DATA  DTHETA  /I 1,8030/,  DALPHA  /11.B030/.W  / 33576. 14/ 

DATA  BIXX  /23644.  /,  BI YY  /181847. /,  BIZZ  /199674. /,  BIXZ  /-30B6.  / 

DATA  CZA  /-7.  04976E-2/,  CXA  /t. 5095276E-3/,  CMA  /9.  5741 1BE-3/ 

DATA  CZQ  /0.  /,  CXO  /0.  /,  CMQ  /-.  16951603/ 

DATA  CZU  /-l. 06551597/,  CXU  /~6. 1932E-3/,  CMU  /6. 3942B9E-2/ 

DATA  CZH  /-1.676463E-4/,  CXH  /6. 662777E-4/,  CMH  /1.76622E-4/ 

DATA  CZD1  /-2.63634E-3/,  CXD1  /-l. 552420E-3/,  CMD1  /5. 57696E-3/ 

DATA  CZD2  /-8.  3151  IE-3/,  CXD2  /-2. 749671E-4/,  CMD2  /-l. 02066E-2/ 

DATA  CZD3  /-5. 59102E-3/,  CXD3  /l. 157373E-3/,  CMD3  /B.52107E-4/ 

DATA  CZD4  /-4. 50843E-3/,  CXD4  /9. 421 1093E-4/,  CMD4  /-2. 1 1 1 1BE-3/ 

DATA  CZD5  /l.  896349E-3/,  CXD5  /-3. 1209B9E-3/,  CMD5  /2. 55459E-3/ 

DATA  CZD6  /— 7. 422954E— 4/,  CXD6  /-3. 595656E-3/,  CMDS  /-l. 30123E-3/ 
DATA  CZD7  /1.B96349E-3/,  CXD7  /-3. 120989E-3/,  CMD7  /2. 55459E-3/ 

DATA  CZDS  /-7. 422954E-4/,  CXD8  /-3.  S9565BE-3/,  CMDS  /-l.  30123E-3/ 
DATA  CLB  /-2. 973933E-3/,  CNB  /-5. 5065055E-4/,  CYB  /-l. 637941E-2/ 

DATA  CLP  /-5. 740524E-3/,  CNP  /-2. 3099719E-3/,  CYP  /  0.000000000/ 

DATA  CLR  /  3. 902348E-3/,  CNR  /-9. 6398151E-3/,  CYR  /  0.000000000/ 

DATA  CLD1/1. 0017E-4/,  CND1/-1. 3256E-3/,  CYDl/3. 0606E-3/ 

DATA  CLD2/-1. 14999E-4/,  CND2/5. 1323E-4/,  CYD2/1. 3139E-3/ 

DATA  CLD3/8. 5104E-4/,  CND3/4. 4B37E-4/,  CYD3/-1. 0622E-3/ 

DATA  CLD4/7. 5284E-4/,  CND4/7. 6 1 38E-5/,  CYD4/-1. 5235E-4/ 

DATA  CLD5/6. 9959E-4/,  CND5/0.00/,  CYD5/0. 00/ 

DATA  CLD6/9. 6816E-5/,  CND6/1. 5934E-4/,  CYD6/0. 0/ 

DATA  CLD7/-3. 7B97E-5/,  CND7/1. 8357E-4/,  CYD7/0.0/ 

DATA  CLD8/-9. 6816E-5/,  CND8/-1. 5934E-4/,  CYD8/C. 0/ 

DATA  CLD9/3. 7897E-5/,  CND9/-1. 8357E-4/,  CYD9/0.0/ 

DPR  =  57.2957795 

C************************************************************************ 

C 

C  SCREEN  INFO  PACKAGE 

C 

C  it-**#****##**********#*************************************************** 
WRITE!*, 5) 

5  FORMAT (IX,' ****************************************************’  > 


n  n  o  n 


1 


WRITE**, 10) 

10  FORMAT* IX,’***  STABILITY  DERIVATIVE  TRANSFORMATION  PROGRAM  ***’  > 

i  WRITE  (*,20) 

20  FORMAT < IX, ' ****************************************************’ ) 

;  WRITE**, 100) 

I  100  FORMAT* IX, ’ENTER  BODY  AXIS  (NON-DIMENSIONAL I  ZED)  COEFFICIENTS  ’) 

WRITE**, 101) 

101  FORMAT* IX, ’FOR  TRANSFORMATION  TO  DIMENSIONAL I  ZED  BODY  AXIS’) 

!  WRITE**,  102) 

I  102  FORMAT (IX, ’AND  TO  GENERATE  STATE  AND  INPUT  MATRICES.’) 

WRITE**, 41) 

41  FORMAT* IX, ’NOTEj  ALL  COEFFICIENTS  ARE  REQUESTED  WHEN  COMPUTING’) 

103  CONTINUE 

WRITE <*, 30) 

30  FORMAT  *  1 X, ’ ****************************************************’ ) 

WRITE**, 106) 

106  FORMAT < 1 X, ’ TO  TRANSFORM  ONLY  LONGITUDINAL  DATA  -  TYPE  LONG’) 
WRITE**, 107) 

107  FORMATUX, ’TO  TRANSFORM  ONLY  LATERAL-DIRECTIONAL  DATA  -  TYPE  LAT’ ) 
WRITE**, 10S) 

108  FORMAT* IX, ’TO  TRANSFORM  BOTH  LONG  AND  LAT-D1R  DATA  -  TYPE  BOTH’) 
WRITE**, 111) 

111  FORMATUX, 'KEYWORD  -  ’) 

READ <*, 109)  KEY 

109  FORMAT (A3) 

IF <KEY  .EQ.  ’LAT’)  GO  TO  104 
IF (KEY  .EQ.  ’LON’)  GO  TO  104 
IF (KEY  .EQ.  ’BOT’)  GO  TO  104 
IF <KEY  .EQ.  ’GAM’)  GO  TO  396 
GO  TO  103 

C********************************************************************** 
INPUT  DATA 


104  CONTINUE 

WRITE**, 500) 

500  FORMAT  *  1 X, ’ ****************************************************’  ) 
WRITE**, 510) 

510  FORMATUX, ’Q  (DYNAMIC  PRESSURE  -  LBS/FT**2>  -  ’) 

READ**,*)  Q 
WRITE (*,320) 

520  FORMATUX, ’S  <WING  REFERENCE  AREA  -  FT**2>  -  ’) 

READ**,*)  S 
WRITE**, 530) 

530  FORMATUX, ’C  <WING  MEAN  AERODYNAMIC  CORD  -  FT)  <*  ’  > 

READ**,*)  C 
WRITE**,  340) 

540  FORMATUX, ’B  <WING  SPAN  -  FT)  -  ’> 

READ**,*)  B 
WRITE**, 350) 

550  FORMATUX, ’VT  (TRIM  VELOCITY  -  FT/SEC)  -  ’) 

READ  <*,*)  U 
VT=U 

WRITE (*,560) 

560  FORMAT (IX, ’THETA  (PITCH  ANGLE  -  DEGS)  =  ’) 

READ**,*)  DTHETA 
WRITE**, 570) 

570  FORMATUX, ’W  (WEIGHT  -  LBS)  -  ’> 

READ**,*)  W 


n  n 


WRITE  (#,575) 

575  FORMAT ( iX, *  INERTIAS  MUST  BE  INPUT  IN  BODY  AXIS.’) 

WRITE (*,  580) 

580  FORMAT ( 1 X,  ’  IXX  <SLUG-FT**2)  =  •> 

READ<*, *)  BIXX 
WRITE ( #,  585) 

585  FORMAT (IX, ’ I YY  (SLUG-FT**2>  »  ’> 

READ(*, *)  BIYY 
WRITE (*, 590) 

530  FORMAT (IX,  ' IZZ  (SLUQ-FT**2>  -  ’> 

READ ( *, *>  BIZZ 
WRITE(*, 595) 

595  FORMAT ( 1 X, ’ I XZ  (SLUG-FT**2>  -  ’> 

READ (*, *)  BIXZ 

596  CONTINUE 
WRITE!#, 597) 

597  FORMAT (IX,*  *##*##****##***#*#*##*###**#*##*##*#**#**#**#**###**’ ) 
WRITE!#, 610) 

610  FORMAT (16X, ’AIRCRAFT  PARAMETERS’) 

WRITE!*, 615)  Q 

615  FORMAT ( 1 X, ' Q  (DYNAMIC  PRESSURE  -  LBS/FT**2>  -  ’,G13.6> 

WRITE (*,620)  S 

620  FORMAT (1X,’S  (WING  REFERENCE  AREA  -  FT**2)  ■  ’,G13.6> 

WRITE (*,625)  C 

625  FORMAT ( I X,  ’ C  (WING  MEAN  AERODYNAMIC  CORD  -  FT)  -  ’,613.6) 

WRITE (*,630)  B 

630  FORMAT ( 1 X, ’ B  (WING  SPAN  -  FT)  -  ’,613.6) 

WRITE  (*,635)  U 

635  FORMAT (IX,’ VT  (TRIM  VELOCITY  -  FT/SEC)  «  ’,G13.6> 

WRITE (*, 640)  DTHETA 

640  FORMAT ( 1 X, ’ THETA  -  ’.G13.6) 

WRITE ( *, 645)  W 

645  FORMAT (IX, ’W  (WEIGHT  -  LBS)  «  ’.G13.6) 

WRITE (*, 650)  BIXX 

650  FORMAT (IX, 1 IXX  (SLUG-FT**2)  -’,613.6) 

WRITE(*, 655)  BIYY 

655  FORMAT ( 1 X, ’ I YY  (SLUG-FT**2)  -’,613.6) 

WRITE (*, 660)  BIZZ 

660  FORMAT ( IX, ’ IZZ  (SLUG-FT**2)  -  ’,G13.6> 

WRITE ( *, 665)  BIXZ 

665  FORMAT ( 1 X, ' IXZ  (SLUG-FT**2)  =’,G13.6> 

WRITE!*,  670) 

670  FORMAT (IX, ’ **********#****#*#************#*****#**##***#*##****’ ) 

600  CONTINUE 

WRITE (*, 675) 

675  FORMAT ( 1 X, ’ IS  THE  ENTERED  DATA  CORRECT  ?  (YES/NO)  ’) 

READ (#,680)  DATA3 

680  FORMAT (A3) 

WRITE!*, 685) 

685  FORMAT (IX, ’ #*#***»*******»*##********##*****#*##**#**»*********,  ) 
IF (DATA3  .EQ.  ’ NO  ’ )  GO  TO  104 

IF (DATA3  .EQ.  ’YES’)  GO  TO  686 
GO  TO  600 

686  CONTINUE 
WRTTE (*, 105) 

105*  FORMAT  ( 1 X,  '  ALPHA  (DEG)  =  ’) 

READ!*,*)  DALPHA 

CHANGE  FROM  DEGREES  TO  RADIANS 


’.--V-W  "*• 


V  VHi  *  ,  w  m  v*jT* 


THETA  -  DTHETA/DPR 
ALPHA  -  D ALPHA/DPR 
IF  (KEY  .EQ.  *  LAT’  )  GO  TO  446 
IF (KEY  .EQ.  ’GAM') GO  TO  97 
C 

C  INPUT  LONGITUDINAL  VARIABLES 
C 

WRITE!*,  110) 

110  FORMAT  (IX, ’CZA  -  ’ ) 

READ!*, *)  CZA 
WRITE!*,  1£0> 

120  FORMAT ( 1 X,  ’ CXA  *  ’) 
READ!*,*)  CXA 
WRITE!*,  130) 

130  FORMAT ( IX,  ’  CMA  «  ’) 
READ!*,*)  CMA 
WRITE!*, 140) 

140  FORMAT ( 1 X,  •  CZQ  -  '  ) 
READ!*,*)  CZQ 
WRITE!*, 1S0> 

150  FORMAT! IX,  ’  CXQ  -  » ) 

READ!*,  *)  CXQ 
WRITE (*, 160) 

160  FORMAT ( 1 X, ’CMQ  -  ’) 
READ!*,*)  CMQ 
WRITE!*, 170) 

170  FORMAT ( 1 X, ’ CZU  -  ’) 
READ!*,*)  CZU 
WRITE!*, 160) 

180  FORMAT ( 1 X,  ’  CXU  -  ’  ) 
READ!*,*)  CXU 
WRITE!*, 190) 

190  FORMAT ( 1 X, ’ CMU  -  ’  ) 
READ!*,*)  CMU 
WRITE!*, 191) 

191  FORMAT ! 1 X,  ’ CZH  =  ’ ) 
READ!*,*)  CZH 
WRITE!*, 192) 

192  FORMAT ( 1 X,  ’CXH  «  ’  ) 
READ!*,*)  CXH 
WRITE!*, 193) 

193  FORMAT ( 1 X, ’ CMH  -  ') 
READ!*,*)  CMH 
WRITE!*, 200) 

£00  FORMAT ( 1 X, 'CZD1  *  ’) 
READ!*,*)  CZD1 
WRITE!*,  202) 

£02  FORMAT (1X,'CXD1  -  ’) 
READ!*,*)  CXD1 
WRITE!*,  204) 

204  FORMAT ( 1 X, 'CMDl  »  ’) 
READ!*,*)  CMDl 
WRITE!*,  206) 

£06  FORMAT ( IX, ’CZD2  =  ’) 
READ!*,*)  CZD2 
WRITE!*, 206) 

£06  FORMAT! IX,  ’  CXD2  *  ’) 
READ!*,*)  CXD2 
WRITE!*,  210) 

£10  FORMAT ( 1 X, ’ CMD2  *  ’> 
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*  ) 


»  ) 


) 


> 


»  ) 


1  ) 


-  »  ) 


1  ) 


READ  !*, *)  CMD2 
WRITE (*, 212) 

FORMAT ( IX, ' CZD3 
READ (*, * )  CZD3 
WRITE  !*,  214) 

FORMAT ! 1 X, ’CXD3 
READ (#, *)  CXD3 
WRITE (*,  2 16> 

FORMAT ! 1 X, ’CMD3 
READ!*,*)  CMD3 
WRITE!*,  218) 

FORMAT !  1 X,  ’ CZD4 
READ!*,*)  CZD4 
WRITE!*, 45) 

FORMAT ( IX, *CXD4 
READ!*,*)  CXD4 
WRITE!*,  50) 

FORMAT <1 X, »CMD4 
READ!*,*)  CMD4 
WRITE!*, 55) 

FORMAT ( 1 X,  *CZD5 
READ!*,*)  CZD5 
WRITE!*, 60) 

FORMAT ! 1 X, •  CXDS 
READ!*,*)  CXDS 
WRITE!*, 65) 

FORMAT ! 1 X, « CMD5 
READ!*,  •>  CMD5 
WRITE!*,  70) 

FORMAT < 1 X, >CZD6 
READ!*,*)  CZD6 
WRITE!*,  75) 

FORMAT ! IX, 1 CXD6 
READ!*,*)  CXD6 
WRITE!*,  80) 

FORMAT < 1 X, ’CMD6  - 
READ!*, *)  CMD6 
WRITE!*, 85) 

FORMAT! IX, > CZD7’  > 

READ!*, *>  CZD7 
WRITE!*, 88) 

FORMAT! IX, * CXD7’ ) 

READ!*,*)  CXD7 
WRITE!*, 90) 

FORMAT ! 1 X, •  CMD7 
READ!*,*)  CMD7 
WRITE!*,  92) 

FORMAT !1X,’CZD8 
READ!*,  *>  CZD8 
WRITE!*,  94) 

FORMAT! IX, *0X08 
READ!*,  *>  CXDS 
WRITE!*,  96) 

FORMAT ! 1 X,  * CMD8 
READ!*,*)  CMD8 
CONTINUE 
WRITE!*, 225) 

FORMAT ! 1 X, * ****************************************************'  > 

WRITE!*, 230)  DALPHA 
FORMAT (15X, ’ALPHA  =’,Gi3. 6) 


,  *  ,v 

• n  ' 


«  j 

Si 


•  > 


*  > 


1  > 


V  ’ 


□ 


* ) 


’ ) 


) 


.  •  /V 


) 


n  n  n 


WRITE!*,  345) 

345  FORMAT <6X, •LONGITUDINAL  NON-DIM  BODY  AXIS  COEFFICIENTS ( 1 /DEG) ' ) 

CAL  =  COS (ALPHA) 

SAL  «  SIN(ALPHA) 

COSSQ  »  CAL**2 
SINSQ  =  SAL**2 
COSSIN  =  CAL*SAL 
CTH  =  COS (THETA) 

STH  =  SIN (THETA) 

C 

WRITE (*,  360)  CZA,  CMA,  CXA 

360  FORMAT  <3X,  *CZA  »  ’ ,  G1 3.  6,  8X,  '  CMA  =  * ,  G13.  6,  5X, ’  CXA  =  ’,G13.6> 
WRITE (*,390)  CZQ, CMQ, CXQ 

390  FORMAT <3X, 'CZQ  =  ’ ,  G1 3.  6,  SX,  '  CMQ  *  ' ,  G13.  6,  5X,  •  CXO  »  ’.G13.6) 
WRITE (*,400)  CZH,CMH, CXH 

400  FORMAT  <3X,  ’  CZH  =  ’  ,  G13.  6,  SX,  '  CMH  =  > ,  G13.  6,  5X,  •  CXH  «  ’,613.6) 
WRITE(*,  410)  CZU,  CMU,  CXU 

410  FORMAT  (3X, 'CZU  ■  • ,  G1 3.  6,  SX,  ’  CMU  »  » ,  G13.  6,  5X, '  CXU  «  ',613.6) 
WRITE(*, 370)  CZD1, CMD1 , CXD1 

370  FORMAT  <2X,  'CZD1  -  ’  ,  G13.  6,  7X,  ’  CMD1  =  '  ,  G1 3.  6,  4X,  •  CXD1  »  ’,G13.6) 
WRITE (*, 380)  CZD2, CMD2, CXD2 

380  FORMAT  <2X,  ’  CZD2  -  ' ,  G13.  6,  7X,  '  CMD2  -  ',613.  6,  4X,  ’  CXD2  -  *,G13.6) 
WRITE (*,381)  CZD3, CMD3, CXD3 

381  FORMAT <2X, ’CZD3  »  '  ,  G13.  6,  7X,  ’  CMD3  “  ' ,  G13.  6,  4X, '  CXD3  -  ',G13.6> 
WRITE (*,382)  CZD4.CMD4, CXD4 

382  FORMAT  (2X,  •  CZD4  -  *  ,  G1 3.  6,  7X,  •  CMD4  -  ' ,  G13.  6,  4X,  •  CXD4  -  ',G13.6> 
WRITE (*,383)  CZD5, CMD5, CXD5 

383  FORMAT  <2X,  '  CZD5  -  ' ,  G13.  6,  7X,  '  CMD5  -  ' ,  G13.  6,  4X,  '  CXD5  -  ',G13.6> 
WRITE (*,384)  CZD6, CMD6, CXD6 

384  FORMAT  (2X, 'CZD6  -  ’ ,  G1 3.  6,  7X,  *  CMD6  -  » ,  G13.  6,  4X,  '  CXD6  -  ’,G13.6) 
WRITE(*,385)  CZD7, CMD7, CXD7 

385  FORMAT  (2X,  ’  CZD7  -  « ,  G13.  6,  7X,  »CMD7  *  ' ,  G13.  6,  4X,  •  CXD7  *  ’.G13.6) 
WRITE (*,386)  CZD8, CMOS, CXD8 

386  FORMAT  (2X,  «CZD8  -  ' ,  G13.  6,  7X,  ’  CMD8  =  ' ,  G13.  6,  4X,  ’  CXD8  -  ',G13.6> 
WRITE (*,  310) 

310  FORMAT (IX, '****************************************************’  > 

315  CONTINUE 

WRITE (*, 320) 

320  FORMAT (IX, ’ IS  THE  ENTERED  DATA  CORRECT  ?  (YES/NO)') 

READ (*,330)  DATA1 
330  FORMAT (A3) 

IF (DATA1  .EQ.  ' NO  ’ )  GO  TO  686 
IF (DATA1  .EQ.  ’YES’)  GO  TO  340 
GO  TO  315 

START  THE  CALCULATIONS  TO  BUILD  LONG.  STATE  SPACE  MODEL 

340  CONTINUE 

WRITE (*, 420) 

420  FORMAT (IX,' •*••••••****•***•**••*******•••••****•*••••••*••*•**’  ) 

Z1  -  (0*S*32.2)/W 
A  »  C/(2. 0*U> 

THETA  »  DTHETA/DPR 

ZA  »  Z 1 *CZA*DPR 
ZH  =  ( Z 1 /U) *CZH 
ZQ  =  Z1*A*CZQ*DPR 
ZU  =  2. * ( Z 1 /U) *CZU 
ZD  1  =  Z 1 *CZD1*DPR 


n  n  n 


c 

c 


ZDS  =  Z1*CZD2*DPR 
ZD3  «  Z1*CZD3#DPR 
ZD4  =  Z1*CZD4*DPR 
ZDS  «  Z1*CZD5*DPR 
ZDS  =  Z 1*CZD6*DPR 
ZD7  =  Z1#CZD7*DPR 
ZDS  =  Z 1*CZD8*DPR 

XA  =  Z1*CXA*DPR 
XH  ■=  (Z1/U)*CXH 
XQ  =  Z1*A*CXQ*DPR 
XU  =  S. *(Z1/U)*CXU 
XD1  =  Z 1*CXD1*DPR 
XD2  -  Z1*CXD2*DPR 
XD3  =  Z 1*CXD3*DPR 
XD4  -  Z1*CXD4#DPR 
XDS  =  Z 1*CXD5*DPR 
XCS  =»  Z1*CXD6*DPR 
XD7  -  Z1*CXD7*0PR 
XDS  -  Z 1 *CXDfl*DPR 

Ml  »  (Q*S*C> /BIYV 

MA  -  Ml*CMA*DPfl 
MH  «  (Ml /U) *CMH 
MQ  -  M1*A*CMQ*DPR 
MU  -  2. * (Ml/U) *CMU 
MD1  *  M1*CMD1*DPR 
MD2  -  M1*CMD2#DPR 
MD3  -  M1*CMD3*DPR 
MD4  -  M1*CMD4*DPR 
MD5  =  M1*CMD5#DPR 
MD6  -  M1*CMD6*DPR 
MD7  -  M1*CMD7*DPR 
MDS  *  M1*CMD8#DPR 


WRITE  THE  DERIVATIVES 


WRITE (*, 700) 

700  FORMAT  (5X, ‘LONGITUDINAL  AXIS  DIMENSIONAL  DERIVATIVES’) 
WRITE (*, 70S) 

70S  FORMAT  ( 15X, • BODY  AX  I S  (1/RAD)’) 

WRITE  <«,  710)  ZA,MA,  XA 

710  FORMAT  (4X,  ’  ZA  =  ’ ,  G 13.  6,  9X,  ’  MA  =  ’ ,  G 1 3.  6,  6X,  ’  XA  =  ’,G13.6) 
WRITE  (#,720)  ZQ,  MQ,  XQ 

720  FORMAT  (4X,  ’  ZQ  «  •  ,G13.  6,  9X,  ’  MQ  =  ’  ,  G 1 3.  6,  6X,  '  XQ  =  ',G13.6> 
WRITE  (*,730)  ZH,  MH,  XH 

730  FORMAT  (4X,  ’  ZH  =  ' ,  G13.  6,  SX,  '  MH  =  ’ ,  G13.  6,  6X,  *  XH  »  ’,G13.6) 
WRITE(*,  740)  ZU,  MU,  XU 

740  FORMAT  (4X,  ’  ZU  «  ’ ,  G13.  6,  9X,  *  MU  =  ’ ,  G 1 3.  6,  6X,  ’  XU  «  ’.G13.6) 
WRITE (*, 750)  ZD  1 , MD1 , XD 1 


750 

FORMAT  (3X,  ’  ZD1  =  ' ,  G 1 3.  6,  8X,  ’  MD1 
WRITE (*, 760)  ZD2, MD2, XD2 

S 

’  ,  G13.  6,  SX, 

’  XD1 

’  ,  G13.  6) 

760 

FORMAT  (3X,  ’  ZD2  >  '  ,  G 13.  6,  8X,  ’  MD2 
WRITE(*, 770)  ZD3, MD3.XD3 

= 

’ ,  G13.  6,  SX, 

'  XDS 

= 

’  ,  G13.  6) 

770 

FORMAT  (3X,  ’  ZD3  =>  ’  ,  G 13.  6,  8X,  ’  MD3 
WRITE  (*,780)  ZD4,  MD4,  XD4 

= 

’  ,  613.  6,  SX, 

’  XD 3 

= 

’  ,  G13. 6) 

780 

FORMAT  (3X,  ’  ZD4  =  ’ ,  G13.  6,  8X,  ’  MD4 
WRITE ( *, 790 )  ZDS, MDS, XD5 

£ 

’  ,  G13.  6,  SX, 

’  XD4 

’  ,  G13. 6) 

790 

FORMAT  (3X,  ’  ZDS  =  ’  ,  G13.  6,  8X,  ’  MDS 

= 

’  ,  G13.  6,  SX, 

’  XDS 

= 

1  ,  G 1 3.  G ) 

WHITE ( *,  8013)  ZDS,  MD6,  XD6 

800  FORMAT  <3X,  ’  ZD6  =  '  ,  G13.  6,  8X,  '  MD6  =  '  ,  G13.  6,  5X,  ’  XD6  =  ’,613.6) 
WRITE (*, 810)  ZD7, MD7, XD7 

810  FORMAT  (3X,  ’  ZD7  =  ’  ,  G13.  6,  8X,  ’  MD7  =  ',613.6,  5X,  '  XD7  =  ’,G13.6> 
WRITE!*,  820)  ZDS,  MD8,  XD8 

820  FORMAT  (3X,  '  ZDS  =  ’  ,  G13.  6,  8X,  ’  MDS  -  ’ ,  G13.  6,  5X,  ’  XD8  =  *,G13.6> 
WRITE!*,  830) 

830  FORMAT  < 1 X, ' ****************************************************'  ) 
:  DEVELOPMENT  OF  STATE  MATRICIES 

:  DEVELOPMENT  OF  THE  PLANT  MATRIX  -  A 

VT=U 

AMAT  < 1 , 1)  =  XU 

AMAT (1,2)  =  -VT  *SAL 

AMAT (1,3)  =  XA 

AMAT  (1,4)  *  -2  2.  2*CTH 

AMAT (2, 1)  =  MU 

AMAT (2, 2)  »  MQ 

AMAT (2, 3)  »  MA 

AMAT  (2,  4)  =  0.0 

AMAT (3,1)  -  ZU/VT 

AMAT (3, 2)  «  CAL 

AMAT <3, 3)  -  ZA/VT 

AMAT (3, 4)  =  -32. 2*STH/VT 

AMAT  (4,  1)  -  0.0 

AMAT  <4,  2)  =  1.0 

AMAT  (4,  3)  =  0.0 

AMAT  (4,  4)  *  0.0 

OK,  LET’S  WRITE  THIS  SUCKER  OUT 

WRITE (*, *) 

WRITE (*, 850) 

850  FORMAT (’ 1’ , 5X, ’ LONGITUDNAL  STATE  MATRIX (BODY  AXIS)’) 

WRITE (*, *) 

WRITE (*,  842) 

842  FORMAT ( ’ 0’ , 2X, ’ FOR  STATE1 =U, STATE2=Q, STATE3=ALPHA,  STATE4=THETA’  ) 
WRITE!*,  *) 

DO  855  1=1,4 

WRITE!*,  860)  <AMAT(I,  J>,  J-1,4) 

855  CONTINUE 

860  FORMAT  <’0’  ,2X,  4  < G1 3.  6,  4X)  ) 

WRITE ( *,  *) 

NOW  WE’LL  GET  THE  INPUT  MATRIX  -  B 

BMAT ( 1 , 1)  »  XD1 
BMAT (1,2)  =  XD2 
BMAT (1,3)  =  XD3 
BMAT (1,4)  =  XD4 
BMAT (1,5)  =  XD5 
BMAT (1,6)  =  XD6 
BMAT (1,7)  =  XD7 
BMAT  (1,8)  =  XD8 
BMAT  <2, 1 )  =  MD1 
BMAT (2, 2)  =  MD2 
BMAT <2, 3)  =  MD3 
BMAT (2, 4)  =  MD4 
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END  ' 

6MAT  <2, 5)  -  MD3 
BMAT !2, 6)  -  HD6 
BMAT!2, 7)  -  MD7 
BMAT! 2, 8)  »  MD8 
BMAT  <3, 1 >  -  ZD1/VT 
BMAT  <3, 2)  -  ZD2/VT 
BMAT  <3, 3)  -  ZD3/VT 
BMAT (3, 4)  -  ZD4/VT 
BMAT! 3,  5)  -  ZD5/VT 
BMAT  <3, 6)  -  ZD6/VT 
BMAT 13, 7)  -  ZD7/VT 
BMAT (3,8)  ■  ZD8/VT 
DO  865  1-1,8 
BMAT (4,  I)  -  0.0 
865  CONTINUE 
C 

C  PRINT  OUT  THE  LONG  INPUT  MATRIX 
C 

WRITE!*,*) 

WRITE (*,870) 

870  FORMAT!  » 0* , 5X, • LONGITUDNAL  INPUT  MATRIX’) 

WRITE!*,  *) 

WRITE!*  868) 

868  FORMAT I2X, ’FOR  DELI-CANARD, DEL2-STAB, DEL3-TEF, DEL4-DR  AILERON’) 
WRITE!*  869) 

869  FORMAT I2X,  ’  DEL5-RT  RV,  DEL6-RB  RV,  DEL7-LT  RV,  DEL8-LB  RV’ ) 
WRITE!*,*) 

WRITE!*,  *) 

WRITE!*, 871) 

871  FORMAT  !•  0’ ,  5X,  ’  ROW1  ’ ,  1 1 X,  ’  R0W2* ,  1 1 X,  •  R0W3’  ,  1 1 X,  ’  R0W4’  > 
WRITE!*,*) 

DO  872  1*1,8 

WRITE  !*,  880)  IBMATIJ,  I),  J-l,4> 

872  CONTINUE 
WRITE!*,*) 

875  CONTINUE 

WRITE!*, 873) 

873  FORMAT! IX,’  DO  YOU  WANT  STAB  AXIS  DATA  FOR  LONG? ! Y/N)  ’  ) 

READ! *,874)  STAB1 

874  FORMAT !A1> 

IF  !  STABl  .EO.  ’Y*  )  GO  TO  877 
IF  <  STABl  .EQ.  ’N’  >  GO  TO  857 
GO  TO  875 
877  CONTINUE 
C 

c *********************************************************** 

c*  * 

C*  CONVERT  BODY  AXIS  DATA  TO  STABILITY  AXIS  * 

C*  IFOR  CHECK  WITH  MCAIR  DATA)  * 

C*  * 

c*  * 

c »*••»••*••*••*••••••*••*•••••*••••*•••••**••••**••••••*•••* 

c 

SMU  - !  MU*CAL  +  1MA/U) *SAL*CAL) 


SMH  »(  (  SMU  /  MU  )  *  MH  ) 

SMA  -(  MA  *  COSSQ  -  MU  *  U  *  SAL  )  >  ' 


SXU=XU*COSSQ*  <  ZA/U) *SINSQ*CAL+ ( <  XA/U)  *CAL*ZU> *SAL*CAL 
SXH  »  (SXU/XU) *XH 

SXA  «  XA*CAL**3  -U*ZU*SINSQ  -  !U*XU  -  ZA*CAL )  *CAL*SAL 
SXQ  *<  XQ*CAL  ♦  ZO*SAL  ) 

SXD1  -  ( XD1 *CflL  ♦  ZD1*SAL> 

SXD2  »<  XD2*CAL  ■*  ZD2*SAL> 

SXD3  -  !XD3*CAL  ♦  ZD3*SAL> 

SXD4  =  <XD4*CAL  ♦  ZD 4*SAL) 

SXDS  »  !XD3*CAL  ♦  ZD5*SAL> 

SXDS  «  <  XD6*CAL  ■*  ZD6*SAL> 

SXD7  *  !XD7*CAL  ♦  ZD7*SAL> 

SXDS  -  ! XD8*CAL  ♦  ZD8*SAL) 

SZU*ZU*COSSQ- <  XA/U) *SINSQ*CAL- (XU- ( ZA/U) *CAL> *SAL*CAL 
SZH  -  (SZU/ZU)  •  ZH 

SZA»ZA*CAL»*3  *  U*XU*SINSO  -  !U*ZU  *  XA*CAL) *CAL*SAL 
SZQ  -  ( ZQ*CAL  -  XQ*SAL) 

SZD1  -  IZD1*CAL  -  XD1*SAL> 

SZD2  -  <ZD2*CAL  -  XD2*SAL) 

SZD3  -  <ZD3*CAL  -  XD3*SAL> 

SZD4  -  <ZD4*CAL  -  XD4*SAL) 

SZDS  -  !ZD5*CAL  -  XD5*SAL) 

SZD6  -  <  ZD6*CAL  -  XD6*SAL) 

SZD7  »  !ZD7*CAL  -  XD7*SAL> 

SZD8  -  !ZD8*CAL  -  XD8*SAL> 

WRITE  !*,  701 ) 

701  FORMAT  !’  0’  ,  SX,  ’ LONGITUDINAL  AXIS  DIMENSIONAL  DERIVATIVES* > 
WRITE!*, 702) 

702  FORMAT  !15X,’  STABILITY  AXIS  < 1/RAD)  •> 

WRITE!*,  711)  SZA,  SMA,  SXA 

711  FORMAT  !4X,  *  Ztt  -  * ,  G13.  6,  SX,  »  MA  -  • ,  S13.  6,  6X,  •  XA  «  ’.G13.6) 
WRITE!*, 721)  SZQ, SMQ, SXO 

721  FORMAT  !4X,  *  ZQ  -  ’  ,  613.  6,  9X,  •  MO  -  ’  ,  G1 3.  6,  6X,  •  XQ  -  *,G13.6> 
WRITE!*,  731)  SZH,  SMH,  SXH 

731  FORMAT  !4X,  *  ZH  -  *  ,  G 1 3.  6,  SX,  ’  MH  -  • ,  61 3.  S,  6X,  •  XH  -  ’,613.6) 
WRITE!*, 741 )  SZU,SMU,SXU 

741  FORMAT  !4X,  ’  ZU  *  ’ ,  G1 3.  6,  9X,  ’  MU  -  ’ ,  G13.  6,  6X,  ’  XU  «  *,G13.6) 
WRITE!*, 731)  SZD1, SMD1, SXD1 

751  FORMAT  (3X,  ’  ZD1  -  ’  ,  G13.  6,  8X,  ’  MD1  «  * ,  G13.  6,  SX,  •  XD1  »  ’,613.6) 
WRITE!*, 761)  SZD2, SMD2, SXD2 

761  FORMAT  <3X,  ’  ZD2  «  ’  ,  G13.  6,  8X,  ’  MD2  «  ’  ,  G13.  6,  5X,  ’  XD2  -  ’,G13.6> 
WRITE!*, 771)  SZD3, SMD3.SXD3 

771  FORMAT  !3X,  •  ZD3  «  ’  ,  G13.  6,  8X,  ’  MD3  *  ’ ,  G13.  6,  SX,  ’  XD3  -  ’.G13.6) 
WRITE!*, 781)  SZD4, SMD4, SXD4 

781  FORMAT  <3X,  ’  ZD4  *  '  ,  G13.  6,  8X,  ’  MD4  -  ’  ,  G13.  6,  SX,  ’  XD4  -  ’,G13.6> 
WRITE!*, 791)  SZDS, SMD5, SXD5 

791  FORMAT  <3X,  ’  ZD5  «  » ,  G13.  6,  8X,  '  MD5  »  ’ ,  G 1 3.  6,  SX,  ’  XD5  »  ’,G13.6) 
WRITE!*, 800)  SZD6, SMD6, SXD6 

801  FORMAT  I3X,  ’  ZD6  -  '  ,  G 13.  6,  SX,  ’  MD6  «  ’ ,  G 1 3.  6,  SX,  ’  XD6  »  ’.G13.6) 
WRITE!*,  811)  SZD7,  SMD7,  SXD7 

811  FORMAT  !3X,  ’  ZD7  *  ’  ,  613.  6,  8X,  '  MD7  *  • ,  G13.  6,  SX,  ’  XD7  -  ’.G13.6) 
WRITE!*, 820)  SZDS, SMD8, SXDS 

821  FORMAT  !3X,  •  ZD8  =  ’  ,  G13.  6,  BX,  *  MD8  «  ’  ,  613.  6,  5X,  ’  XD8  *  ’.G13.6) 
WRITE!*, 830) 


880  FORMAT  <2X,  4(G13.  6,  2X)> 

C 

C  CALCULATE  'A*  MATRIX  ELEMENTS 

C 


AMAT (1,1) 

m 

sxu 

AMAT (1,2) 

m 

0.0 

AMAT (1,3) 

m 

SXA 

AMAT (1, 4) 

m 

-32. 2*CTH 

AMAT (2, 1) 

m 

SMU 

AMAT (2, 2) 

m 

SMQ 

AMAT  (2,  3) 

m 

SMA 

AMAT (2, 4) 

m 

0.0 

AMAT  (3,  1) 

m 

szu/u 

AMAT  (3,  2) 

m 

1.0 

AMATO,  3) 

m 

SZA/U 

AMAT  (3,  4) 

m 

-32. 2*STH/U 

AMAT  (4,  1) 

m 

0.0 

AMAT (4, 2) 

m 

1.0  - 

AMAT  (4,  3) 

m 

0.0 

AMAT (4, 4) 

m 

0.0 

WRITE!*, 851) 

C851  FORMAT  !’0',3X, ’LONGITUDNAL  STATE  MATRIX  (STAB  AXIS)*) 

C  WRITE (*, *) 

C  WRITE (*, 842) 

C  WRITE!*, *> 

C  DO  836  1-1,4 

C  WRITE (*,860)  (AMAT(I,  J),  J-1,4) 

C836  CONTINUE 

857  CONTINUE 

IF  (KEY  .Ed.  ’BOT*  )  GO  TO  446 

IF  (KEY  .EQ.  ’ GAM'  )  GO  TO  1463 

421  CONTINUE 

WRITE (*, 430) 

430  FORMAT ( 1 X, ’ IS  ANOTHER  PROGRAM  RUN  DESIRED  ?  (YES/NO) * > 
READ!*. 440)  RUN 
440  FORMAT (A3) 

WRITE!*, 443) 

445  FORMAT ( 1 X,  •  ******************************************** 
IF (RUN  .EQ.  'NO  •)  GO  TO  450 

IF (RUN  .EQ.  ’ YES’ )  GO  TO  103 
GO  TO  421 

446  CONTINUE 
C 

C  THIS  IS  WHERE  THE  LATERAL  DIRECTIONAL  STARTS 
C 

WRITE!*, 1110) 

1110  FORMAT (lX.’CLB  (i/DEG)  -  *  ) 

READ<*,*)  CLB 
WRITE!*, 1120) 

1120  FORMAT ( 1 X, * CNB  (1/DEG)  -  •> 

READ!*,*)  CNB 
WRITE!*, 1130) 

1130  FORMAT! IX, ’CYB  (1/DEG)  «  ’) 

READ!*,*)  CYB 
WRITE!*, 1140) 

1140  FORMAT ( 1 X, • CLP  (1/DEG)  -  •) 

READ!*,*)  CLP 
WRITE!*, 1 150) 

1130  FORMAT ( IX, ’ CNP  (1/DEG)  -  ’> 

READ!*,*)  CNP 


WRITE (*, 1160) 

1160  FORMAT ( 1 X, *CYP  (1/DEG)  -  • 
READ  <*,  *)  CYP 
WRITE<*, 1170) 

1170  FORMAT (1X,*CLR  (1/DEG)  -  • 
READ!*,*)  CLR 
WRITE  (#,  1180) 

1180  FORMAT (IX,  » CNR  (1/DEG)  -  * 
READ (*, *)  CNR 
WRITE(*,  1190) 

1190  FORMAT ( 1 X,  »CYR  (1/DEG)  -  » 
READ (•, *)  CYR 
WRITE  (*,  1200) 

1200  FORMAT ( 1 X, *CLD1  ( 1 / DEG )  - 
READ (*, *)  CLD1 
WRITE  (*,  1210) 

1210  FORMAT (IX, »CND1  (1/DEG)  - 
READ  (*, *)  CND1 
WRITE (*, 1220) 

1220  FORMAT ( 1 X, *C YD 1  (1/DEG)  - 
READ (*, •)  CYD1 
WRITE)*, 1230) 

1230  FORMAT ( IX, •  CLD2  (1/DEG)  - 
READ)*, *)  CLD2 
WRITE)*, 1240) 

1240  FORMAT (IX,  ’ CND2  (1/DEG)  - 
READ)*,*)  CND2 
WRITE)*,  1230) 

1250  FORMAT ( 1 X,  >CY02  (1/DEG)  - 
READ)*,*)  CYD2 
WRITE (*, 1260) 

1260  FORMAT  (IX,  *0.03  (1/DEG)  * 
READ)*,*)  CLD3 
WRITE)*,  1270) 

1270  FORMAT (IX, *CND3  (1/DEG)  - 
READ (*, *)  CND3 
WRITE (*, 1280) 

1280  FORMAT ( 1 X,  •  CYD3  (1/DEG)  - 
READ (*, *)  CYD3 
WRITE!*, 1290) 

1290  FORMAT (1X,»CL04  (1/DEG)  - 
READ!*,  *>  CLD4 
WRITE <*, 1300) 

1300  FORMAT ( 1 X, ’CND4  (1/DEG)  - 
READ  <*, *)  CND4 
WRITE(*,  1310) 

1310  FORMAT (IX, ’  CYD4  (1/DEG)  »  ' 
READ)*,*)  CYD4 
WRITE!*, 1320) 

1320  FORMAT  (IX,*  CUDS  (1/DEG)  ** 
READ!*,*)  CLD5 
WRITE!*,  1330) 

1330  FORMAT ( 1 X, *CND5  (1/DEG)  • 
READ!*,*)  CND5 
WRITE!*,  1340) 

1340  FORMAT ( 1 X, *CYD5  (1/DEG)  *  ’ 
READ ( *,  *)  CYD5 
WRITE!*, 1350) 

1350  FORMAT (1X,*CLD6  (1/DEG)  »  * 
READ!*,*)  CLD6 


WRITE (*, 1360) 

1360  FORMAT ( 1 X, *  CND6  (1/DEG)  -  •  ) 

READ (*, *>  CND6 
WRITE (*, 1370) 

1370  FORMAT ( IX, •  CYD6  (1/DEG)  -  ’> 

READ (*i *)  CYD6 
WRITE (*, 1380) 

1380  FORMAT (IX, * CLD7  (1/DEG)  -  ') 

READ (*, *)  CLD7 
WRITE(*, 1390) 

1390  FORMAT (1X,’CND7  (1/DEG)  -  •) 

READ  <*, •)  CND7 
WRITE ( *, 1400) 

1400  FORMAT (IX, ’CYD7  (1/DEG)  •  *  > 

READ(*,*>  CYD7 
WRITE (*, 1410) 

1410  FORMAT (IX,* CLDS  (1/DEG)  -  ’) 

READ(*, *)  CLDS 
WRITE (*, 1420) 

1420  FORMAT (1X,*CND8  (1/DEG)  -  •> 

READ  <*, *)  CND8 
WRITE!*,  1430) 

1430  FORMAT ( IX, *  CYD8  (1/DEG)  -  ’  > 

READ!*, •)  CYD8 
WRITE (*, 1440) 

1440  FORMAT ( 1 X, ’CLD9  (1/DEG)  -  •> 

READ!*,*)  CLD9 
WRITE!*, 14S0) 

1450  FORMAT (IX,’ CN09  (1/DEG)  -  ') 

READ (*, *>  CND9 
WRITE (*, 1460) 

1460  FORMAT (1X,’CYD9  (1/DEG)  -  •> 

READ!*,*)  CYD9 

1465  CONTINUE 

WRITE (*, 1470) 

1470  FORMAT ( ' 1 ’ , BX , ’ LAT-D I R  BODY  AXIS  COEFFICIENTS*) 

IF (KEY  .EQ.  ’LON*)  GO  TO  1490 
IF (KEY  .EQ.  • BOT* )  GO  TO  1490 
WRITE(*, 1480)  DALPHA 

1480  FORMAT (15X, ’ALPHA  «  *,G13.6) 

1490  CONTINUE 

WRITE (*, 1500)  CLB, CNB, CYB 

1500  FORMAT  (3X,  'CLB  *  ’ ,  G1 3. 6,  8X,  *  CNB  -  ’ ,  G13.  6,  5X,  ’  CYB  -  *,G13.6) 
WRITE  (*,  1510)  CLP,  CNP,  CYP 

1510  FORMAT  <3X,’ CLP  *  ’ ,  G1 3.  6, 8X,  ’  CNP  -  ’ ,  G13.  6,  5X,  ’  CYP  »  *,G13.6) 
WRITE  (*,  1520)  CLR,  CNR,  CYR 

1520  FORMAT  <3X,  ’CLR  *  ’ ,  G13.  6,  8X,  ’  CNR  -  » ,  G13.  6,  5X,  ’  CYR  «  *,G13.6) 
WRITE!*,  1530)  CLD1,  CND1,  CYD1 

1530  FORMAT  (2X,’CLD1  -  ’ ,  G13.  6,  7X, ’CND1  -  ’  ,  G13.  6,  4X,  •  CYD1  -  *,G13.6) 
WRITE!*, 1540)  CLD2, CND2, CYD2 

1540  FORMAT  (2X,’CLD2  -  ’  ,  G13.  6,  7X,  ’  CND2  -  ’ ,  G13.  6,  4X,  ’  CYD2  -  *,G13.6) 
WRITE!*, 1550)  CLD3, CND3, CYD3 

1550  FORMAT  (2X,’CLD3  -  ’ ,  G13.  6,  7X, ’CND3  *  ’ ,  G13.  6,  4X, ’CYD3  »  ’.G13.6) 
WRITE!*, 1560)  CLD4.CND4, CYD4 

1560  FORMAT  <2X,*CLD4  -  ’  ,  G1 3.  6,  7X,  '  CND4  -  ’  ,  G1 3.  6,  4X,  ’  CYD4  »  *,G13.6) 
WRITE!*, 1570)  CLD5.CND5, CYD5 

1570  FORMAT  !2X,  '  CLDS  -  ’ ,  G13.  6,  7X,  '  CND5  -  ’ ,  G13.  6,  4X,  ’  CYD5  -  *,G13.6) 
WRITE!*, 1580)  CLD6, CND6, CYD6 

1530  FORMAT  (2X,'CLD6  »  ’  ,  G1 3.  6,  7X,  ’  CND6  »  ’ ,  G13.  6,  4X,  ’  CYD6  -  *,G13.6) 
WRITE!*,  1590)  CLD7,  CND7,  CYD7 


1590 

1600 

1610 

1620 

1625 

1630 

1640 

1645 

1646 

1647 

1648 


C 


C 


FORMAT  <2X,  •  CLD7  -  ’  ,  G13.  6,  7X,  ’  CND7  -  ' ,  613.  6,  4X,  '  CVD7  -  ',G13.6> 
WRITE (*, 1600)  CLD8, CND8, CYD8 

FORMAT  <2X,  •  CLD8  -  ’  ,  G13.  6,  7X,  '  CND8  -  ’ ,  G13.  6,  4X,  ’  CY08  -  ’,613.6) 
WRITE (*, 1610)  CLD9, CND9, CYD9 

FORMAT  !2X,  *  CLD9  -  ’  ,  G13.  6,  7X,  ’  CND9  -  • ,  613.  6,  4X,  •  CYD9  -  ',613.6) 
WRITE!*,  *> 

WRITE!*, 1620) 

FORMAT ! 1 X, ' ***************************************************** 

CONTINUE 
WRITE!*, 1630) 

FORMAT !1X, ' IS  THE  ENTERED  DATA  CORRECT  ?  1YES/NO)') 

READ!*, 1640)  DATA2 
FORMAT! A3) 

IF  !  DATA2  . EQ.  'NO')  60  TO  446 
IF  !  DATA2  . EQ.  'YES'  )  60  TO  1645 
GO  TO  1625 
CONTINUE 
WRITE!*, 1646) 

FORMAT  < 1 X, • DO  YOU  WANT  STAB  AXIS  DATA  FOR  LAT-DIR?  IY/N>’> 
READ!*, 1647)  STAB2 
FORMAT (At) 

IF  I  STAB2  .EQ.  ’N'>  GO  TO  1801 

IF  <  STAB2  .EQ.  •  Y' )  60  TO  1648 

GO  TO  1645 

CONTINUE 

BSOL  PH— ALPHA 

CSA*COS IBSALPH) 

SSA«S I N  <  BS AL  PH) 

CS*CSA*CSA 

SS*SSA*SSA 


SCLP-CLP*CS  *  CNR*SS  -  (CLR  ♦  CNP)*CSA*SSA 
SCLR«CLR*CS  -  CNP*SS  ♦  !CLP  -  CNR)*CSA*SSA 
SCLB-CLB*CSA  -  CNB*SSA 
SCLD1-CLD1*CSA  -  CND1*SSA 
SCLD2-CLD2*CSA  -  CND2*SSA 
SCLD3-CLD3*CSA  -  CND3*SSA 
SCLD4-CLD4*CSA  -  CND4*SSA 
SCLD5-CLD5*CSA  -  CND5*SSA 
SCLD6»CLD6*CSA  -  CND6*SSA 
SCLD7-CLD7*CSA  -  CND7*SSA 
SCLD8-CLD8*CSA  -  CND8*SSA 
SCLD9*CLD9*CSA  -  CND9*SSA 


SCNP=CNP*CS  -  CLR*SS  ♦  (CLP  -  CNR)*CSA*S3A 
SCNR»CNR*CS  ♦  CL°*SS  *  (CL’  *  Cf.P)  *C3A*3Im 
SCNB»CNB*CSA  +  CLB*SSA 
SCND 1 “CND 1 *CS A  +  CLD 1 *SS A 
SCND2»CND2*CSA  ♦  CLD2*SSA 
SCND3-CND3*CSA  ♦  CLD3*SSA 
SCND4*CND4*CSA  ♦  CLD4*SSA 
SCND5«CND5*CSA  ♦  CLD5*SSA 
SCND6»CND6*CSA  ♦  CLD6*SSA 
SCND7-CND7*CSA  ♦  CLD7*S=A 
SCND8»CND8*CSA  ♦  CLD3*S3A 
SCND9=CND3*CSA  +  CLD3*SSA 


SCYP=CYP*CSA  -  CYR*SSA 
SCYR=CYR*CSA  +  CYP*S3A 
SCYB=CYB 


F-17 


WRITE!*, 1471) 

1471  FORMAT (SX, *  LAT-DIR  STAB  AXIS  COEFFICIENTS’) 

WRITE!*, 1501)  SCLB, SCNB, 5CV9 

1501  FORMAT  (3X,  ’  CLB  «  ’ ,  G1 3.  6,  8X,  •  CNB  =  ’ ,  Gl 3.  6,  5X,  ’  CYB  -  ’.G13.6) 
WRITE!*, 1511)  SCLP, SCNP.SCYP 

1511  FORMAT  !3X, ’CLP  -  ’  ,  G1 3.  6,  SX,  ’  CNP  =  ’ ,  G13.  6,  5X,  ’  CYP  =  ’,G13.6> 
WRITE!*, 1521)  SCLR, SCNR,  SCYR 

1521  FORMAT <3X, ’CLR  »  ’ , G13. 6, SX, ' CNR  =  ’ , G23. 6, 5X, ’ CYR  »  ’,313.6) 
WRITE!*, 1531)  SCLD1, SCND1, CYD1 

1531  FORMAT !2X, ’CLDl  -  ' , G13. 6, 7X, • CND1  =  ’ , G13. 6, 4X, * CYD1  *  ’,G13.6> 
WRITE!*, 1541)  SCLD2, SCND2, CYD2 

1541  FORMAT <2X, ’CLD2  «  ’ , G13. 6, 7X, ' CND2  »  ’ , G13. 6, 4X, ’ CYD2  -  ’,G13.6> 
WRITE!*, 1551)  SCLD3, SCND3, CY03 

1551  FORMAT  !2X,  •  CLD3  -  ' ,  G 13.  6,  7X,  '  CND3  *  •  ,  G13.  6,  4X,  *  CYD3  «  *,G13.6> 
WRITE!*, 1561)  SCLD4, SCND4, CYD4 

1561  FORMAT  !2X,  ’  CLD4  «  ’  ,  G13.  6,  7X,  ’  CND4  -  ’  ,  G 1 3.  6,  4X,  •  CYD4  »  ’,Gl3.  6) 
WRITE!*, 1571)  SCLD5, SCND5, CYD5 

1571  FORMAT  (SX,  ’CLD5  -  ’ ,  GI3.  6,  7X,  '  CND5  *  * ,  GI3.  6,  4X,  *CYD5  •  ',G13.6> 
WRITE !*, 1561 )  SCLD6, SCND6, CYD6 

1581  FORMAT  <2X,  •  CL06  -  '  ,  G13.  6,  7X,  '  CND6  ■  ' ,  G13.  6,  4X,  ’  CYD6  *  ’,G13.  6) 
WRITE !*, 1591 )  SCLD7, SCND7, CYD7 

1591  FORMAT  <2X,’CLD7  -  ' ,  G13.  6,  7X,  '  CND7  -  • ,  G13.  6,  4X,  ’  CYD7  «  ’  ,G13.  6) 
WRITE!*, 1601)  SCLDS, SCND8, CYD8 

1601  FORMAT!  2  X,’CLD8  -  ’ ,  G13.  6,  7X,  •  CND8  ■  ’  ,  G13.  6,  4X, '  CYD8  -  »,G13.  6) 
WRITE!*,  1611)  SCLD9, SCND9,  CYD9 

1611  FORMAT  !2X,  ’  CLD9  -  '  ,  G13.  6,  7X,  '  CND9  -  ’ ,  Gl3.  6,  4X,  •  CYD9  »  ’,613.6) 
WRITE!*,*) 

C 

SIXX=BIXX*COSSQ  ♦  BIZZ*SINSQ  -  BI XZ*SIN !2*ALPHA) 

SIYY-BIYY 

sizz=bizz*cossq  ♦  bixx*sinso  ♦  bixz*sin!2*alpha> 

SIXZ«BIXZ*COS!2*ALPHA)  ♦  .5*!BIXX  -  BIZZ) *SIN(2*ALPHA) 

C 

SN  -  DPR*  <Q*S*B) /SIZZ 
SL  -  DPR*  <Q*S*B) /SIXX 
SB  -  B/(2. 0*U) 

SY  -  DPR* !Q*S*32. 2) /W 
SNB  »  SN*SCNB 
SNP  -  SN*SB*SCNP 
SNR  *■  SN*SB*SCNR 
SND1  -  SN*SCND1 
SND2  »  SN*SCND2 
SND3  »  SN*SCND3 
SND4  «  SN*SCND4 
SND5  =  SN*SCND5 
SND6  =  SN*SCND6 
SND7  *  SN*SCND7 
SND8  »  SN*SCND3 
SND9  »  SN*SCND3 
C 

SLB  »  SL*SCLB 
SLP  ■  SL*SB*SCLP 
SLR  -  SL*SB*SClR 
SLD1  «  SL*SCLD1 
SLD2  «  SL*SCL32 
SLD3  =  SL*SC_D3 
SLD4  *  SL*SCw3‘* 

SLD5  =  SL*SCL35 
SLDb  =  SL*SCL36 
SLD7  *  SL*SCL37 


SLD8  »  SL*SCLD8 
SLD9  ■  SL*SCLD9 

SYB  «  SY*SCYB 
SYR  »  SY*SB*SCYR 
SYP  »  SY*SB*SCYP 
SYD1  *  SY*SCYD1 
SYD2  *  SY*SCYD2 
SYD3  «  SY*SCYD3 
SYD4  *  SY*SCYD4 
SYD5  ■  SY*SCYD5 
SYD6  *  SY*SCYD6 
SYD7  ■  SY*SCYD7 
SYD8  *  SY*SCYD8 
SYD9  *  SY*SCYD9 

WRITE (*, 1661) 

1661  FORMAT !5X, ’LAT-DIR  STAB  AXIS  DIMENSIONAL  DESIVATI VES ( 1 /RAD)  ’  ) 
WRITE (*, 1671)  SNB, SLB, SYB 

1671  FORMAT !4X,»NB  -  •  ,  G1 3.  6,  9X,  '  LB  =*  • ,  G13.  6,  EX,  •  YB  «  *  ,B13.  6) 

WRITE (*,  1681)  SNP,  SLP.SYP 

1681  FORMAT  !4X, '  NP  -  * ,  G1 3.  6,  9X,  *  LP  *  '  ,  G13.  6,  EX,  *  YP  «  »,G13.6> 
WRITE!*,  1691)  SNR,  SLR,  SYR 

1691  FORMAT  !4X,  *  NR  -  •  ,  G13.  6,  9X,  '  LR  -  * ,  G13.  6,  EX,  ’  YR  »  ’,Gl3.6) 
WRITE!*, 1701)  SND1, SLD1, SYD1 

1701  FORMAT !3X, ’ND1  ■  • ,  G13.  6,  8X,  »LD1  -  '  ,  G13.  6,  4X,  '  YD1  -  »,Gl3.6) 
WRITE!*, 1711)  SND2,  SLD2, SYDS 

1711  FORMAT  !3X,  •  ND2  -  •  ,  G13.  6,  8X,  •  LD2  «  • ,  G13.  6,  4X,  •  YDS  «  »,G13.  6) 
WRITE!*, 1721)  SND3, SLD3, SYD3 

1721  FORMAT  !3X,  •  ND3  -  • ,  G13.  6,  8X,  »LD3  -  • ,  G13.  6,  4X,  •  YD3  -  *,G13.6) 
WRITE!*,  1731)  SND4,  SLD4,  SYD4 

1731  FORMAT !3X,*ND4  -  *  ,  G13.  6,  8X,  *  LD4  -  • ,  G13.  6,  4X,  *  YD4  »  ’.G13.6) 
WRITE!*, 1741)  SNDS, SLD5, SYDS 

1741  FORMAT  !3X,  •  ND5  -  '  ,  G13.  6,  BX,  '  LDE  -  '  ,  G13.  6,  4X,  *  YD5  *  *,G13.6> 
WRITE (*, 1751 )  SND6, SLD6, SYD6 

1751  FORMAT  <3X,  »  ND6  *  •  ,  G13.  6,  8X,  '  LD6  -  * ,  G13.  6,  4X,  *  YDS  -  *,G13.6> 
WRITE!*,  1761)  SND7, SLD7,  SYD7 

1761  FORMAT  !3X,*ND7  -  » ,  G13.  6,  8X,  •  LD7  -  '  ,  G13.  6,  4X,  ’  YD7  ■  *,G13.6> 
WRITE!*, 1771)  SND8, SLD8, SYDS 

1771  FORMAT  !3X,  •  ND8  -  ’  ,  G13.  6,  8X,  ’  LD8  «  * ,  G13.  6,  4X,  *  YDS  *  »,G13.6> 
WRITE!*, 1781 )  SND9, SLD9, SYD9 

1781  FORMAT  !3X,'ND9  *  '  ,  G13.  6,  8X,  ’  LD9  *  ’ ,  G13.  6,  4X,  ’  YD9  »  ’.G13.6) 
WRITE!*, 1650) 

1650  FORMAT !  1  X,  ’ ****************************************************’  ) 

1601  CONTINUE 

N  *  DPR*  <Q*S*B> /BI ZZ 
L  »  DPR* !Q*S*B) /B I X X 
BB  *  B/I2. 0*U) 

Y  -  DPR*  <Q*S*3£. 2) /W 
BNB  »  N*CNB 
BNP  »  N*BB*CNP 
BNR  *  N*BB*CNR 
BND1  =  N*CND1 
BND2  =  N*CND2 
BND3  =  N*CND3 
BND4  =  N*CND4 
BNDS  -  N*CND5 
BND6  =  N*CND6 


BND8  -  N*CND8 
BND9  *  N*CND9 
C 

BLB  -  L*CLB 
BLP  »  L*BB*CLP 
BLR  -  l*BB*CLS 
BLD1  ■  L*CLD1 
BLD2  -  L*CLD2 
BLD3  *  L*CLD3 
BLD4  -  L*CLD4 
BLD5  *  L*CLD5 
BLD6  «  L*CLD6 
BLD7  m  L*CLD7 
BLD8  «  L*CLD8 
BLDS  *  L*CLD9 
C 

BYB  ■  Y*CYB 
BVR  «  Y*BB*CYR 
BYP  »  Y*BB#CYP 
BYD1  -  Y*CYD1 
BYD2  -  Y*CYD2 
BYD3  -  Y*CYD3 
BYD4  -  Y*CYD4 
BYDS  -  Y*CYD3 
BYD6  -  Y*CYD6 
BYD7  -  Y*CYD7 
BYDS  ■  Y*CYD8 
BYDS  -  Y*CYD9 
C 

WRITE  <  *,  1660) 

1660  FORMAT  !5X, ’LAT-DIR  BODY  AXIS  DIMENSIONAL  DERI VATIVES ! 1/RAD)  ’  > 
WRITE!*, 1670)  BNB,BLB,8YB 

1670  FORMAT !4X,’NB  *  * ,  G13.  6,  9X,  «  LB  *  ’ ,  G13.  6,  5X,  ’  YB  *  ’,G13.6> 
WRITE!*,  1680)  BNP,  BLP,  BYP 

1680  FORMAT  !4X,  ’  NP  -  *  ,  G13.  6,  9X,  •  LP  -  * ,  G13.  6,  5X,  •  YP  *  ’  ,G13.  6) 
WRITE!*,  1690)  BNR,  BLR,  BYR 

1690  FORMAT  !4X, 'NR  »  ’  ,  G13.  6,  9X,  '  LR  »  ’  ,  G13.  6,  5X,  ’  YR  »  ’  ,G13.  6) 
WRITE!*, 1700)  BND1, BLD1, BYD1 

1700  FORMAT  <3X,  ’ND1  «  '  ,  G13.  6,  BX,  '  LD1  -  '  ,  G13.  6,  4X,  ’  YD1  »  ’,G13.6> 
WRITE!*, 1710)  ENDS, BLD2, BYDS 

1710  FORMAT !3X,’ND2  «  ’  ,  G13.  6,  8X,  '  LD2  -  * ,  G13.  6,  4X,  *  YD2  «  ’,G13.6) 
WRITE!*, 1720)  BND3, BLD3, BYD3 

1720  FORMAT  !3X,  •  ND3  -  ’  ,  G13.  6,  8X,  ’  LD3  »  ’ ,  G13.  6,  4X,  ’  YD3  =  ’,G13.6> 
WRITE!*, 1730)  BND4, BLD4, BYD4 

1730  FORMAT  <3X,  ’  ND4  -  ’  ,  G13.  6,  8X,  ’  LD4  =  ’ ,  G13.  6,  4X,  ’  YD4  *  *,G13.6) 
WRITE (*, 1740)  BND5, BLDS, BYDS 

1740  FORMAT  <3X,  ’  ND5  =  ’ ,  G 1 3.  6,  8X,  ’  LD5  ■  ’ ,  G13.  6,  4X,  ’  YDS  =  ’  ,G13.  6) 
WRITE!*, 17S0)  BND6, BLD6, BYD6 

1750  FORMAT  <3X,  ’  ND6  -  '  ,  G1 3.  6,  8X,  *  LD6  »  ’  ,  G13.  6,  4X,  '  YD6  =  ’,Gl3.6) 
WRITE!*,  1760)  BND7,  BLD7,  BYD7 

1760  FORMAT  !3X,  ’  ND7  -  '  ,  G13.  6,  SX,  *  LD7  *  '  ,  G 1 3.  6,  4X,  ’  YD7  =  ’,313.6) 
WRITE!*, 1770)  BND8, BLDS, BYDS 

1770  FORMAT  !3X,  ’  NDB  «  '  ,  G13.  6,  8X,  '  LDS  »  ’ ,  G13.  6,  4X,  ’  YDS  *  ’,313.6) 

WRI TE  ( *,  1 780)  BND9,  BLDS,  BYDS 

1780  FORMAT  (3X,  ’  ND9  *  ’  ,  G13.  6,  8X,  ’  LDS  =  ’ ,  G13.  6,  4X,  ’  YDS  =  *,G13.6) 
WRITE!*, 1 7S0) 

1 7S0  FORMAT (IX,' *********■*************************************»*****'  ) 

WRITE!*, 1800) 

1800  FORMAT ( IX, • ****•*****•*****************************************' ) 
C 


CONVERSION  OF  DATA  INTO  STATE  SPACE  FORM 


'.  A  A.  1  a ivi'vuii - 


C 

c 
c 
c 

D  -  1.0  -  ( <BIXZ*BIXZ) / <BIXX*BIZZ> ) 


R1  •  BIXZ/BIZZ 

R2  »  BIXZ/BIXX 

PBNB  » 

<BNB  + 

R1*BLB)/D 

PBNP  • 

( 

(BNP  ♦ 

R1*BLP> /□ 

PBNR  » 

(BNR  ♦ 

R 1 *BLR) /D 

PBND1 

(BND1 

4 

R1*BLD1)/D 

PBND2 

(BND2 

♦ 

R1*BLD2) /D 

PBND3 

(BND3 

4 

R1*BLD3) /D 

PBND4 

(BND4 

4 

R 1*BLDA) /d 

PBND5 

(BND5 

4 

R1*BLD5)/D 

PBND6 

(BND6 

4 

R1*BLD61 /D 

PBND7 

(BND7 

4 

R1«BLD7> /D 

PBND8 

(BND8 

4 

R1*BLD8)/D 

PBND9 

(BND9 

4 

R1*BLD91/D 

PBLB  • 

(BLB  ♦ 

R2*BNB) /D 

PBLP  » 

<BLP  ♦ 

R2*BNP)/D 

PBLR  « 

(BLR  ♦ 

R2*BNR) /D 

PBLD1 

(BLD1 

4 

R2*BND1)/D 

PBLD2 

(BLD2 

4 

R2*BND2)/D 

PBLD3 

(BLD3 

4 

R2*BND3)/D 

PBLD4 

(BLD4 

4 

R2*BND4) /D 

PBLD5 

(BLD5 

4 

R2*BNDS)/D 

PBLD6 

(BLD6 

4 

R2*BND6) /D 

PBLD7 

(BLD7 

4 

R2*BND7)/D 

PBLD8 

(BLD8 

4 

R2*BND8) /D 

PBLD9 

(BLD9 

4 

R2*BND9)/D 

PBYB  - 

'  BYB/U 

PBYP  » 

'  SAL 

PBYR  - 

-CAL 

PBYPHI 

»  32. 2*CTH/U 

PBYD1 

BYD1/U 

PBYD2 

BYD2/U 

PBYD3 

BYD3/U 

PBYD4 

BYD4/U 

PBYD5 

BYD5/U 

PBYD6 

BYD6/U 

PBYD7 

BYD7/U 

PBYD8 

= 

BYD8/U 

PBYD9 

BYD9/U 

C 

C  LATERAL  DIRECTIONAL  STATE  MATRIX 

C 

C 


DO  1805  1=1,5 
DO  1806  J=t, 5 
1806  DIRMAT ( I, J) =0.  0 

1805  CONTINUE 

DIRMAT(1, 31  =  1.0 
DIRMAT  <£,  1 1 =PBYPHI 
DIRMAT ( 3, 2 1 =PBYB 
DIRMAT (2, 3) =P&YP 
D I  RMAT  < i,  <.)  =PBYR 


DIRMAT  <2, 51 =32. 2»57H/U 


F-21 


/.  -.Y.  . 

J-V-V.'. 


-----  £1 


DIRMAT(3,2)=PBLB 
DIRMAT  <3, 3) =PBlP 
DIRMAT  (3,  4)  =PBlR 
DIRMflT  (4,  2)  =P8N8 
DIRMflT  (4,  3)=PBNP 
DIRMflT (4, 4) =PBNR 
DIRMflT  <5,4)=1,0 


OUTPUT  THE  STATE  MATRIX 


WRITE  <*, 830) 

WRITE  <  *, 1810) 

1810  FORMAT <’ 1’ ,2X, ’LATERAL  DIRECTIONAL  STATE  MATRIX’) 
WRITE!*, 1820) 

1820  FORMAT!’ 0’,5X,’ STATES  -  PHI, BETA, P,  R,  PS I’  ) 

-  WRITE!*,*) 

WRITE!*, 1825)  (DIRMAT (1, I), 1=1, 5) 

WRITE!*, 1825)  !DIRMAT!2, I) , 1=1, 5) 

WRITE!*, 1825)  !DI RMAT (3, I), 1=1,5) 

WRITE!*,  1825)  (DIRMflT !4,  I),  1  =  1,5) 

WRITE!*,  1825)  (DIRMflT !5,  I) ,  1  =  1, 5) 

WRITE!*,*) 

1825  FORMAT  (’  0’  ,  2X,  5  (G1 1.  4,  4X)  ) 

LATERAL  DIRECTIONAL  INPUT  MATRIX 


DO  1830  1  =  1,  <3 
DIRBMAT  (1,  I)»0.  0 
DIRBMAT (5, I)«0. 0 
1830  CONTINUE 

DIRBMAT (2, 1)*PBYD1 
DIRBMAT (2,  2) =PBYD2 
D I RBMAT (2,3) =PBYD3 
DIRBMAT (2, 4>=PBYD4 
DIRBMAT (2, 5) =PBYD5 
DIRBMAT (2, 6)=PBY06 
DIRBMAT <2,  7)=PBYD7 
DIRBMAT (2, 8) *PBYD8 
DIRBMAT (2, 9) =PBYD9 
DIRBMAT (3, 1 ) =PBND 1 
D I RBMAT  <  3, 2) »PBND2 
DIRBMAT (3, 3) =PBND3 
DIRBMAT (3, 4>=PBND4 
DIRBMAT (3, 5) “PBND5 
DIRBMAT (3, 6) =PBND6 
DIRBMAT (3,  7>=PBND7 
DIRBMAT (3, 8) =PBND8 
DIRBMAT (3, 9>=PBND3 
DIRBMAT (4, 1>=PBlD1 
DIRBMAT  <4, 2) =PELC2 
D I RBMAT ( 4, 3) =PBLD3 
DIRBMAT (4, 4>=PBLD4 
D I RBMAT (4, 5) =PBLD5 
DIRBMAT (4, 6)=PBLD6 
D I RBMAT (4, 7) =PELD7 
D I RBMAT (4, 8) =PBLDS 
D I RBMAT (4,3)=PBLD3 


non 


?.  t 


PRINT  OUT  THE  INPUT  MATRIX 
WRITE (*, 1850) 

1850  FORMAT (’ 0’ ,2X, 'LATERAL  DIRECTIONAL  INPUT  MATRIX’) 

WRITE  <*, I860) 

1860  FORMAT  <’ 0’ ,  4X, 'FOR  INPUTS:  DEL  1 -RUDDER, DEL2-D IFF  CAN’) 

WRITE (*, 1870) 

1870  FORMAT (  6X, ’ DEL3=DIFF  STAB,  DEL4-DIFF  AIL,  DEL5=DIFF  TEF’ ) 
WRITE!*, 1880) 

1880  FORMAT !6X, ’ DEL6  TO  9  ARE  REVERSER  UANE  PORTS’) 

WRITE!*, 1890) 

1890  FORMAT!’  0'  ,5X,  ’  ROW1’  ,  11X,  ’  R0U2’  ,  UX,  '  R0W3’  ,  1 1 X,  ’  R0W4’  ,  11X,  ’  R0W5’  ) 
DO  1900  1-1,9 

WRITE!*,  1825)  !DIRBMAT  <  J,  I),  J-1,5) 

1900  CONTINUE 
GO  TO  421 
450  CONTINUE 
END 


k 


L 


Appendix  G 


To  simulate  the  response  of  the  STOL  aircraft  the  system 
in  Figure  V.  3-2  was  constructed  using  the  CAD  package 
MATRIX  XC253.  MATRIX  X  allows  the  designer  to  construct  the 
system  several  ways,  including  state— space  and  transfer 
functions.  Transfer  function  representation  for  simulation  is 
a  "natural "  since  the  required  transfer  functions  already 
exist  from  the  QFT  design  method.  Using  the  "System  Build" 
option  the  individual  transfer  functions  can  be  entered  as 
"blocks"  of  a  larger  'Super  Block".  These  Super  Blocks  can  be 
nested  as  parts  of  a  larger  Super  Block.  This  is  how  the 
STDL  simulation  is  constructed.  Figure  G-l  is  a  connection 
diagram  of  the  system.  The  highest  level  Super  Block,  ST0L15, 
has  two  inputs (commands)  and  four  outputs (two  output 
variables  and  two  equivalent'  surface  deflections)  and 
contains  six  other  Super  Blocks: 

1)  prefil  contains  both  prefilters  and  the  feedback 
loop  summing  junctions. 

2)  servo  contains  two  Super  Blocks (both  shown  as 
transfer  functions,  but  actually  are  Super  Blocks) 
servol  and  servo2  both  of  which  contain  servo 
transfer  function,  rate  and  deflection  saturations 
as  outlined  in  Super  Block  stuff  . 

3)  compen  contains  both  g  and  g  along  with  the 

1  2 

(Block  'Recof ' ) . 


reconf igurable  terms 


4)  vanesp  contains  the  plant  transfer  functions 
relating  vane  deflections  to  the  outputs. 


5)  SUPLTS  contains  the  plant  transfer  functions 
relating  aerodynamic  surface  deflection  to  the 
outputs. 

6)  out sum  contains  the  summing  junctions  that  add 
together  the  outputs  of  the  plant  transfer 
functions,  constructing  the  system  outputs. 

The  MATRIX  X  CAD  package  used  is  hosted  on  a  VAX 
computer  with  a  VMS  operating  system.  For  simulation  the  CAD 
package  takes  the  transfer  functions,  algebraic  equations, 
and  non-linear  elements (such  as  saturations)  and  builds  a 
state-space  model  to  represent  the  entire  system.  For  this 
simulation  the  state-space  vector  is  of  dimension  64. 

Needless  to  say,  simulation  is  very  slow.  In  order  to 
decrease  the  computation  time  change  the  default  integration 
routine  from  the  Variable  Step  Kutta-Merson  to  the  Implicit 
Stiff  System  Solver.  Doth  are  variable  step  methods;  however, 
the  MATRIX  X  manual  states  that  making  this  change  will 
decrease  computational  times  for  systems  containing  algebraic 
loops  and/or  dynamic  systems.  This  it  does  markedly, 
sometimes  by  an  order  of  three  or  more  depending  on  the  VAX 
load. 

MATRIX  X  significantly  reduced  the  time  spent  on  the 
simulation  part  of  this  thesis.  Previous  theses  used  CAD 
packages  that  are  not  as  powerful  as  MATRIX  X,  or  they  wrote 
their  own  simulation  routines.  In  the  authors  opinion  MATRIX 
X  is  the  "best"  control  system  design  CAD  package  at  the  time 
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of  this  writing.  AFIT  administrators  should  expedite  purchase 
of  the  CAD  package  to  host  on  AFIT  computer  resourses.  Doing 
so  will  reduce  time,  frustration,  and  stress  level  of  future 
AFIT  Controls  Sequence  students. 


APPENDIX  H 


The  loop  transmissions  developed  in  Chapter  IV  are 

overdesigned  since  they  do  not  lie  on  the  bounds  at  each 

frequency.  This  results  in  the  magnitudes  of  the  loop 

compensations  g  and  g  being  greater  than  necessary  over  the 
1  2 

entire  frequency  spectrum.  This  is  especially  noticable  at 

higher  frequencies.  In  order  to  reduce  the  overdesign  and 

associated  problems  with  wide  bandwidth  such  as  unmodelled 

pole  excitation  and  noise,  both  loop  transmissions  are 

redesigned  to  be  as  close  as  possible  to  the  bounds  and  to 

decrease  in  magnitude  faster  at  high  frequencies. 

The  redesign  of  loop  one  starts  with  the  plotting  of 

plant  templates  for  frequencies  of  200  and  400  rad/sec.  These 

extra  bounds  are  required  when  the  new  L  is  shaped  so  no 

1 

penetration  of  the  maximum  desired  M  contour (Forbi dden 

m 

Region)  occurs  in  the  range  of  80  -  1000  rad/sec  .  The  UHFB 

is  again  approximated  by  using  the  bound  at  (*)  =  1000.  The 

resulting  bounds  for  the  modified  loop  transmission  are  in 

Figure  H-l.  The  initial  form  for  L  includes  just  the 

lo 

unstable  pole 

L  =  K (0.3851)  (H-la) 

lo  (s  -  0.3851) 

For  K  >  0  the  system  must  remain  stable  and  must  drop  in 
magnitude  as  fast  as  practically  possible.  To  keep  the  loop 


H-l 


Fi a. H-l :  Loop  One  bounds  and  the  reshaped  Loop 
One  Transmission 

transmission  as  close  as  possible  to  the  -180  degree  line 
without  crossing  it  (causing  a  conditionally  stable  system), 
and  to  increase  the  rolloff,  a  pole  is  added  at  -0.3851.  To 
meet  the  bound  at  =  1  the  gain  is  increased  by  200(46  dB) . 
Using  trial  and  error  along  with  experience,  the  rest  of  the 
loop  transmission  is  shaped;  however,  this  time  the  bound  at 
CJ  =1000  is  crossed  at  a  much  lower  frequency  to  increase  the 


The  -3  dB  bandwidth  is  34  rad/sec ,  10  rad/sec  greater  than 
the  -3  dB  bandwidth  of  the  earlier  design.  The  increase  in 
bandwidth  is  due  to  the  loop  bounds  for  distrurbance 


rejection  shifting  upward.  The  required  loop  compensatii 


tain 

mifi 


IBIS*! 


miu 


l.H-3: 


and  Respective  Bounds 


60742 (s+15) (s+0.7631) (s+3.715) <s+24.31> 


(5+0.3851) (s+0. 7124±j0. 2019) vs+232.6> (s+22.49) 
*  (s+36.25) <s+0.6217ii0.2283> (s+1.979) 


(s+60. jll9) (s+1. 126) (s+5.677) (S-U0.74) 


*  <s+20. 34)  (s+33.09) 
(s+41.71) (s+220) 


Figure  H-4  shows  the  difference  between  the  old  and  new  g 

2 

Notice  the  decrease  in  high  frequency  gain,  30  dB  at  10 
rad/sec  ,  14  dB  at  100  rad/sec,  and  over  50  dB  at  1000 
rad/sec . 


Fiq.H— 4;  Comparison  of  Old  and  New  Loop  Two 
Compensat i on 

The  prefilter  f  does  not  have  to  be  changed  to  give  the 
22 

desired  tracking  performance,  thus: 

f  =  4  (H— 6) 

22  (S+0.B) 

These  compensations  and  prefilters,  leading  to  loop 
transmissions  having  much  smaller  high  frequency  gain  than 
the  previously  designed  ones,  are  used  in  the  simulations  in 


Chapter  V. 


Appendix  I 

Derivation  of  Expanded  Plant  Determinant 
Equations 


The  determinant  o-f  the  2X2  aircraft  plant  P'  is: 
detC  P>  =  [p*  P'  -  P*  P'  1  (I  — 1  > 


'>  = Ip '  P'  -  P'  P' 

1 1  22  12  21 


From  Eq(2.3-10)  (A-l)  can  be  expanded  to 


(P  +  n  P  )  (P  +p  P  )-<P  +  //  P  >  (P  +  n  P  )  (1-2) 

11  21  12  22  12  21  12  12  11  21  21  22 


Multiplying  this  expression  out  and  collecting  the  terms 
results  in: 


( 1  -  /j  fJ  MP  P  +  P  P  ) 
12  21  11  22  12  21 


(1-3) 


Ignoring  the  first  quantity  since  its  just  a  constant,  th 
terms  in  the  second  quantity  are  expanded  using  Eq(2.3-9) 


(  4  P  +  4  P  *  A  P  )<  A  P  +  A  P  > 

1  11  2  12  3  13  4  24  5  25 

-  I  A  P  +  4  p  +  4  P  )  (  4  P  +  4  P  )  (  1-4) 

1  21  2  22  3  23  4  14  5  15 

When  (A-4)  is  multiplied  out  and  grouped  according  to  the 

^  and  4  values  it  can  be  expressed  as: 

4  5 


4 

(P 

P  - 

P 

P 

>+  4 

4 

(P  P  - 

F 

P 

4 

11  24 

14  21 

1 

5 

25  11 

15  21 

4 

(P 

P  - 

P 

P 

)+  4 

4 

(P  P  - 

P 

P 

4 

12  24 

14  22 

2 

5 

12  25 

15  22 

4 

(P 

P  - 

P 

P 

)+  4 

4 

(P  P  - 

P 

P 

4 

13  24 

14  23 

3 

5 

13  25 

15  23 
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Colorado  is  used  to  design  the  control  laws  for  a  Short  Take  Off  and  Landing  (STOL) 
aircraft.  Compensators  are  presented  for  two  longitudinal  variables,  angle  of  attack 
and  forward  velocity,  which  are  controlled  via  the  use  of  five  separate  control  sur¬ 
faces:  canard,  stabilator,  ailerons,  upper  and  lower  thrust  reversing  vanes.  The 
Jfinal  design  must  exhibit  robust  qualities  over  three  flight  conditions  despite 
-'■jrface  failures, 
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The  state-space  matrix  representation  of  the  aircraft  is  developed  from  pertubation  • 
equations  using  linearized  aerodynamic  data.  Transfer  functions  relating  servo  input 
signals  to  aircraft  outputs  are  obtained  from  the  state-space  equations.  The  original 
output  set  included  the  flight  path  angle  and  velocity;  however,  the  non-minimum  phase 
characteristics  of  the  flight  path  angle  precluded  its  use  by  the  type  of  Quantitative 
Feedback  Theory  used  in  this  thesis  since  unstable  plants  can  arise.  'Instead,  the 
minimum  phase  variables  angle  of  attack  and  velocity  are  controlled.  The  ten  separate 
transfer  functions  relating  the  two  output  variables  to  the  five  input  commands  form  a 
5X2  plant  transfer  function  matrix.  These  separate  transfer  functions  are  combined 
using  a  weighting  vector  into  a  2  X  2  minimum  phase  plant  matrix  for  each  flight  con- 
cition/failure  combination.  Quantitative  Feedback  Theory  is  applied  to  the  resulting 
plants  to  yield  robust  control. 

A  single  set  of  fi::ed  compensators  and  prefilters  are  designed  to  handle  the  entire 
riant  set,  consisting  or  three  single-surface  failures  and  two  dual-surface  failures  at 
each  flight  condition.  lor  these  failures  neither  Fault  Detection/Identification,  nor 
scheduled  compensation,  is  required.  Surfaces  are  assumed  locked  at  zero  degrees  deflec¬ 
tion  after  failure,  generating  no  net  moment  after  failure.  Digital  simulations  have 
shown  the  control  to  be  robust  over  the  three  flight  conditions  and  surface  failures. 

Loop  banawidths  for  the  velocity  and  angle  of  attack  loops  are  35  and  12  rad /sec 
respectively.  Control  surface  rates  and  deflections  are  shown  to .saturate  only  for  the 
double  failure  cases. 

Quantitative  Feedback  Theory  effectively  controls  the  aircraft  despite  large  un¬ 
certainty  due  to  flight  condition  changes  and/or  control  surface  failures  without  identi¬ 
fication.  Application  of  QFT  eliminates  the  use  of  identification  to  achieve  robustness  . 
and  the  associated  false  alarm  and  missed  detection  problems.  Efforts  to  expand  upon  ^ 
the  base  of  flight  control  design  using  this  method  are  recommended,  especially  direct 
design  in  the  discrete  domain.  Research  should  also  continue  on  developing  a  computer- 
aided  design  program  to  expedite  the  synthesis  of  controllers  using  QFT. 
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